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ANNOUNCEMENT | | 
| INTERNATIONAL CONGRESS OF MATHEMATICIANS, 1954. 


First Communication. 


The International Congress of Math&maticiuns 1954 will be held in Amsterdam 
from September 2nd to September 9th under the auspices of “Het Wiskundig 
Genootschap” (The Mathematical Society of the Netherlands). It is the sincere hope of 
the ‘‘Wiskundig Genootsehap” that the Congress 1954, which will be open to all 
mathematicians from all parts of the world, will be a fertile international gathering. 

The Organizing Committee has invited a number of outstanding mathematicians to 
deliver one-hour addresses, hopiog that in this way a survey or the recent development 
in the whole field of mathematics may be furnished. 

There will be seven sections. viz.: 
1. Algebra and Theory of Numbers. 2. Analysis. 
3. Geometry and Topology. 4. Probability and Statistics. 
5. Mathematical Physics and Applied Mathematics. 6. Logic and Foundations. 
7. Philosophy, History and Education. 

Those who wish io attend the Congress are requested to communicate their name 
(with degrees, qualifications etc.) and full address to the secretariat as soon as possible. 
They will receive a more detailed communication m due course. 


Amsterdam, 
Od Boerhasvestrast 40. 


Price per volume Rs, 12/- (India) and Rs. 13/- (Foreign) 


The Organizing Committee. 


ON A GENERALISATION OF BATEMAN’S k-FUNCTION 


By 
N. K. CTAKRABARTY, Krishnagar, West Bengal 


(Recered— January 10, 1952) 


. 

1. In a paper communicated for publication, Mitra and Srivastava have obtained a 
generalisation of Bateman's k-funelion. A similar generalisation has also been obtained 
by Chakrabarty (1952) in a recent paper published in the Proceedings of the Academy 
of Sciences, the Netherlands, He has introduced the notation 4! (x) which represents 
the generalised k-function of Bateman, viz., 


Ae 
Loi = 2 J 2! cos'0 cos (z tan «—n0) d0, 17-—1, 
T 
0 


In the two papers mentioned above, various properlies of the &k'-functions have 
been obtained and certain relations with the other known special functions have been 
deduced. In the present paper, further avtempts have been made in obtuining a number 
of properties of this function. 

2. Let and m be positive integers (including zero) and z > m. 


Since f(p) = h(z), we have 
d p e zy 
—— . R — qd. - 5 
p{(p) = iz h(x); provided A(0) = 0, an FESTE pus (1) 
SEH Die 4 1)p"51, d" 


(pa Dar dee OM) "om (2) 


Let us write m =n —l—1 and z = n and interpret with the help of 


edebi (şa) ei nc P s (5) 
” Bə (L+ p)#-1+1 z š 
We get E 9] (—1)*-!-1 d n—1—1 TT 
a key (3a) “nar Ep (a ERTL (4) 
In Jacobi’s ıdentity 
I'(n4-1-4«) «(1— t8. F ( — +1+4a4+ +v: t dd i*t.(1—1y) T8 
T +a) t (1 ) 2 i n, n a, p, ) ) dö" { ( ) ie 


let us writen =a, n4-14- +Ñ b and 14 = c and replace tby s/b. Then making 


b—co, we get the relation 


F(a; c; x)= en at) erraten), (5) 


Oo comparing the right hand sides of (4) and (5) we obtain easily 


a iz) MAUS +1, 2? +2; 22) (6) 
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` and 


(22) "e" karya) = (Hal 


—. + : ( 
Pu per z)”, [b 


By using a result due to Howell (1987), viz 


ü Ç Ta- B+) Bn ; 
La (z) = cs rl Pi = B) Laena) 


we get, after slight reductions, $ 


Din +27-+2) < o ri D(Q7-B+1) 


Putting 8 = + in sucession, we can deduce the ns 


eI thy ir) - 5070000005 00 Lie (8) 


Qd ağ. LU er, pi) mp 
2(n-Fa-1) 3x) = QHT(924 > Dr (n +27. +2). rina)! h)  Ə “257 ) 
and 
z? Thee, ze ren 2 m (Certa D öz +, 


whence we obtain the integrals 


® kətana (127) Diss: d: (123)D £2, de _ (—1)*7?51 (22m + 1)! 12: +n—m+t#) (Gei Pa < n) 
en O maA —m)lnl Dia +3) In + 922 +): — 


mm 0, m > T : (10) 
TF), = Arm! m) Ts tn mo (am) 


ES qii de Fn —m)| m ET (2x + 8) D (n +22 + 2) (msn) 
=O (m). (11) 


8. Wecan write cu in the form 


p(l-py) ` 
pepe - LY (: m. (1 (+ pyrtrti- Se) 


and then interpreting term by term with the help of 


“1-1 
də: 5. 


nal] 
Qo tgn — > ` a 12 
Tin + 1j Wee < | + ss aa bp (z). , ( ) 


Writing (12) in the form 
-1-1 


HM 
Qng 72v, Za > 7 23581 Š 
————. 6” lote 14:1) (z ) 
y 
I'(n-4 1) < 


and interpreting both sides with the the help of 


we obtain from (1) 





Mem amannan aan Me 





a This connection of Laguerre polynomial with the Generalized Bateman's Function leads us to surmise 
that the kl function can be utilised in the study of wave mevhanical treatment of the (heory of hydrogen atom, 
This connection has also been obtained by Mittra and Srivastava. 


$ Obtained by Srivastava independently, 
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COT .. Deut. s, (4p) (Howell, 1987). 


and oe Kd) s (— 1) DR + s ho Por ap?) HE W Ls copia pages Gp?) 
( Chakrabarty, 1952, Srivastava ) 
We derive, after slight simplifications 


n—- 2-1 
Ont) n+ LIT (90 3) Wr -2.- N) 
D 3, 2. 1 2y1r+f 
p D 34 (3p) Jalan Tj —” y (3p?) x 


¢ 
x W air-aan, -r-i Cp?) (18) 
Again, we have the relation 





des and the well-kuown formula (Mclachlan and Humbert, 1041) that if f(p) hf), 
a 


then he - DI am Së Lu. 


We obtain 
gom zs ici (E. ) T. (E a a 0. (14) 


Let us next write the relation (Chakrabarty) 


KU — (1 + lest DEEL) [t] < 1 
n 0 
in the form 
° 2... 
Ima nz?) = (1+ t)4'¢ +i” 
In wm () s 


and interpret both sides by known operationul representations. 


We obtain after alight reductions 
(L+t)j%ttaltieiet Hy , _ (1 -+ ta} 


€ Tili” S (HTW 44 2-4, porter, (2) ` (15) 
nm w Ü) 
where the left hand side can be expressed in terms of the “Error Function” by the rela- 
tion erlelz)=4x "tet" W Lu +, (z2) 


4. Let us now write 
21+1p(1 — pyis-d-l z olaq (1 — p)ià-1a-1 (1— pyir-de-i 
Sed. p)is++l+1 7 p — (1 + p) tat! pe (1+ p)prdeti 


where A+r=n and a+s+2=!, we obtain, ou interpretation with the help of the 
product theorem, 


Hz) = 2 | Bay) AC) dy (16) 
The following particular cases can be e at onco 


kit (ar) = d Kay) ky (z) dy 
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9 
kita) = 2 | e s) dy 
f 0 
| Mu (s) = 2 f e rta) dy, C1 (17) 
S 0 
Again we have ei 
_ 21 Qab X—1)— H r 
e~it ka,(šz)=(—1)#"! Ta pp 
and B 
Therefore writing pk-m+4 1 ph-l pratt 
(1 Epis D ° (1 + peti ° (p +1) ater! 
where n+ u = k and m =l+v+1 and interpreting with the help of tho prcduct- Aheorem 
we obtain 
e em ANA iz) = (—1)9-!-1T' (2m +1) 


m-1-8/3%2 (ie — Mis l (18 
Pam- T) Pia (4e — 8) x Mi-s s -i-i() dy 118) 
Putting k = n, | = 0 and m 1-3 and observing that 
ne) = (7 1971120 + 2) en (be), 
Mon) = 2T (n + 1)ekI, (hz), 
we geb f 
J ais Di Ry = CHE Kae), 120 (10) 
0 TE 
Putting k = n and m = 1+ in (18), we also get , 
J Fe) sinh dy dy = ese (20) 
0 
Again writing ` 
pr! E 1 ph- —m +f ' 
(nitt) m 


p- att 
p Epiteti (1 +p)Atatr 
where k+À =n, m y = l and proceeding similarly as before 


. we get 


J y m-31(z—y)” My ml — y n-h,l DÉI dı 
0 


| 2l—-2m+1>0 
n—-l-1 >z ' 12 21 
(—1) | I(2m + DI 2l — 9m +1) kaalt) Ed m (21) 
B. From the definition of the Generalised k-funetion, it is easy to verify that $ 


1 Yrs X " 


21—1 : 
Ziel = kü ca) — cl) 2 Q 
5 This integral has been obtained Ly Mitra and Srivastava independently 
§ This recurrence relation was also obtained by them 
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and d? 1-9 
Xə Kante) = ir) AN la). (23) 


We obtain after considerable reductions and muking use of the formula 


erh (x) = = 942 D)" prc 


pri 
z əə 
E (exa) = seki), (>) (24) 
By the Principle of Induction, we can write 
| DIN jetski (a) = of “ki “M (42). (25) 


for every positive integral value of N. 
We can extend the definition (25) to non-integral values of N *, vis. 
| Deietekäi Dell = ef hays (42) (26) 
when 8 is a positive integer, the expression on the left of (26) ) represents the atb, differen- 
tial co-efficient of the expression within the bracket. 
6. Let us, 3, ROW consider the ə. of the two absolutely convergent series 


he 


nre PL, cl: =F ien) 


Yur Yay YR 3 T YR, 22.3 = 
aP, í ; a’) = aP,( j= btaa) 
Öl, ds, 2 e " d 


The coefficient of æ% in the above product 


Dice (z,), - x” mir Joan Rad pA Ui: 


TNT (~A) I (8), ALTE) 
where : II(vr)y = (a), (25) a la) 
the above can de written a8 
sët, EE Tei, Tü P 
m. ie +1, Wd), Va 


8 I Léa ns -B.—ktl, yp 2 ” 
p! IK Ba) y -9pm pel; ër : ər 


Consequently, after some reductions**, we have 


my sələn (s ce) 


QN He» .. ~Be-p+l, yn 2) | 
ITB) Le Ge He ee ə. 


* See Proc Edin. Math. Soo, (2) 6 11 (1999) p. 75, a paper by C. G, Lambe or Ferar—bPioo, Roy, Boca 
Edin, 48 (1997-28), 92. The argument is on the same lines, 
** Tt 1g unnecessary to give the details of calculations; 
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Since ta)” 


_ _ Í -— 
J (2) mə Ton” D 12 ), 


we obtain from the above 
‚Fila; 8; E J, (bz) 


= nn (öz) 2 TT ə... IR —a —. (zz) ) (28) 


next 





Pub a = 6, weet 





ve] (b dk Ze 2 
6 (ba) = Tv +1) EST (30 Co "H: vtl; 4a? (20) 
The following particular cases are at once deducible 
2 1 
ra Tan (bz) = (242) 2 Cmm kənə gö) (80) 
" T R 
erg cos Ge — Dri LO =) E — Z JD al Ze. (31) 
y Za 
"E ; dar I(-4-u) b 
GH? ba! gin bz = 75 Hi: a'z?) eD, LZ ) 02 
1-14 nl i *) 2a-+1 V 9d ( ) 


Lastly writing 
F (a; o; $e") in the form 


P (z; p; iz?) = nn) "ett da) ate 
Best ” 


and observing that 





(4z)*** = Flu +n +1) > ə (je)? Jutntp(z), (Watson, 1948) 
p=0 
we obtain, in view of absolute convergence of the repeated series involved ( u being non- 
negative 3, the relation 


F(a; p, bü) 
23 a... > E TAG Wei 


=$ Ab ( en 1i) H (da) tr A (s): 


maf 


the series ,F, being terminating ; ER we obtain finally 


ex(ay)tr-9-1 ks, (y) 
5--:"” pti, —m 1) x (J, 2/3) - dü 
TL+ 2) 21-21 m 


which is the formula expressing Ai in the form of nn infinite series involving Bessel 


= (— 1)"-1-1 





functions, 
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My best thanks are due to Dr, 8. C. Mitra for the keen interest he has taken in the 
preparation of this papor. 
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ON THE BENDING OF AN ELASTIC PLATE--IV 


; Dy : 
H. M. SgNauPTA, Krishnagar, West Bengal. 


( Received —January 18, 1982 ) 


1. The following pages contain an expression for the deflexion of thg central plane 
of a thin elliptic plate made of isotropic-elastic material, clamped at the edge, carrying 
a load of weight W spread uniformly over a rectangular region on the upper face of the 
plate. Prior to writing this paper the author had solved the corresponding problem for a 
circular plate. But when he was about to send his results for publication he received a 
paper by Prof. Tuyosi Sekiya and Atusi Barto (1949) where they have given a formal 
expression for the deflexion of the circular plate under arbitrary distribution of load in an 
integral form and they worked out a particular case where the load is concentrated ata 
number of discrete points. ‘The author then turned his attention to the problem discussed 
in this paper. He introduces elliptic co-ordinates for the solution of the problem and 
ultimately converts the plane harmonic and binarmonic functions appearing in the solution 
into rational ıntegral functions of z and y thus preparing the results into forms easy for 
numerical calculation, Biharmonic functions in eliiptic co-ordinates appear first to have 
been introduced by Póschl (1921) for solving elastic problems with elliptic boundaries. 
Pöschl (1921) was followed by Happel (1921) and Galerkin (1928) who gave two sets of 
bibarmonic functions in elliptic coordinates suitable for some particular problems which 
they solved. An exhanstive list of biharmonic functions was given by Timpe which have 
been used by a large number of workers in the field. 

2. Let + = , (a > b > 0) be the equation of the edge line of a thin elliptic 
plate made of isotropic elastic material. Suppose it to be clamped at the edge and let it 
carry a load of weight W spread uniformiy over a rectangular region X, < z < Xa 
Y, < y < Y, contained entirely within the plate. | 


As usual, we introduce elliptie co-ordinates by the transformation 
f z+iy=ccosh(£+in (c > 0). 
Then £ = a, represents the bounding ellipse, provided a = c cosh z and b = c sinh a, 

Let £ = 8 be an ellipse confocal with the bounding ellipse and lying in its interior, 
such that the rectangle X, < zz Xa Y, < z = Ya is composed entirely of the interior 
points of the ellipse é = 8. Let (x y) be an arbitrary point in the open region X,<zg<X, 
Y, <y < Y, 

Writing R? = (z—2z,)?-(y—3y,), we may establish the following result provided K 
denotes, the rectangle X, zzz < X, Y, sy zz Y, and X(n) and Y(n) stand for the 


n VR n vn " 
numbers i-i and rar respectively. | 
d 


2—1840P- 1 


10 H. M. SENGUPTA 





J [9 a Rae, ay, = EFES tt cos wn + ol sin m] 
where fol) = ELE) + LI? exp(2é)+L® + LO ent 28) 
III = 1 + Li’ exp é+ Le exp( - 0 * Lf? exp( — 8£) 
P11) = £M.) + Ml oxp P exp( — £) + MÜ) exp( — 8£) 
f) = £L,($ + LY + L© exp( — 26) +L exp(—4£) 
e£) = MP + Mf? exp( — 2E) + Mf? exp( —4£) 


and for n > 8 
$) = Lf? expl — (n —2)£] + LO exp( —n£) +L expf — (n + 2)8] 


delt) = MO expf—(n DE + MO exp(—n£) + MO expí — (n + 2)£] 


and 





L,(é) = 2X(1)¥(1) cosh 264 [X(8)Y(1) + X(1)Y(3)] 


at t 
r (0 = LAY) 


ara cosh £ 





u =D i ç 


L4(£) = 2X(DYU) 
a) — X(1)Y(1) > b’ 
LY — log n 





LG =—* (ig FO oYryayym xey 
0 2) 5 2 [X(8) Tn (1) DU 


3. 42 - 
LƏ = a (tog FE +2)+ [X(3)¥(t) —X(1)Y(8)] 




















q? a 
2 pa 
LPD =- EE 
(B) ei 
Le = == [Ú | — 2j) + 2 — 





(8 — 2 ais r 
Li Z= = le IXY) XFN - X moi 








MY) = ees (10g EE „jora 


Dd 
Ei 


(9) 
(10) 


(11) 


(13) 
(14) 
(15) 


(16) 


1X(4)Y(1) + X(2)Y ey lam 


(18) 


(19) 
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| nP = es [@ CC? @)¥(@)+X()¥(4)}] (20) 
Mt = wë EU (8) ¥(2) — X(1)1 wh x()Yo)| (21) 
LP = (log ey (1) -xq(9] (22) 
LO m sen OD i5 128Y0)-X00)Y(8)1 (23) 
o) "gp (5) Y Q1) - OX(B) VB) + X COY GJ] (24) 
— (1) - X())Y(3)] +4X(DY (1) (25) 

MM = X(2)Y (2) 
M® = “geyə [X(4)Y(2) + (27) 
MP ges [X OY) - -— END (28) 
x — = 7 Gë | > 
2572 
MP ^u ziya" ppens A ei D 
"= ern Mna em ge 


where 


Z=X+iY, Z=X-iY P,(Z) = (2 


_n(n—4) (n — 5) (a E b?) (2Z)n-* + 


^ (qa payor? RE (qs — poe 


21 


öl 
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The last six of the above formulae are symbolic. ‘he expressions on the right are 
first to be simplified taking X and Y as algebraic symbols and in the final result X” and Y? 
are to be replaced by X(n) und Y(n) respectively for all positive integral values of m. 


(XY p| 55 || A 
In particular erden 1 | 
E X*y -8X3Y3 -x+¥| 
` V a* — b? L3(a" — b?) 
2 9 | — M 
Qe 2 [ 2 gun -8X(2) Y(3)) — X(8)) İl 
Bo Lf = rla XY) -8xevey-xevo 


A process of term by term integration has been effected above. The process is valid 
on account of the uniform convergence of the series involved in (z, yı) so lond as iz, y) 
lies in the annular space 8 < £ < and Üz, yı) lies in the closed rectangle X, < 2 = X,, 


Y,sysY,. It may be proved that e fe log R may be obtained by term by term 
K 


differention of the above series for all £ > 8 and so in particular for £ = a. 
8. Let us now consider the expression 


TV II 
EE E a 36 
i 8:DX()YQ) / 7 --—.. i 


"where (Timpe, 1923) 
w = A,’+A, cosh £ + (4,7 cosh é + A, cosh 3£) cosh q +(B,’ sinh £ + B, sinh 8£) sın y 


+ > [fán cosh (n —2)£ + Ay cosh n£ + An cosh (n +2)E} cos ny 


n=2 


+ {Banua sinh (n—2)£+ B,” sinh n£ + B, sinh(n + 2)€} sin ny] (37) 








The boundary conditions for clamped edge, viz., o = > = 0 over £ = « are formally 
satisfied if we take z 

I (f(x) . cosh 22 

A, As = Q 7 (38) 
fo (a) sinh 2a 

A, A, = Q fola) p (39) 

” f,(x) cosh 8zİ ` | 

A, A = Q (40) 
Lol 3 sinh 82 
f?) asha 

A: A, — Q (41) 
f (x) sinh o 
" | z (a) sinh 32 | 

D, Ay = (42) 


ie, (a) 8 cosh Ba 
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R d sinh e 
B, — (49) 

GR cosh a 

fala) cosh (n + 2)« | 
Ay An =—A„_,(n sinh 4a + 2 sinh 2na) + Q | (44) 
fn (e) (n+2) sinh (n +2)a 

fanla) cosh na 

Ån Ân = Ans Ar (45) 
fn (e) mn anh na 3 
S qa (a) sinh (n + 2)a 
By A, = Baal sinh 4«—2 sinh 2na` + Q (46) 
— | I q (z) (m + 2) cosh (n + 2)e 
= əz oa (a) sinh na 
B,A = Bas 24-2— Q (47) 
Pn‘ (x) n cosh na 

where /(a2—h3 

Des IE eds (48) 
SD X (DY (1) 
( "A, = sinh (2n + 2)a + (n + 1) sinh 2a (49) 
and m 

A, = sınh (2n+2)¢—(n+1) sinh 2a (60) 


Explicit formulae for An, BA, Ay’ and B,” may be easily written down. 
We may prove that the series I A, cosh n£ cos ny is absolutely and uniformly 
convergent within a closed ellipse that contains the ellipse É = a in its interior. Itis also 


P„(z) where 





possible to show that it is the real part of the series > ^ 


P,(2) = ann (ar lët HEH enger: i4 (8 Š “ 5) (a? — 3) (22) 


ue es 


Since the series may be proved to be uniformly convergent in a closed ellipse con- 
taining the ellipse £ = « in its ınterior, its sum is an analytic function throughout the 
interior of the bigger ellipse, and denoting the sum by u(z, y) + io (z, y), we have 

I u(z, y) = > A, cosh n£ cos ny 

Since u(x, y) is the real part of a function analytic in the-region mentioned above, 
it is finite and continuous throughout the region and its partial derivatives of all orders 
with rospest to x and y are also finite and continuous throughout the same region. Again 
Vu = 0 throughout the interior of the region, Without entering into compex mathe- 
matical details we shall content ourselves by asserting that the formal expression for w’ 
piven by formula (87) represents a function of z and y finite and continuous throughout the 
plate including the boundary and that it yields continuous partial derivatives with respect 
to z and y throughout the same -— and that the equation Vo” = 0 is satisfied 
throughout the interior of -ellipse 0 x $ «2«. ‘The expressions contributed by o” to the 
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das functions G, H and N are therefore finite and continuous throughout the plate 
including the boundary. It is also possible to demonstrate that the value of Aw’/äf for 
£ = a. may be had by term-by-term differentiation of the series for o” with respect to € 
and then putting £ = x in the resulting series. This justifies the formal process that bas 
been employed to get the values of the constants Az, Bn, dy’, Br’. 


4. Now writing 


we may shew that 
8U(z, y) Las y) + 1, ,(z, y) = 1,,,(z, y) — L,,,(z, y) 


+B (X, = X DG yılı Y (e@—X,)°-(@ -X p°} 


X z—X 
+ (ym Y ft -1 T 3 — tan”! d 
W d ER y—Y, i Yy-Y, 


U (z, y) = [fe log Rdw dy, 
K wer 


-X z—X 
(y — Y,)* < tan ey. y tan pre Y, 


1 


+ (a -xy Stan Sa SE — tanti Y, zi 


-(z— X, (ten Ua 5 Än — tan”! Za =- (D 

where Tim, y) = tz— Xi) - Yə - Xy + (y — id log f(z- XF Yu) (52) 
La, Y) = (z— Xə — YA — — (u— Y,)*} log f(e— X, + — . (58) 

Lg, y) = (z—X,)(e— YJ- X + (y— Y, Gi log f(z— X,)* + (y - Yl (54) 

J, (z, y) = (z—Xj)(z — Y,M(z — Y, + (y — YƏ" log f(z— X): + (y -Y;Y] (55) 


the inverse circular functions involved have their principal values. 
We now replace the formula (86) for w by 


W os 
T" , 56 
e=- Ta px Yay mY te Ge 


where U(x, y) is given above and o” is given by the formula (87) of section 8. 


I, (e, y) is not defined for (X,, Y,), the lower left-hand corner of the rectangle 
Xi gra An Y, Sys Y, Butit tends to zero as (x, y) tends to (X,, Yı). We ascribe 
to 1,,, (X,, Y,) the value zero and similarly we set I, (X,, Y d = To (Xa, Ya) = Jes 


(Xa Ya) = 0. On the other hand the-inverse function tan”! 24 


722 is not defined for 
Y, 

: ' : 4 -1 z — X; 

y = Y, Itis however bounded throughout the plate, The function (y—Y,)* tan DE 

n a 

therefore tends to zero as y Lends lo Y,. We ascribe to this function the value zero for 

y= Y, Similarly we set (z— Xə)” tan! = — equal to zero for z = X,. Similar speci- 
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fications for other terms involving inverse functions are made. This makes U(x, y) and 
80 the function o finite and continuous throughout the ellipse 0 < £ = € including the 
boundary. We may shew that derivatives with respect to z and y up to the third order 
are finite and continuous throughout the plate on our making necessary adjustments for 
the values of certain functions appearing in them for z = X,, z = X, y = Y, and y = Y, 
Thus the stresses G, H, N which depend on derivatives of not more than the third order 
with respect to z ane y are finite and continuous throughout tne plate including its 
boundary. | 


8. Coming now to the difierential equation, the proposed deflexion function w is 
expected to satisfy, we can establish that 


W 








4. 4 
Vie Dry VH 
S W |( — z - X, ra) 
an DADİ” yx Oy, 








+(ton™ = tani £71: ) 





1 z—X, 
_ —1 z—X, —ı ül. | 
(tan oy ten all (57) 


at all points of the plate excepting nt those that lio on z = X,, z = X, y= Y,, y= Y,. 
These lines divide the interior of the ellipse into nine compartments. We call them 
D,, Da D,, Dy, Di De, Dy, D, and D, 
MESE A, Y, < y < Y,] 
Diz ES , [z > Xa, y > Ya 
£— < z < X, wy Y. 
D,zBls x Xu y> Y] ` ) 
D,zE[e < X, Y, c < Y] 
Diz Ele < X, yx Y] 
ODE <a < Xn y<Y,] 
Di PE > X,, y < Y,] ` 
2... e z< X, Y, <ç < Y,] 


Evidently D; represents the interior of the rectangle X, < z < X,, Y, «y < Ya: 
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R, therefore contains in its ınterior an arbitrary piece of the line y = Y, to the right 
z= X, 


Similarly let R, contain an arbitrary piece of the liue. z = X, lying abovo y = Y, 
within it. This will bo ensured by X, < ğı) < X, < & and Y, « 1, 
Similarly lot l 
for R, ë LXK Aa Y,< h 


for H. Y, <, < Y, < ia & < X, 
for R, Ny «Y i= N; < Ya £s < X, 
for Hh Š, < X, < £, <A, Ha <Y, 
for R, X, < £, < X, < Es, Ne < Y, 
and fo Re neuen «Y, sh 

We shall further consider the rectangles R,, Rio Ri and R,, all lying entirely 
within the plate such that R, contains in its interior an arbitrary piece of the line z = X, 
lying between y= Y, and y = Y, . This will be ensured by X, < £, € X, < é, and 
Y, <n, < n, < Y, Similarly for Ra, we take X, < &, < £4 < X, and Y, <q, < Y, <, 
for By) we take ğı € X, < £, < X, and Y, < n, «y, < Y, and lastly for R,, we have 
X, < i < &, < X, and q, < Y, < 1, < Y,. 

Finally we consider four more rectangles Riy Rio Ris anad Rie alllying entirely 
within the ellipse, containing the four corners of the rectangles X, < z < X, Y,sysY, 
in their interiors respectively. We take R,, to contain the corner (X,, Yə) in its interior. 
To ensure this, we take 


for Ry, ALE LX, Lén Y, <, < Y, <4, 
for Eu, & < X, Lé <L Xn Y,<n,<Y,<"n, 
for Rs & < X, <, < X,, m, < Y, < n, < Y, 
and for Rio X, << <<, w < Y, < n, < Y, 


We now consider the value of the integral J Nds along the boundaries of the rect- 


angle R,,(m = 1, 2, 8,...... 16) where 8,, denotes the boundaries of the rectangle Fo. 


By direct integration, we can shew that 


INN 
Fee Ofk —X Y — Y) I e 
/ Nds = X Cam ps Yj (59) 
where Y, 
İc > pra+rt+s — Y,) |: n7162 Ag "an. xil (60 
(— ) Di X4) (7 r D He — Y, $ — X ) 


BO 
where each of p, q, r, 8 assumes the valuos Land 2 and where the summation extends 
over all possible permutations of 1 and 2 in the four positions p, 9, r, 8, any integer being 
permitted to be repeated any number of times. 
Now consider the rectangle R, where X, < ğı < ğa and Y, < n, < Y, « 14. 
8 -1$40P—I 
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It is easy to verify that , 
= Ae LE X4) (m~ Y- (&- A) Y)-(&—-XJ(-Y2-*(6-7X)00 7 Y2] b 


m 0 
[nae = 0 


Bo 
This proves that there is no concentrated point load or line load along the line y — Y,, 
v > X. Similarly we can show that there is no concentrated point load or a- line load 
along any of the lines y= Y, s «X,; y—Y, 22 X, y= Y, z < Ei ëss X, 
y > Y,; z = Aa y < Y,; z = X,, yo Y, and z =X, z < Y.. 


“Taking the integral Jas nlong the boundaries of the rectangle E, viz. Eet < Š, 


7); S& S 1, where X, an doy Y, <n, < n, < Y,, E 


Ww 
[ae = = a, -X F, XY-Y (X, — £ı)(65- m) 


Now; smee every point in the interior of X, < € XA, Y, « y < Y, is subject to a 
uniform pressure TV /(X, — X )(Y, — Y,) and every point of the region Y, < y < Y, z < X, 
is free from load and since £,, Ea: 11s yə are perfectly arbitrary subject only to the condi- 
tions prescribed, it follows that the line s = X,, Y, «y < Y, is free from any point load 
or a line load. In like manner, we can shew that the lines 

x= X,, Y, < < Y, 

y= Y, X, <a < X, 

y= Y, X, < € X, 
are from any concentrated point load or a line load - 


We now consider the integral il Nda taken nlong the boundaries of the rectangle R 


84 


vis. SCSEs 1 SY Sta where X, < £ < X, < éan Y, «n < Y, < a, we find 


UNE ao C 
J was = X, Y) re ?1) 


which considered with the results alrendy obtained proves that the corner (X,, Y,) is free 
from any concentrated point load. In like manner we prove that the other three corners 
via. (Xa, Yu), (Xa Ya) and (X,, Y,) of the rectangle X, < z = XA Y, «y Y, are also 
free from any concentrated point load. 


13? 


We have therefore established that the function 


V 
o EBEN U(x, 3) +w 


where U(z, y) is given by (51) and ə” is given by (37) and the constants- 2 zu Ba Di 
are given by the formulae (88), (89),...... (50) read along with the formulae (2), (8),...... (85) 
represents the deflexion of the central plane of a thin elliptic plate. clamped at the edge 
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2 2 
“til = Ü carrying & load of weight W spread uniformly over the rectangle 


X, <z < X,, Y,« y < Y, on the upper face of the plate. 


Note.—It may be mentioned here that the deflexion of the central plane of a thin 
2 2 
elliptic plate clamped at the edgo-; u = 0, curıymg a load of weight IV along u 
line parallel (o the minor axis, say ulong z = X, Y, «y < Y, may bo written down 
from the expression giving tho deflexion in this papor by making X,>A,. Tho limiting 
function obtained by effecting the abovo process, term by torm, may be shown to satisfy 
all the requirements of tbe problem, 


KRISHN\GAR COLLEGE, 
KRISHNAGAR. 
NADIA, WEST BENGAL. 
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ON THE GENERAL THEORY OF DOWNWASH BEHIND 
A FINITE WING IN SUPERSONIC RANGE 


By 
M. Z. V. KnzvwoBLOoOKI, Urbana, Illinois, U. S. A. 


(Communicated by Prof B. R. Seih— Received November 15, 1951) > 


INTRODUOTION 


The purpose of the present paper is to derive the most general theory of the 
downwash phenomenon behind a finite wing in a nonviscous isentropic fluid. The 
integro-differential equation obtained in the last step of the theory is solved by means 
of Picard’s iteration process. As is known, there are four items which may cause 
discrepancies between the numerical results based on the linearized theory aud the 
experimental data: distortion and displacement of the vorticity sheet, volume distribution 
of vorticity and influence of viscosity. The first two items jointly with the influence 
of a high angle of attack are taken mto account in the present note. The third item 
seems to be practically of small importance and the fourth one has not yet been 
determined. The author assumes that the reader is acquainted with the papers of 
‚Previous writers on the subject, and consequently, only a few items will be cited from the 
previous literature. 


Part I. Basic Notions AND EQUATIONS. 


1.1. Potential Equation and Its Transformation. 


Assume a steady, irrotational, three-dimensional flow of a nonviscous, isentropic 
fluid past a body. The velocity potential equation has the form: 


0730, Da Pik = 9, ul = e (1.1.1) 
where the ordinary tensor notation is used and the system of coordinates (z,, Y,, 2,) 
ig looated as follows. the z,-oxis is in the direction of the free-stream flow, the 
Ti. = y,-axis lies in the horizontal plane perpendicular to z-axis and the 2,, = Z-axis is 
vertical upwards. Let us decompose the velocity components in the sense that 


u — U +P, e, = Uu, zz U tu, u; = Pre, = Q Us, d = 2, 9, u; =v, us =w. From 
the energy equation c?+ BV? = ez 4.8V2, B = 4(y-1), one may calculate 


c? = cz2[14+Z,(1-Z,)"], Z, = Bez?Y,, (1.1.2) 
where 


Y, = V2— U? = 2Uu, + Siu)? = 2Ug., + Dip). P JJ S44 
Introduce the affine transformation : 


€ = T4 Y = ui 8 MEL, Hz | - (12 M23)]17, 
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Ccmbining (1.1.1) and (1.1.2) results in: | 
Ges Py Piss = [Pp  P(yiz, y, eil = Giz, y, a), (1.1.3) 
Hiel = [(1.— M2) (c5 —8Y,)] "I, G, = BMEY pm, 
Y, = 209.4 (P,e)? + KEP), G, = Deg, 19,49. 


where i, j, k denote now the coordinates 2, y, 2 and the quantities e and y are: 


: lfı=k=z, 
(1 if i =g, 
ux cT | |= y fi=z, k= y or g, if £ — y or z, k= z, 
B^ 1 1 = y or 2, 
f p* if i= y or e, k — y or z. 


` For some future purposes introduce the operator Dp = 9,4—9,y—9,;, which applied 
to a scalar is analogous to a gradient and applied to a vector is analogous to a divergence. 
Fundamental solution associated with (1.1.8) is 


. pes il f J r Gd£dgdi, 
K 


where r” = (a-2,)’-(y-y2)’-(#-2,)*, and K is the domain in (z,, Ya, 2.)-space defined 
by the inequalities a, = z, (z—z,)* > (y—9,)! --(«—2,)*. This solution represents the 
potential at P(x, y, s) of the space-(z,, Na #,)-mass distribution of the variable density G. 
In ease Hg = 0, the fundamental solution is 9 = 2-1, Of course, the variables of 
integration are the (£3, Yas 2,)-coordinates. | 


1.2. Modified Green's and Gauss's Theorems. 


If n denotes the inner normal to S, of D in a three dimensional space, then 
Green’s second formula is: 


ili [uv -vyv'u]de,dy,ds, = - f f (v. —».18,. (1.2.1) 
D _ 1 8, 
Ihe co-normal v to 5, is defined by v, — —n, ='cos (r— 2) = —cosa, v, = n, = cos f), 


v, = n, = cosy, where m i = 1 to B, are direction cosines of n. Repeating the procedure 
applied to derive Green's second formula, one easily obtains  ''modified" Green's 
formula (Heaslet and Lomax, 1948, a, b): 


| [Jo -vorulda,ay,dz, = J [Io aas. (1.2.2) 
D I g | 


Thus, the operator Y? transforms into -0° and 8/8, > 8/Əv. ` 


In case v = +-1, there is a singularity at P. On the right side ın (1.2,1) one must 
add the term pu, where p = 4r, 2r, or 0 for P inside D, on 8, or outside D, respectively, 
An analogous situation is in (1.2.2), i.e., on the rıght-hand side of (1.2.2) one has- to add 
ihe term —pu, The sign — is due to the fact that passing from (1.2.1) to (1.2.2) we 
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multiply the entire equation by —1. Let zl, or D*v = O(r-!), u = p, then (1.2.2) 


transforms into: 
pp = - f J 7. 9672,348,* f f aap, (1.2.8) 
8, D 
+ 


+ + + + + + 
Let in (1.2.1) u = 1, V*v = div grad v = div a, ora = grad v, and a = ia kia t ka, la | 


++ 
x Olio, Hence, v,, = grad,v = a, = u.n. Thon from (1.2.1) one gets Gauss's theorem: 


+> 
f J J ir Der f f aids. (1.2.4) 
D S, 


The right-hand side of (1.2.4) may be written in the form 


- f J (eos v -Fecs B + cos y)d S,. 
8, 


> 
Similarly, put in (12.2) u = 1, oO?» = m. grad v = Da, and obtain the "modified" Gauss 


theorem : 
> 
J J [paap = f fods, (1.2.5) 
D 5, 


NL + + ++ 
As is seen, diva > —-Da and a.n — a.v. Transform the second integral in (1.2.8) by 


means of the modified Gauss theorem. Evidently it must be: 
> 
[la = a,,,—a,,—a,,, = TIG = M(a, y, g; 2 Yo, 92). 


2 Img 


This is Lagrange’s linear equation whose subsidiary equations (equations of characteristics) 
are: 


döş f1 = dy,/—1 = dz,/-1= da/M. 
Applying a set of multipliers one gets : 
da = 4M(dx,—-—dy,—dz.). 


from which it is obvious (since da = a,,da,+ a,, dy, +a,z dz,) that Gan, İM, ayy, = de, 


=—4M, and hence, 
1 y E 
a= sÍ i Mdz,— d Mdy,— il Mds, , (1.2.6) 
Va 


fa 


where the integrals extend from a point (2,, Ya, 2,) on S, to the point P(z, y, 2) in 
K or D, Thus, the final form of (1.2.3) is 


po = - f [ipn 967.4 ds, (1.2.7) 
8, 


with a, = a(— cos a + cos 2 + eos y= ah,. Since the function a has a singularity at P, 


one has to investigate it in the neighborhood of P. Assume a sphere C, of radius e around 
P as the origin, 
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L-I 


Since dC, = ei sin edgd0, the last integral in (1.2.7) takes the form : 


Thus, ` 


Za: Un £5 OD Ce 27 
lim (dz,-dy,—dz,) | də İ «Gh. sin gdp 
0 


c-Ü 0 0 


i = lim 4re(Gh,h.);.. Yas Za on Oe = 0, 
c +0 


with h, = (z,—3,—2,). In all the cases (P inside X, on S, or outside K) the integration 
over C, yields zero, and there is no effect of the singularity at P (due to a,) upon the 
form of the integro-differential equation (1.2.7). In the linearized theory a, = 0. ` 


1.8. Picard’s Iteration Process. 


The solution of the integro-differential equation (1.2.7) will be achieved by 
Picard’s iferation procedure. The most general form of (1.2.7) is 


p = İ [ras, F == F(a, h., 12. 9, Prvyy ATE Prik). 
8, I 


The proof of convergence of the iteration process is the same as in (Courant and 
Hilbert, 1987) or (Tamarkiu and Feller, 1941). Consequently, only a few deviations 
from the pattern given in (Tamarkin and Feller, 1941) will be marked. The first 
part of the proof is unchanged with F, = P (z; ha r-!, fal, TÖ, o, gs Ports Porth) 
v, being an identically zero function and f(?, f, denoting certain given initial function 
of p or p,,, which should be used only once at the beginning of the procedure. Instead 
of (8.40) in (Tamarkin and Feller, 1941) one has the relations: 


Pn = f İP.-as, Puy = Pn = f [7.-..28, (1.8.1) 
8, 8, 


Pri = In SUE Qs, = in = f Parad, 
. 8, I 8, 


Assume that there exists n number K > 0 such that for all æ;, ha, 171, 9, ete., for which 
max (|e, İpi, lg |, ° 0 fe), ° .) < K, Ps: ere i, 8— exist and max (HET İF. 8— 


|F, |... M. With max (İga-il, İpa-il, |n- - .) < K, one obtains: 
| Py, | <= MÜ, hz < MÜ, | qal < MÜ, |En | < MT. 
With MÜ < K, we can be sure that Pa ..., are defined. Similarly, as in (Tamarkin 
and Feller, 1941) apply the law of the mean and with 
I, =f f | Qn — Pn—-1 | + | Pa—Pn-1 E lan — In-ı |+ İsa — Ba Dd 8,, 
S 


and I = BE ye; Fən — F q q |+ F 2818, etc., obtain : 
| ai Pn | < Mİ, | Pati Pn | < MII, cer ı "In | € MPI, 
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8544 —85 |  M?P(1--MI?1,. . - . (18.8) 
Denoting by S, the maximum value the integrand of I, éan take “we obtain as the sum 
of the absolute values in (1.3.8) the expression: : — 


MPQ-SMP.MMOS, = 8. - (1.8.4) 


The last part of the proof of convergence with « « 1 and the proof of uniqueness are 
identically the same as in (Courant and Hilbert, 1937) or (Tamurkin and Feller, 1941). 
In order that the solution exists, the condition « « 1 must be preserved. This holds for 
any M if S, — P is sufficiently small, - or for any S, if M is sulficiently small. This 
latter occurs in the method of small perturbations. 


4.4. Transformation of Multiple Integrals. 


Let the transformation z = plu, v), y = Yu, v) represent a continuous one-to-one 
mapping of the closed region R of the zy-plane on a region R’ of the uv-plane. Let the 
function a = f(x, y) = x(u, v) represent a ‘surface which lies above the region E (or R”) 
of the zy (or uv)-plane. Let all the functions p, V, g, x possess continuous first derivatives 
and lət the Jacobian D, = Ir, y) /Ə(u, v) be everywhere positive, whereas J acobıans: 
D, = Id, x)/O(u, v), and D, = O(x, p)/Əlu, v) may have any values or may even vanish. 
Then, as is known, (Courant, 1987), the area of & curved surface is given in the 


(z, y)-system by the formula ` 
— Ae f f+ alt + donas, | (1.4.1) 
R 


with dS = dzdy and R representing the projection of the surface on the zy-plane, and 
in the (u, v)-system by the formula 


2 J J rei, (Dp T'548', 020. (14.3) 
.. m 
with dS’ = dudv and R” representing the projection of the surface on the uv-plane. 


A double integral over a curved surface 2, J J FdS may be represented either by means 
of (1.4.1) ss f [r^ ...]%dzdy in the (z, y) system or by means of (1.4.2) as 
R | , “r 
J [riz ..]%dudv in the (u, v)-gystem. >; MEL 
R' 


1.8. Hadamard's Improper Integral. 


The integral. I f f 
: 
1- İ ugs üz | v uo sip 

a I j ] 


4—1840P—1. 
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is meaningless (p being any integer). To have somehow a means of calculating the 
integral (1.5.1) Hadamard (1923) introduced a new kind of improper integral (the ''finite 
part” of the integral in question with the symbol[ ), Without going into too many 
details, we cite here only the final result. 


(1.2.5) 





A(z)de _ A(b) ,  —_(—1)'AG-5(b) “fa | A,(z)dz. 
(b= 2) PR  (p-po-ay- i (p-0H6-a)B^ J Gaerne 


where 

A (z) = A(zr) — 14(b) A (b)(b —z)-- ... -(—-1)*4G-P(b)(b—z)*7]. (1.5.20) 
One may introduce another representation of a finite part of an improper integral (Hada- 
mard, 1928). The indefinite integral (1.6.1) with p = 1, say, can be written in the form 
F(b)--C—F(a)—C. Develop the function A(z) in Taylor's series, 1.e., A(x) = A(b) 
t(z-—b5)A' (b)-- ..., and caiculate the first integral of the expansion; i.e., 


i Teo _ 240) Ë 
: ( 


S T n. 1.5.3) 
b-a >a] ( 





Inserting the upper limit b into the integral (1.5 8) gives a meaningless value. One can 
show that such a case occurs only for the first term of the expansion. Hence, in order 
to obtain a finite value of that integral one must subtract that expression from the 
original value of the integral. In other words, 











1, =| ° A(z)dz 
° J (b= a)? 
| a 


It follows that if C is chosen so that 


WEI 2Alb) _ | 
c= tim ii 
then the expression (1.5.4) for I, may be written 
I, = — [F(a) + C]. 


When C is chosen in this manner, the notation for the caloulation may thus be modified 
to the form 


T -Fla e 24(b) 
= lim İTF(z) +C— id — C s|. (1.5.4) 


— f* A(z)de 
> J (bay 


where the asterisk indicates that the upper limit is not substituted into the indefinite 
integral, F(z)--C, C being calculated as explained above. 


x —[F(a)+C], 





Part II. Wina Tugony AND Vorticity ŠHEET. 


0.1. Pressure. 


Assume a supersonic flow past a lifting surface of symmetrical pianform. The 
plane of symmetry lies in the z,2,-plane, with +z, directed downstream, +y, sideward 
and te, upward. From the Bernoulli’s equation, one gets 
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y 
p = p—- ba =~ J eVdV. 
Since Y, = V?-U?, VdV = 4dY,. With e = e«(1—2,)/09. one may solve the integral in 
question by elementary method obtuinmg : 
p= H(Y,) = encéy t [(1—Z,* -1], < = yGy— 1)”, (2.1.1) 


Let the subscripts u, l, denote the conditions on the upper and lower parts of the lifting 
surface, respectively, then, with q = 4o4U?. 


Ap[q = (pi- pa^ = (P-P) = pochy [A -Za -(L-Z,4)]. (2.1.2) 
From (2.1.1), one easily obtains: 
u, =—-U[1-(1-U“By?], B= o, tpa — 688711 (14-pye2 ez )'^]. — (8.1.8) 


Applying twice the binomial theorem in (2 1 3) and preserving only one term, one obtains 


p/q = —2U-lu,, a known quantity from the linearized theory. 


2.2. Supersonic Wing Theory 


Assume, for the sake of simplicity, that a supersonic wing is substituted by a flat 
plate of zero thickness with the trailing edge parallel to y,-axis. Denote the angle of attack 
by 8, and introduce the additional system of coordinates (ti, ir, ni) where t, denotes 
the tangent and n, the normal to the plate respectively. Then, 

t, = z cos ô,— s, sin ô N, = x, sind, +2, cos ÖL, 
and 
ui = u! cos Öl — 45 sin 8,, un = ni gin $, + 8 cos 6,. (3.2.1) 


With i I 
u = U cos $ +u,  u,U sin Ó, + tt, 


the components u, and u, may be easily calculated. Moreover, it is obvious thut 
V? = (u)! +07 + (uy). 


The following conditions are superimposed upon the system: the upand u, = v-velocity 
components are antisymmetric and the u,-component is symmetric about the plane of the 
lifting surface. ‘Thus, the velocity components on the top and bottom surface of the 
lifting surface are related by: 
Uy = Uy, VF — Vy, Un Un 

88 n — surface. The discontinuities in u, and v constitute a vorticity sheet. Preserving 
all the notions and definitions introduced into the supersonic wing theory by previous 
authors, we obtain the velocity potential from the integral 


-> + 
g x Hds, 


l 
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and the circulation from ER. sy 
Ite il V-ds, the integral being taken round the circuit 


where 


+ + + + 

V = iu+jv+kw, 
and -> > 

ge io, + jy + en. 


Assuming that the circulation 18 equal to the j jump in 9 when s. to the lifting 
surface, one obtains: 


eT Yir | “ir 
I = Ag, =| (us ude, + | (ou vi)dy, f (w„-wide,, (2.2.2) 
Tir M | Be 
where the subscripts L and T denote the points on the leading and Lrailing edge, 
respectively. Equation (2.2.2) can be written as 


sch 
Ag, = | grad pu- 9-42, 


To obtain the circulation per unit span, one has to integrate along the line y, = const., 
i.e., this is given by (2.2.2) without the second integral. For comparison purposes, one 
may reduce the conditions to a linearized flow: 8, ə 0°, n, öp, Wy = W; and Ag, is 
given by the first integral in (2.2.2). Preserving only one term in (2.1.3), the difference 
(ua — uj) is equal to U(2q)"! (p; — py), which is known from the linearized theory. To have 
a lift, there must be: p> 0, Pu < 0, i.e., the function H(Y,) (2.1.1) is an odd function 
of n,, since it changes its sign when approaching to the lifting surface in the n,-direction. 


2.8. Worticity Sheet. 


The vorticity system behind a lifting surface consists of a vorticity sheet and two 
vortex cores except just at the trailing edge, where it is only a vorticity sheet and no 
cores, and at infinity, where there are only two cores and no vorticity sheet. Although 
the fluid is nonviscous, let us assume a transfer of vorticity from the center part of the 
vorticity sheet towards the cores. The total strength I' of the vortex-lines trailing from 
eich half of the lifting surface remains constant at all stations. Following (Spreiter and 
Saaks, 1051) assume that the strength of the vortex-lines I’, between two cores 18 given by 

P Doz, Y) = Mes YP), OS P < 1, | 
where I'(z,, dl refer to the trailing edge. Hence, 


ə 
I' ei Dody + I. 
where Yy} and 1) denote the variable lateral position of the ‘venter and the strength of a 
core, respectively. Since ! 


Frau =T,-T, 
0 
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one has I', = T,+P(T,-T). For z, o co, the semi-span s, of completely rolled-up cores 
is given in (Spreiter and Sacks, 1951) and is equal to 1C (cc), where Cr is the lift 
coefficient of the lifting surface of area S, and (c:c)y,=o is the product of the section 
lift coefficient and the chord at y, = 0. Similarly, the function P(z,) is calculated in 
(Spreiter and Sacks, 1951) and is equal to: 


hə 
Pis) = İns 911n4e9- | OT. (2.8.1) 
0 

An approximate determination of the diameter of the core, as given in (Spreiter and 
Sacks, 1951) may be preserved in the present case with two changes: instead of o,, the 
quantity o,, and instead of a”, the quantity y, should be used, respectively, where o, 
may be assumed to refer to the center of the core. Assuming that the diameter a, of the 
core may not be much less than 8, we may calculate a, from the original equation in 

(Bpreiter and Sacks, 1951) instead of from the simplified equation, thus obtaining: 


8 (8,)* = 2 exp «[exp (22) 2 1]77, (2.8.2) 
I . a = 1610981 Cp, (0,801)? — 1/4. | 8.8.20) 


Preserving the analogue with the lifting surface, let it be assumed that behind the surface 
one has: Apm = Le or Lv. To find the shape of the vorticity sheet at any distance 
behind the trailıng edge, we have to caleulate the down-and side-displacement ofa 
vortex-line. Similarly, to know the value of P(z,) (2.8.1), we have to know the function 
Ji» which refers to the vortex-line originated at the tip and next displaced in both 
vertical and transverse directions. 


Any displacement in y,-or z-direction may be calculated by means of expressions 
; MS 
dz, = —uldt, i = 2,8, ul = uç, us, 


(since, in general, the vortex lines approach to the z,-axis with time, hence, dy, « 0). 
The differential equations of a vortex filament trailing behind a lifting surface are 


(with dt = (u])"1da,): i 


qi Ë 
_ 1 Tos : 
iz Zy = ~ | u; (ul) !dg,, 1 = 2, B, uj = us, U3, (2.3.8) 


EZ T 


= İ Pele) den i= 2,8. (2.8.4) 


Cit 


or 


‘Transforming onto (a, y, #)-or (zz, Və, 6,)-system results in: 


x, I . : 
Yor Yor =, OY 34: = Sdt GU t 9,2) "dën: = Yor a. (2.3.5) 


To find Y,., one has to insert into the left-hand side of (2.8.5) the expression y,,~8,, and 
on the right-hand side to apply the limits of integration from zz up (o z,. From equations 


dn, badi, 1 = 1 to 3, a, = Pay 
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aj dz, = —a, dy, = —u) de, 
or E d B - | 
Yı =— ja a,dz, +y (bÜ, 2, —— [ar a,dz,+2,,(t), (2.8.6) 


where the parametric representation of the initial curve (trailing edge) is given by the 
expression : | 


g, (t) = e, (z, (t), y, (1) (Courant and Hilbert, II p 55). 
The parameter t may be eliminated and the function z,, representing the vorlicity surface 
behind the trailing edge, may be given in the form z, = 2. (z,, yu). Assume 
g, (t) = X (Yi10(t)), 
then: z a 
zg = 2, =- IE adz, + X (y, + IE a,dz,), (2.8.7) 


and in (z, y, z)- or (z,, Ya 2,)-8ystem, since 2, = ës the same surface is given by: 


8 = ps, = -a| İüseəyəəda t xsl at [ (0*9 7.02]. 


Let: 
z,,( = at+d, y(t) = bt+e, göt) = otf, 
where small letters denote constants. ‘Then, 
Sa =cb (fi — e) +f, 
and the second term on the right-hand side of (2.3.8) transforms into 


ara f (U+ pre) Pdea] +i. 


When the trailing edge is parallel to (z,, y,)-plane, then c = 0, gi, = f, and equation 
simplifies enormously. An element of area of the surface z, = 2,(z,, Yi) 18 given by 
(1.4.1) with dS = dæ dya. The partial derivative z,,., is equal to 

(— DÉI, ag cbe, )(U + quts 


and 2,,,, may be obtuined by means of differentiation under the integral sign. 


Parr 111. FINAL Form or Egvarion, 
9.1. Integro-differential Equation. 


The problem of the down- and side-wash behind a finite lifting surface may be easily 
reduced to that one discussed in Part I. The same procedure is applied as in (Heaslet 
and Lomax, 1948, a, b; 1047). A point P(x, y, s) located above the lifting surface 1s 
assumed. It hes on the surface of the volume of integration which volume is out by 
the lifting surface into two parts. Hence, the surface is also cub into two parts, upper 
and lower, respectively : 


8, = Sout Sat. 
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Ba consists of the forecone I’ of P, “front surface” A (envelope of a family of Mach cones 
‘having vertices on the leading edge) and the upper area of the lifting surface za, The 
Mach cones in the (z, y, s) or (£a, Ya, 2,)-8ystem of coordinates have their axes coinciding 
with or parallel to the z,-axis and their vertex angles are 90°. ‘Thus, a conormal lies 
on A. The equation of I’ is r = (z — z,)” — (p — yə)” — (z — 22)” = 0. Since P lies on 
Sau p in (1.2.7) is equal to 2r. Since P does not lie on Ba p = 0 in that case. 
Subtraction of those two equations of the form (1.2.7) gives the result: . 


2rp = -$ | İs, - [ [5455 f = gm — gir Tu Dan (8.1.1) 

8, 7 88 I 
withi =u, l. Asis known, r represents a unit source and (r-1),, a unit doublet. The 
quantity r”! becomes infinite along F. To overcome this difficulty, Hadamard’s method 
of the finite part of an integral is applied (Hadamard, 1918). Regarding the integration 
over A, in (Heaslet and Lomax, 1948a, p. 80), it was shown that p= 9,,=0onA. Since 
uny derivative or may be represented as g,,v,; it is seen that G(1.1.8), M,a(1.2.0) and 
a,(8.1.1) vanish on A and the integration over A is excluded. Thus, only the integration 
over Tau and 74; should be performed. Since the surface of discontinuity consists of the 
lifting surface and of the vorticity sheet behind the lifting surface, the final form of 
(8.1.1) is: 


Ang = — 


Ta 





where the subscript w denotes the vorticity sheet and rs in all the four terms are given 
by the second form in (8.1.1). 


8.2. Boundary Value Problem. 


Usually, it may not be possible to satisfy the boundary conditions in an exact, but 
only in an approximate way. For this reason, ib has to be discussed seperately in each 
particular case. Two kinds of boundary value problems are the most common: (a) The 
load distribution over the lifting surface of zero thickness is supposed to be known and 
given in the form of Ag, or (u,—) (2.2.2), Briefly, the jump in the velocity potential 
over the area of the lifting surface is known. (b) The lifting surface is of a finite 
thickness and shape of the airfoil is given. The slope of a cross-section (airfoil) must be 
equal to the slope of the streamlines at the surface, 1.6., 


T £s = İs, A gələm, Be, = non Urn . (8.2.1) 


where £ = &(2, yi), or £ = &(z,. yə) denotes the function desoribing the airfoil section, 
including the angle of attack. As is known, in the linearized equation vie = 0 or 
D’P=0, the function p may be assumed to represent either acceleration or velocity 
potential or any of the three perturbation velocity components. The explanation of this 
fact is simple. Differentiating Die = 0 partially with respect to z,. y, or s, and letting 
Ga = P, we see that Dec Q does not change its form. Similarly, from Prandil’s 
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definition of the acceleration potential (Heaslet and Lomax, 1948, p. 15,) (5,4, = Vo, —. 
Gy, = V Q, ete.), we see that E = 9,2, and hence, again D’p = 0, may be used to find i e 
In the present case, it is obvious that this procedure is much more complicated, 
Fundamenlally,. equation O’p = G, refers only to the velocity potential. In case one 
would like to put g = 9,,, = u, the form of (he function G,., is not identical with that of 
G. If one performs the indicated differentiation G,,, and puts 9,2, € v, one can use the 
general method* presented in the present work. Of course, on the righthand side of (1.1.8) 
there will appear new integral forms since, for example, 9,,, must be represented ag 


il P, ye, dz, == il P, y,d2,, etc. 


Part IV. SAMPLE OF APPLICATION. 


4.1. Lifting Surface. | 

For illustrative purposes, the theory presented above will be applied to the simplest 
possible case of a lifting surface. Assume in (zı, Yı, #,)-Or (Za, Ya, 22)-aystem a flat 
lifting surface 7,(r,) of finite dimensions, zero thickness, with the leading edge coinciding 
with, and irailing edge parallel to the y,-(ys}-axıs. Let in İz İ-system the induced 
velocity component normal to the surface be a continuous function across 7,, i.e., if 
subscript 8 denotes the surface: 


Apin, = (Pu- Pdi, = Aue Dan APs y ni I Apunys = 0, (4.1.1) 


with nyu = nu = N (nii nj, nia). Moreover, vw = vi = vi(vi Via, Vis). From “the 
normal form of the equation of a plane, i.e., v, cos%,+y, cos B, +2, cosy, = Pı, with 
a, = 7/2—-8,, B, = 7/2, r, = Bach, p, = 0, and ni, = ru = sind, f, = v, = 0, 
his = Viz = 608 8,, one easily obtains the equation of the lifting surface 
g = s, = fe, ya) = — z tan ö,, Filt, Yu, 2) = zi tan öy +8, = Ü, (4.1.2) 
Hence, (4.1.1) may be written in the form 
Ads, +Â, z, tan 8, =0, 


where | 
` Ag, = A921, y, % = —2, lan SÉ) 


since the value of o, should refer to points on the lıfting surface. Since the applied 
transformation z, = 2, Ya = HÜ), Za = u2,, denotes stretching in y-and $-divections, the 
corresponding surface in the des system is again a flat one passing through the 1/,-axis 
whose equation is 


a= ps, = f,(z;, ya) = —p a, tan d,, F, = x, tan ö, +g, = 0. (4.1.8) 

By comparison, one easily obtains the relation between the angles of attack : 
tan 8, = p tan ô, sind, =pBsind,, cos 8, = B cos ô, — (4.14) 
B = [1+(p!-1) sin? 5,177, (4.1.4) 
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where the positive sign of the square root should be taken. With m, = nangi, nagi naş), 


Va = vivi, Vaz: Vay): Ba = — Vay, Ban = Vag, Tan = vas, One has 


Ə 8 Oz, Om, Ə dy, dys, O Os, Or, 


8v, Öz, Oz, Ov, By, Oy, Ov, Öz, Bz, Ow, 
which in the present case (on the surface), gives: 


a (ə) » 9 . 
a: x. ó. + — L Š o 


8v, Be, (4.1.5) 


Hence, cne may easily derive the formula 
Ae, = BAY, Bi = —(u *l) Ben, (4.1.6) 
with Ag, referring to points on the surface. Since 
17 = (x— x) (y — Ya) — (a - əy)” (2-2)? ey year), 


after a fev transformations, one has: 


(rt), = —7?Byu[(z—2,) sin 6, + (s — ə)) cos 8,]. . (4.1.7) 
The functions a,, v=u orl, may be represented by means of (1 2.6). 
x= Y £ 
a, = 7 rG dr — il T G,udy, — il SIM (4.1.8) 
T, y: f, 
G, = (cà -BY (F QT PF. rə = Ze, + 21(9,, (4.1.8a) 
FÉ, = BMEY Prazı (4 1.8b) 
H, Kom 2Up,, Prz T Dei, ak. . (4.1.8c) 


With z = g, = —z, lan, = nz, tan 0,, D, = u, D, = p tan 9,, D, = 0, one gets 


7 2 ds, = T. .. ull+tan?d,)Y?dz, dy,, 


which may be verified by applying (1.4.1) and (4.1.8) and transforming the so obtained 
result [1+p-? tan78,|"de,dy, onto the (rı, Yı, z,2space by means of (4.1.4) and 
dy, = udy,. Thus, the formula for o over the lifting surface area takes the form: 


ane = -E İ fäert sin 6, + (g — a,) cos 8,]Ap, 


+ Br Aşı, I (du T aj) de, du, (4.1.9) 
K = allen, yn = (e@—x,)*~p?[(y—y,)? + (s 2)?], (4.1.90) 


with B given by (4.1.4a), B, by (4.1.0) and a, by (4.1.8), and Ag, = Ag,(z,, Vi, 2, 
=-r,tand,). Ihe area of the integration 7, 18 the projection on (z,, y,)-plane of the 
area contained inside the curve of intersection of the plane of the lifting surface with the 
forecone I. In {z,}-space, the flat surface is given by (4.1.8) and the forecone by 


T = (z—z,)*+ ---. Eliminating s, by means of (4.1.8) results in 


| Lë ell —(y — ya)” — (z +z, tan öy)? = 0. (4.1.10) 
5—1840P-—1 


34 M. Z. V. KRZYWOBLOCKI 

Equations (4.1.8), i.e., 24 = g,(z,), and (4.1.10), i.e., y, = y,(z,), represent the curve 
of intersection in the space {z,}. The projection of that curve on the (z,, z,)-plane 
(y; = 0) is a straight line, and ils projection on the (z,, y,)-plane (g, = 0) is a hyperbola 
(4.1.10). In {x,}-space that curve may be found in an analogous way or by means of a 
direct transformation from {z,}-onto {x,}-space and by use of (4.1.4). 


4.9. Yorticity Sheet. 


Assume thut the vorticity sheet rj; or Tz, extending from the trailing edge of the 
lifting surface {air} or {zar} up to infinity is a curved surface of zero thickness, that” 
Ap = 0 and the induced velocity component normal to the vorticity sheet is a continuous 
function across 76. From (2.8.7) with gi, = $17. f = const., the equation of that 
surface 1s: 


FP, = 2, + J +o.) da -f = 0. 


The direction cosines of the inward normal to a surface n,; = Fiz Ki K, = [>,(P Ë 2] 12 
can be easily found by means of differentiation under the integral sign in (4.2.1). The 
condition of a continuous induced normal velocity component (4.1.1) gives 


Aus, = — a. Aan e, = bi guo, yə (4.2.2) 
s -1 —1 ~ Mei? 
with — 9. =—7,,nj5 = vv and b, = ninis = vivis, or 
As, = Ass, (vai — Hla Vas) T gus y Vaa — biva) ut. (4.2.8) 


The condition Ap = 0 leads to Ziu = Zi or Yiu = Y, on the surface, which gives 
QU AP mom, + Su rw m Pi, = 0, i = l to 8, (4.2.4) 


where it? wn = (Put Blees, Insert the value Agana, from (4.2.2) into (4.2.4), calculate 
from the resultant equation Agsy,y, and insert it into (4.2.8), thus obtaining: 


EE: Ais, — Agro, [Vor m KT’ Va; wi DU a bv,s)u 1] , (4.2.5) 
= [2U + off), a, Te, Pl] (4.2.50) 


The vorticity surface in the {z,}-space is given by the formula. 
P, = aşya" | (U +o.) dz, nf = 0. (4.2.6) 


t 
Dy means of partial differentiution one can find 


Pus = BP vel ay, = Ge = F s Üş = Po)” İb EM, = 24,945.24, 
= [1+ (a? — UP KI -t, 


which gives the formulas: | 
The, = — Va = 008 z; = Fam, K, = pMn,, = -eMvın. (4.2.7) 
Han = Va, = cos B; Fa, K, = Mn, = Mu (4.2.70) 


© š m Nas = Vas = COB y, = Fa, K, = Man, = Mvis. (4.2.7b) 
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For a flat vorticity sheet parallel -y,-axis, the above equations’ reduce to those derived 
in section 4.1. Inserting (4.2.7) into (4.2.5) gives : 


Ağı, == Dees H ` vlg? + DM, (4.2.8) 


which for a flat vorticity surface transforms into (4.1.6). In the expression 
Aga = Ags (2.9, 21), the variable z, must be calculated from (4.2.1). Next, 


(r7), = TuMİ(z— vi —(V— Yiv — (z — e, hun ə (4.2.9) 


and (a, — dı) is given by (4.1.8). Using the results of previous sections, by means of 
(4.2.1) and (4.2.6), one easily obtains D, = p, D, = — az), D, = —2,5,, and 


il i dë = | f Lu? nta siele, (4.2.10). 
T Ow T] 


where 2,4, =— Fims 2oy =— Fiy, (4.2.1). To verify this result, use (1.4.1), (4.2.6) 


: mene -1 mr — —1 — -1 ` ` 
with ET ww H S 215, SES TECH 2,y, ass H day, I H Bəy, and obtain D 


f f de = il J s [Le (u, + (p, ] dz dy, (4.2.11) 


(EI Tw - š 
which transforms into (4.2.10), 

The area of integration 7 wel? al 18 again the area of the projection on z,kor Jet 
plane of the area contained inside the curve of intersection of the voticity surface with 
forecone T belonging to P(z, y, z). In fæ,}-space the vorticity surface is given by (4.2.6) 
and the forecone by rz = (z — x,)? — (q — y,)? — (z — ə)? = 0. After performing the indicated 
integration in (4.2.6) one has to calculate the value of 2; and to insert 16 into r Saa to 
eliminate 2,. That equation gives the desired projection of the curve. In {z,}-space a 
similar procedure must be used. Calculation of ®, from (4.2.1) and mserting ıt into 
ro = (@ — 2) IN “yy? Lë si)” = 0, gives the projection of the desired curve on the 
(z,,¥y,)-plane. Of course, those curves are not necessarily hyperbolas of the second 
degree. Hence, one gets as the final result: 


279 = — Jf 1— Mui? [(e- 2) - (y - 9a 7 (8 7 #1) vis] Aw 
Tim 
+ Mp7 (uh vur” Agam, — (04 — a). la” + (ss? + (21m) tdr, dy, (4.2.12) 
with M given above (4.27), K, below (4.2.1), Fua, Foy, Fon by (4.2.1), v i= 
1 to 8, by means of n, below (4.2.1), a, v = w or l, by (4.1.8), pa, and Ziy, by 
(4.2.1). For a flat vorticity surface (4.2.12) transforms onto (4.1.9). By means of 
(2.8.5) one may obtain the magnitudes of the displacements of the axis of the core: 


e —1 
J, = 8,— (U t 9,2) Pryde, (4.2.18) 
TITR 
i 
Sue = SiTR— f (U c 9,2,) 7 9,4 d2,, (4.2,13a) 


TITR 
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which define in a unique way the shape of the space curve, the axis of the core. The 
quantity Ag,, on the vorticity sheet is equal to Dez TP, which quantity is given by the 
equation 1,„=I'P(2.3.1). and was widely discussed in section 2.8. 


4.3. Boundary Value Problem for the Lifting Surface. 


Ihe quantity p™ calculated by means of Picard’s iteration process, should satisfy 
the boundary condition, i.e., the solid surface is a stream-surface (consists of stream- 
lines). This condition will be satisfied in an approximate way by requiring that the 
quantity Ag,', to be used on the right-hand side of (4.1.0) to calculate g(t) should 
satisfy the boundary condition on the lifting snrface. To calculate Ay", n>1, the 
formula (2.2.2) will be used with ul) = Pl, vn zy 


} (n) — go us 
Wr we) z oi 2701 L and 


2, = —z, tan ö, De, m0). The values of gəl arc defined from the condition that 
the solid body be a stream-surface, 1.6., 

go = (reih, (4.3.1) 
with e, t) = g(x, 31 2) = — 2, tan Š,), 


and h, = h,(z,, y,, 2,) being the equation of the lifting surface. In the present case, 
hig, = tan 8,, and with dz, = — (tan 3,)dz,, equation (2. 2. 2) takes the form: 
zITR yiTR 
Aq) = (1— tan? i) f Ag") da, + f Ag dy, (4.8.2) 
ILE yITE 
Since 99 is identically equal to zero, one has to apply at the begining some forms of func- 
tions fü) fÜ, discussed in Section 1.8. One may use known resulta based upon the 
lineariszed theory and a few particular cases are cited below. 


(a) Unswept Wing of Infinite Span. 
For thıs so-called Ackeret's type of lifting surface 
iTR 
Ag, = (U2q)7! f Apdz,, with (Ap)q7! = du”!tand, = du 18, 
Ting 
so that when the leading edge coincides with the y, axis i.e., J eed dz,, the circula- 


0 
tion per unit span is equal to: 


Tz) z Ag,G! z 2Vu-!z, tan 8, ws 2Un- le 8), (4.8.8) 
and at the trailing edge or in the wake: 
T(z) = Ap?) = 2Upu-7!c, tan 0, == 2Un70,8,. (4.8.4) 


where c, denotes the uniform chord length. 
(b) Tringular Lifting Surface with Leading Edge Swept Back of the Mach Cone 
from the Vertex. 


Ag, was calculated in (Heaslet, Lomax and Jones, 1947), is cited in (Heaslet and Lomax, 
1948b) : 


De) = Ag) = 2U8,(E,u) (Gat ya” (4.8.5) 
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where E, is the complete elliptic integral of the second kind with modulus 

k, = (1—6,7)*, with 6, = p tan k, 
k being the semivertical angle of the triangle. At the trailing edge and in the wake the 
value of Ag is given by (4.3.5) with z, = c, (root chord). This represents the circu- 


lation per unit span at the coordinate y, 


(c) Tringular wing with leading edges swept ahead of the Mach cone from the 
vertex. ° 


Equations (4.8.8) and (4.3.4) may be used with c, denoting the variable length 
of the chord. 

(d) Rectangular, Trapeaotdal, Triangular and some other Plunform. 

These are discussed in (Heaslet, Lomax and Jones, 1947). 


In ease of an unswept wing, the total force perpendicular to the lifting surface is 
P = 4778 tan ö,, the lift is L = P cosö,, ep 4” sin ö,, ep= Cr tan 9,. 


AA First Approximation of Picard’s Process (lifting surface). 


e 


The following assumplion are made: Ap”) is given by (4.8.8, a, SÜ, v—uorl, 
the leading edge coincides with the y,-axis, the area of integration r, is the part ofa 
parabola contained between the following values of the coordinates: 

x-direction: x, 0, and with y, = ye: 

$14 = (1—p* tan? 8,) (z + a?ə tan 8,) - p[e*(1 — A? tan? 8,) +(z ez tan 60718, — (4.4.1) 

y,-direction: y, =y+e, and: 

Via = y tp İ(z — s)? — p?(s +z, (an 8,?]*. (4.4.2) 
In ease the (z,,2,)-plane passes through the point P(x, y, s), one has y = 0 in (4.4.2) 


Since there is a symmetry with respect to the vertical plane passing through P, with 
(4.1.9) and (4.8.3) one obtains: 





— a (KU)-1gU! = T, +I, (4.4.8) 
I 14 SI A 
‚= 2B tan 8,(z sin 8, +2 cos ô) lim f z (dz; = J. : r-8dy,, (4.4.80) 
€ 0 pre 
14 V14 
I, = 9B,u-1tan 8,lim | de, 11 r-idy,. (4.4.8b) 
ek 0 Ute 


The calculate the finite part of the integral in (4.4.84), introduce the symbol 
paz, tan’,)?—{x—2,)’= —F,'(z,, apply the procedure described in Section 1.5., and 
by means of elementery methods obtains the indefinite integral : 


[ext *4y. = 20e, + b)0X9)7 (4.4.4) 


where Zo = a+by,+cy,", a = Ri, b = 2g*y, c = —u*,q = 4ac—b*, or 
I, mag F(z) +c = (y- yi) - F,r) c. (4.4.5) 
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Since the function: A(y,) (1.5.8) is equal to unity, the Tayor s series reduces to this oné 


SS and (1.5.4) gives: 
C = lim [(-F’n)!(y-ya-y+y)]= 0, (4.4.6) 
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as one may verify by means L “Hospital”s rule. Thus, one gets from (1.5.5) with 
yı = yo as the lower limit: 


r * E 
ie il tr dy, = —eF,-?[(z—2,)* — p26? — p3(£ +z, tan 8,)?] 74. (4.4.7) 
Ute sc? 


14 | 
The integral f. A r"3dy, may be calculated by means of elementary methods sinoe it is 
St, - | 


finite at both limits. Hence 


^X-idy, = — (c) sin “"İ(2ey, +b) (-q)-3]*14 = spoz- .... > En 


te yt. : 
Inserting (4.4.7) and (4.4.8) into (4.4.8), one obtains the formulae: 
ef! = MO, + M,jIT,, (4.4.9) 
M, = 2Uu(l-+tan?3,)t[x sin?8, (cos 8,) 1 + z sin 8,] M, (4.4.9a) 
M,, = a[1+ (u”—1) sin”ö, |t, (4.4.9b) 
M, = 2U(p"+1)(1+tan?d,)* sın? 6, (cos 8,) 1 (à! M,,)7 1, (4.4.9c) 
eins are, ia, | (4.4.94) 
e>0 0 e = 
IL, = lim geb —sin- (ue, İda, ` (4.4.90) _ 
+0 J t 2 | 
where `. AE f , : : I 

Gy = (z—2,)* ps? — n (e +, tan ô)’. (4.4.9f) 


For 8, small, sin’, = tan"6, eat, sind, = $), and (4.4.02) and (4.4.96) transform 
into the expressions M, = 2Uud,e/r, M, — 0, which are known from the linearized 
theory (Heaslet and Lomax, 19488, p. 16, eq. (16)). To calculate II,, decompose F,’ into: 


Pä = (1-p tan *8,)[2,—(e©—pa)(L+p tan )-*] [w-(e+pe)(1—ptand,)-]. — (4.4.10) 


Applying the formula f e,(z, —a')-!(z, —57)- G-!dz,, one obtains II, = lim e - (L-1), 


where 1,,, La are solved by means of formulas f (vG)-tdz,, and sin-!z—sin-iy. With 
e—0, one gets II, = III, — IIL,, and a m 

III; = A, sin"! [C,(1— D3)& ~D (1— C3 t]. (4.4.11) 

A, = F (t pa) [2p (8 +z tan 5) +p tan 8,)] 1, (4.4.11a) 

with Ci: = Du =—1, i 1,2. The upper signs refer to i = 1, and the lower ones to 
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i 2. The value sin~'0= was chossn since the value O gives II, =0. . The final 
result is: : 
MI, = —N,[(z — uz)(1- p tan 6))71 + (£+ u2)(1 — a tan 8,)71], (4.4.12) 
F, = U(1+tan*8,)* sm ôu [1-- (p? — 1)sin?6, 1-5, (4.4.12a) 
For ó, small, sin?8, =s tan?8, =s tan 8, =s O, sin 6, ^&40,, (4.4.12) transform into M I= 
— Uà, n 1[(z + uz) + (z — nz) ] , which is the same as (18) in (Heaslet and Lə nox, 109495) 
If in (4.4.11) sin-!0 20, then M,II,— -N (z+ pz). (1— u tan 8,)^!, and for 8, small, 
M,II,— —U3,u7! (z + uz) (linearized theory). This gives w=) x — UĞ), wh.chis the 
same as (18) in (Heaslet and Lomax, 19485). The second integral (4.4 Ye) is II, 2I,, — I 
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TIA _ 
La = lim Í An d, = 1:(x— us) (L+ pu tan 8,) 1, (4.4.13a) 
a0 
0 l 
TIA 
I,, = lum sin-!(ueF,-!)dz, = 0. (4.4.13) 
«30 
To prove the last result, expand sinte = x-z’/6+-.. Since F, does not contain 


e and the integral itself is bounded for all 0 < € € co, for e0, TL, = 0. Thus M, II, = 
Nə Na. 


Fy, = Ufa? +1)(1+ tan?8,)*sin!8, (eos 6,)71(z — az), (4.4.14a) 
N,, = p?l1- (a? —1) gin?6, |*(1 * a tan 8). (4.4.148) 
Hence, 9 = M,II,-N4,N,,7!, with M,II, given by (4.4.12) and Na Nag by (4.4.14). 


4.5. First Approximation of Picard’s Process (vorticity sheet). 


In the initial approximation one cannot calculate the side-or-downwash since 
g%=0. The value of DÉI is given by I'OZI'P, P —1,i.e., by (4.4.4). Hence, in 
the first approximation the vorticity surface is a flat surface. The area of integration is 
a parabola consisting of two parts: oneon the lifting surface and another on the flat 
vorticity surface. Assuming that the region of iategration is not affected by the Mach 
cones from the tips of leading or trailing edge (where the situation is different), in the 
part on the lifting surface Ap is given by (4.8.8.) and ali by (4.4.3) with integrals 


f z dgr, in both (4.4.8a) and (4.4.8b), where 
0 
c, denotes the chord. Following (Heaslet and Lomax, 19488) write the integral 


| I EZ Apsudy, 


Il) 
Pi eee Piece 
But since the quantity Ag) r-9 < M is a bounded quantity for c, zzz, < To, where £, 
is a certam value however large, then changing the order of integration We may write 


in the form 


ite 
lim M dy, = lim Me = 0, 
«0 e 


V: 
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Therefore, for P situated behind tho trailing edge wave (envelope of trailing edge 
Mach cones) the contribution of the region lying on the lifting surface is given by the 


integrals " "Ne s. ! 2 
Fade. | | r”dın: To) r7dy,. (4.5.1) 
0 


y 

Yı 0 Yı : . 
The application of Hadamard's procedure shows that tho first integral (4.5.1) is equal to 
zero. The second integral is equal to 3u7"zc, and in the result 

° of) = Ulp? + 1)(1- tan? 8,)!? sin? 8, (cos 8,)"10,27", (4.5.2) 

Z = pg 114 (gu? — 1) sint 6,179. (4.5.2a) 

The expression (4.5.2) vanishes for 8, = 0 (linearized theory). To calculate the contribu- 
tion of the part lying on the (flat) vorticity sheet, equation (4.1.9) or (4.4.8) will be used 
with As) = 2U0x”1c, tan 8,, Ar). = 0, (hence, I, = 0). The rogult is: 


—z(KU)-9$) = 2B tan à [x sin 8, 4 2 cos 8,1. 


; TIB 
.]im o, f de, 


ch 


VIA 
J r-*dy,. (4.5.8) 
Yere 
The finite value of the integral was calculated in Section 4.4 Inserting (4.4.7) into 


(4.5.3) and applying elementary formulas (seo above eq. 4.4.10), after some trans- 
formations gives the result: 





d 


ell = — 2x sin 8, +# cos 8,)c,U tan 5,H 1, (4.5.4) 
H = n cos 8,[14 (x? —1) sin? 6, 11?(z + z tan 8). (4.5.4a) 


Higher approximations can be caleulated ın a similar way. 


The paper represents a resumé of investigations carried out for and sponsored by 
Boeing Aircraft Company. 


PANEL ON FLUID AND Borib MECHANICS, 
UNIVERSITY OF ILLINOIB, 
URBANA, Inrtrsoi8, U. B. A. 
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CALCUTTA MATHEMATICAL SOCIETY. 


Report of the Council for the year 1952 to ihe Arana 
General Meeting of the Society m 


The Council of the Calcutta Mathematical Society has the pleasure to submit the 
followıng report on the general concerns of the Society for the year 1952 as required by 


the provisions of Rule 25. 

The Council: The Council of the Society for the year 1952 consisting of the officers 
and other members elected at the last Annual General Meeting and co-opted thereafter, 
together with the Editorial Secretary, was constiluted as follows: 

President 


Professor H. D. Bagchi. 


Vioe-Presidents 


Prof. N. R. Sen Dr. T. Vijayaraghavan 
Dr. R. N. Sen Prof, A. C. Banerjee 
Dr. S. K. Chakraborty 


Treasurer 


Mr. 8. C. Ghosh 


Secretary 
Mr. U. R. Burman 


Editorial Secretary 
Mr. P. K. Ghosh 


Other Members of the Council 


Dr. G. Prasad Prof. 8. K. Basu 
Prof. B. R. Beth Dr. B. Ghosh 

Mr. B. C. Chatterjee Dr. B. S. Ray 
Prof. B. B. Sen Dr. N. L. Ghosh - 
Dr. N. G. Shabde Prof. N. M. Basu 
Dr. A. C. Choudhury Dr, Ram Behari 


Dr. S. 8 Bhatnagar (Co-opted, representative of the Government of India). 
General: The various activities of the Society have been carried’ on. throughout the 
year in much the same form as in the past few years. d 
Membership: The Council records its profound sense of sorrow und loss at ihe 
death of Prof, Elie Cartan of France, who was an honorary member of the Society and 
an outstanding malhematician of the present century. 
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The Council also reports that there have been some additions to the list of members 
of the Society, eight ordinary and one life-member being enrolled during the year. 


Meetings during 1952: The council held six meetings during the year and there were 
five ordinary general meetings devoted to the reading of original papers communicated to 
the Society for publication in its Bulletin. 


Publications: During the year under review the Society has published four 
numbers of the Bulletin namely, vol. 48, Nos. 2-4 and vol. 44,No. 1. The Council regrets 
very much that on account of circumstances beyond its control the publications of the 
Society are being very considerably delayed and desires to report that with a view to 
expedite the publications the Society has found 16 necessary to appoint a trained compositor 
with effect from July last on a salary of Rs. 75/- per month. This compositor is now 
engaged at the University Press exclusively in the work of printing the Society’s Bulletin 
and it is expected that if funds be available to enable the Society to maintain this arrange- 
ment, the arrears in publications will be wiped off before the current year is out. The 
Council takes the opportunity to record here the Society’s indebtedness to the authorities 
of the Calcutta University for printing the Bulletin free of charge and the Officers and 
members of the staff of the Calcutta University Press for their valued services. 


Exchange of Publications: The transmission of the Society’s publications to various 
countries in the world has been carried on regularly during the year and some new 
exchange relations have also been established. 


Library: It is gratifying to note that there has been no break during the year under 
review in the run of periodicals which the Society has been subscribing for the last few 
years, out of grants received from the Central and State Governments. These granis 
have also enabled the Society to add to its library a: collection of some recent standard 
books on various branches of Pure and Applied Mathemetics. 

Ihe Couneil desires to report here that a book case has been constructed for the 
library at a cost of Rs. 500/- to provide space for the current issues of the periodicals, 
prior to their being sent out to the bindery. 

Finance: The annual accounts of the Society for the year under review have been 
presented to the Council in the standardised form by the auditors Dr. H. M. Sen Gupta 
and Mr. M. C. Chaki. The Council gratefully acknowledges its indebtedness to them for 
their honorary services. | 


The Society received the following grants during the year under review :— 


(i). Government of India Se ön Rs. 1,000/- 
(ii). Government of W. Bengal .. s Bs, 1,000/- 
(ii). Nat. Inst. of Sciences of India ... Rs. 400/- 


and the Council takes this opportunity to offer its grateful thanks to the Central and State 
Governments and to the National Institute of Sciences for these grants, 
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VARIATION OF TEMPERATURE DUE TO SMALL STEADY 
DISTURBANCES IN A COMPRESSIBLE FLOW 


Br 
M. Ray, Agra 


(Recsıved— April 7, 1968) 
Summary 


The problem of temperature distribution due to small steady disturbances in a uniform stream in two- 
dimensions has been discussed. The fluid has heen taken slightly compressible, the effects of viscosity and 
conductivity have been taken into account and the condition has been taken as non-ediabatic. Usual method 
of neglecting squares and products of the components of disturbed velocity has been adopted. It has be&n 
fcund that variation of temperature is proportional to that of velocity and that the disturbance depends only 
on one parameter, that 18, Prandtl number. 


1. Introduction: An infinite mass of fluid has a uniform temperature and flows 
with uniform velocity U in a given direction. In such a fluid small steady disturbances 
are imagined to propagate, which give rise to small variations of velocity and temperature, 
Taking into account the effects of viscosity and conductivity, the disturbed velocity 
will be rotational in character but if this disturbed velocity is small, tha equations of 
motion can be simplified and under some plausible assumptions, can be solved. Temperu- 


ture distribution i8 given by the Energy Equation and the above assumptions arise 
when we go to solve the Energy Equation. ^ 


4. Disturbance Equations: The direction of the uniform undisturbed stream is 
taken as x-axis and the flow is supposed to be non-adiabatie and two-dimensional. If the 
disturbance is small and the fluid is very slightly compressible, then the variations of 
velocity, temperature and pressure will be all very small. We take pressure as constant. 
Let (U+t, v) be the components of disturbed velocity and if e, be the dansity in the 
undisturbed state, let (e+e) be the density when disturbed. The quantities u, v, o 
are all small so that thei squares and products can be neglected as well as the squares 
and products of their first partial derivatives. 

The equation of continuity is 


. fleo + 0)(U +u)} +s {lea rock = 0, 


which simplifies into, with the above assumptions, 
Gu w) 9e 
— +— )+ U— = 0. l 
i (= Oy/ ” “Ər ü) 
The equations of motion are 


(ea +0) (Uu) m = IRIS , 
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ü leo +e) (Uv) CMN È pm 
ED 
" | Pw = dE EM Si 
With the above approximations, these equations simplify into 
| ee -—- 2 
and 
et zr He) "as (sv) e 


8. ‘Energy Equation: IfI be the enthalpy or heat content per unit mass in the 
undisturbed stream and (I-i) ig the enthalpy in. the disturbed state, the equation of 
energy I8 

o Erd) = CH SERA 4 


where o = ucy„/«, is the Prandtl number, assumed constant; x, the thermal conductivity 
and c, the specific heat at constant pressure, assumed constant. 

It is to be noted that in the right hand side of equation (4) the dissipation function 
has been neglected in the present approximation as it consists of terms which are 
squares and products of partial derivatives, so that no heat is genora ed by dissipation 
for such a small disturbance. | 

‘Lhe equation (4) simplifies into 


iB - 2 (2) +2 (8) s 


Here again i and iis partial ned have been taken as small. 
4. Solution of the Equation: If we take v= 0, tbe equations (1), (2) and (8) 

become ; 
Du zz 7: ux = E r (6). 


t 


Qu 4 Of ðu 9, CS) 
ou m 2 


and TET n A " 
SM 8 dy V Gs, (8) 


{ 


VARIATION OF TEMPERATURE DUE TO SMALL STEADY DISTURBANCES ETC. 47 


The equation (8) is satisfied 1f we put ` 


Me et potaa OV (9. 
"5s Se ^ —— "Sy $ By’ 
then the equation (7) gives m 4 90) 2 56 (10) 


"0705 8 On Ge 8 öy Oy 
Now assuming ı to be a function of u alone, the energy equation (5) gives 
ren USL. du Su 8 („A BY) EE 
° "dn Ser Əz õu V du öz/ ` öy Su V du 5)” 


and therefore from (9) we get 


vğbdi . ə 9 (ək di), EYED m 


öz du ` Ge Ou Vöz du 9 8y du Ou du 


To solve the equation, «e furlber assume that öy /Əz and Oy /Əy are functions of u alone 
Then the equation (10) becomes 


y du oy 
—- eat = 484.4 (99), 4 99 d (ov) 
. "Ge 8 özdü Vöz/ 9 8y du \dy 
-2 x (Gb 2 d ST 
8 du NOx 9 du Vöy 
"$a +] 
9 neg * ay) VK 
Therefore : 
4. (EN - 9, ov a3) 
du Öz Oy 2 Ox 


Now to soive equations (11) and (12), we inust make some plausıblı sassumptiuns 
regarding Ov/Oz and Ov/Oy. Since these are functions of x alone, the simplest and 
most suitable assumptions are to put, 


I + = AU, > = bu, (13) 
where a, b are constants. 
Using (13) we get from (12) 
12a? 403 = 9g, Ua. (14) 
Again using (18) we get from (11) 
dt 2 ALS d'r 2 
or : 
d?i di _ 
Wa + ..” mul. i (15) 


48 M. RAY 


where 


% = 70, Ua] (a* + $b). (10) 
The solution of (15) is | 


1= Cüsla, — (17) 
taking i =Owhenu=0, so thata >0. ` | 
Now for a perfect gas we have, | 
ie = yply-D, | 
where y = Cy/C8, c, being the specific heat at con-tant volume. This law in the present 
case gives, since p is constant, 





b (I+ 1)(0, +e) = constant = Io, 
or neglecting ie, we get 
0i t Io = 0, 
therefore : | 
OL -Op 
—+l==0. 
Coa m | . 
Hence with the help of equatıon (6) we get | 
vi ou 
Öz Öz | 
Since i is a function of u only, this gives, | 
di e 
U— =I 
du 
Therefore from (17) we get | 
UCw- = I. 
This can be true only if « = 1, so that | 
C =1/U . 
and hence from (17), we get 
/I=u/U. | (19) 
Now since « = 1, we have from (16), | 
9a? +46? = 900,0a. ` (20) 
We are now in a position to find the values of the constants a and b. 
From (14) and (20) we get, | 
a = Be U (l-o), i 
and 
f | (21) ` 
b = 235 aU viele J. | 
Hence we must have o <1 but > à, which is true practically in all cases 
Turning back to the equations (9) and (18) we get 
Qu Qu 
ə MP LCS 
SS au an ech 5 u 
so that 
Ön ' 
t ADU. o (22) 
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where 
A = Ba[2b. 
Using the values of a and b from (21), we get 
À = S$B(L-o)/(do —8)]*. (28) 
Now the solution of equation (22) is B 
u = f(y-F Az). (24) 


This gives the disturbance. R 

Thus if a disturbance represented by equation (24) is propagated throughout a 
uniform steady stream, the variation of temperature will be given by equation (19). 

For instance, 1f we take the disturbance as 


u = u, Bin (y+Az), (25) 
the lines of no disturbance will be y + Az = nr where n is zero or integer and the lines of 
greatest disturbance will be y + Az = $nz, n being odd integer. On the first set of lines 


temperature wıll remain constant whereas on the second set, variations of temperature 
will be greatest. 


6. Result: When a small steady disturbance is propagated into a steady uniform 
flow of a slightly compressible fluid, tbere will be a small variation of velocity and 
temperature. As a first approximation pressure will remain constant and there will be no 
generation of heat by dissipation and variation of temperature will be proportional to 
that of velocity. In particular we can take these variations to be maximum or minimum 
_ along parallel lines cutting the direction of the uniform stream at an angle, the slope 
depending upon Prandtl number. 
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TWO DIMENSIONAL STEADY FLOWS SUPERPOSABLE ON A 
SOURCE, SINK OR DOUBLET 


Ry 


Prem Prakasa, Lucknow. 


(Communicated by Prof. N R. Sen. Recewed—Aprıl 7, 1959) 


Introduction .—The present paper aims at investigating the general type of steady 
flows that can be superimposed on a two dimensional source, sink or doublet. It is seen 
that only an irrotational motion can be superimposed on a flow due te either of these; 
superposability being defined as earlier (Ram Ballabh, 1940). 


1. First we take the case of a source or sink. Let a two dimensional sou'ce of 
strength m be situated at the origin. The velocity potential 9 is given by 


pg = — m ]og ? = =" log (z^ + y?) 


If (u,, v,) be the velocity, we have 


TMZ m 
v = y 


pcs 2 HE: 2 
cz Ty cry 





Let (u, v) and Ç be the velocity and vorticity for a two dimensional rotational flow 
superposable on (u,, pu. The condition of superposability [Ram Ballabh, 1 c , (1. 8)] 


gives 
Of mz pU ci my \=0 
alge ty! QyNz! + ui d 





which gives 


Gç oC - 
+ — Q 1.1 
of which the general solution is 
= f(y/@), (1.2) 


“TT” being an arbitrary function. 


Ihe condition of integrability of the equations of rotational motion, assumed steady, 
is (Lamb, 1945) 


(ue +o orl = 0, 


v being the kinematic viscosity. 
Substituting for £ from (1.2), 


f'(y[z)(z?v —a^yu) = vi(z” - y?) - f yje) + 2zy - f(y/2)], 
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Putting öy 
dün 


ÖV + ory? = 200 +9) +2ayl. 
Öz Oy f | 
The auxiliary equations to solve this are 


dz _ dy ay, 


ək” nı — 


° g? göy Ga (gà + ui D dary 


two solutions of which are 


| 
y lz = constant, and 


—v 1+ 2) +2yje log x = constant. 
f i 


H.ence the complete solution is given by | 





=vi(l+y? at) OZ i 3 [a lo z xU İz), 1.8 
pt 2ulz) log sexs) (1.8, 

| 

where X is another arbitrary function. | 
The value of d as obtained in (1.8) has to satisfy the requirement 

V^) = t = fla) 


| ay 1%, 1 Oy _ 
ENKER ise? 


transformed into polar coordinates. 


Substituting for y from (1.3) and writing F tor 


fa + tan? T Gah + 2 tan 0 İ 


this becomes 
[ 


v log 7 0) [ F" sect 6+ 2F! sec? 0 tan 0] e Do per 6+ 2y' sec? 0 tan 6 
T 


—vF sec? 0 — 2vF' gec? 0 tan 0] = f( tan 0). 

Since there is no r on the right hand side, the,Jeft hand side must be a function of 

0 alone which is possible only if it vanishes identically, t.e., if f( tan 6) = 0, or the motion 
(u, v, 6.) is irrotational. | 


9. We next take the case of a two dimensional doublet. 


Let a doublet of strength o be situated at the origin, its axis coinciding with the 
axis of x. The velocity potential ọ 18 given by 





TWO DIMENSIONAL STEADY FLOWS SUPERPOSABLE ON A SOURCE, ETC. 53 


The velocity at any poınt of the plane due to the doublet is given by 


u aby") y SEN 
(+g 7 (ex 
If u, v, be the velocity components and £ the vorticity of the flow which can be 
superimposed upon (u,, v,,), the condition of superposability gives 


p(z? -9” of 2 ; 
ə GD ETH ER 


es 


Lë + y") Oy Uz” + y? . 
which reduces to 
ten —,2)95 + 22y = 0 2.1) 
Ox Ou |. 
of which the complete solution is 
= Xiy/(z" +y), (2.2) 


X being an arbitrary function. l 
The condition of integrability of the equations of UTİ motion assumed steady 18 


ae elt 
( Ba Oy vV C E 0, 
or substituting for t from (2.2) i 
VX! fyi (æ rulle —y?) -uX iyi (z* rant day = vX fy (1 +y’), 
and substituting for w and v in terms of the stream function d 
ON. ogy Sh X Juli +97) 
00 ss yy Cm 
ÖV xöv/(s”ag") > 
To solve (his, the auxiliary equations are 
de du _Xy/(at-+y%)} dy 
sy? ey vy" dy[( + ya} 


of which two solutions are 





Y _—_ 
(a? +y") = constant, and 
xy i t y?)] 
xs İy /(z” +y? jut .. 2 = constant. 
Hence the complete solution is 
X” ! ae 3 x 
POS EE = Giy/ (+y), (2.3) 


Xeiya mü müay” 
where @ is another arbitrary function. 
As before this value of y has to fulfil the requirement 





Vy = 0 = XİN /(z yt. (2.4) 
Transforming to polar coordinates the value of d as obtained in (2.8) we get 
sin 57 cos 0 (ee " 
Tv = š 
d Fİ T T x r / ES? 


2--1840P—2. 
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where "(S A 
F SN 


Writing (2.4) as 








a nd substituting y from (2.5), we get 
G" voos6 Hi = x (en 2) 





pt 
Since the only term containing cos 0 is 


_ v cos 0 Fi 


= , we have 


aka = 0, and thereafter G'' = 0, 


so that x = 0, or the vorticity of the flow superposable on the doublet is zero, that is, 
the flow is irrotational. 

My thanks are due to Dr. Ram Ballabh for his guidance and the Scientific Research 
Committee, U. P. for the monetary help given to me. 
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ON THE STRESSES DUE TO A SMALL SPHERICAL INCLUSION 
INA UNIFORM BEAM UNDER CONSTANT 
BENDING MOMENT 


By 
Sısır CHANDRA Das, Chandernagore, West Bengal 


(Recewed—March 23, 1952) 


The concentration of stresses ın the neighbourhood of a spherical inelusion under 
uniform tensıon has been obtained by Goodier (1988). The distribution of stresses round 
spherical and spheroidal inclusions in the case of an uniform state of stress infinitely far 
from the inclusion has been obtained by Edwards (1951). Stresses have alsc been calcu- 
lated by Das (1952) in the case of similar inclusions on the axis of a large shaft which has 
been twisted by torsional stresses. 


The object of this paper is to find the concentration of stresses due to a small 
spherical inclusion situated on the neutral axis of uniform beam under constant bending 
moment. The particular cases of a perfectly rigid inclusion and when the Young's 
modulus of the inclusion is one-quarter that of the surrounding material have been 
considered in detail. 


4. Method of Solution. 


The line of controids of the normal sections of the beam in the unstressed state 
has been taken as the z-axis. The z-axis is taken at right angles to this line in the plane 
in which the line of centroids is bent and the y-axis is drawn perpendicular to the axes 
of z and <. 


When the beam is subjected to a constant bending moment only we can assume 
all the components of stresses to be zero except Z, which is given by 


Zi = Az (1.1) 


where A is a constant depending on the dimensions of the section and the curvature 
of the bent neutral axis. Taking the centre of the inclusion which is situated on the 
neutral axis as the origin and using polar co-ordinates such that 


z = T gin $£ cos p 


y == + gin 0 sin 9 
and 
z = r cosb 
we find that 
Z, = Ár sin 6 cos ç. (1.2) 


I 
l 
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The stress components in polar co-ordinates correspondmg to the above state can be 


put as 
... a " 
üə [8Pi(cos 6) + 2Pl(cos os 6)], 
> zer 157 
00 — = [6P:( (cos 0) — Pa(cos 6)] cos 9, (1.8) 
= BAT xm (cos 6) _ 6dPı (cos Josi 
= din (COS 0) 

Il də K. en 

and x sön. 05 


where Ph (cos 6) is a Legendre’s associated function of the first kind of the n-th order 


^ 


and the m-th degree. 


The equations of equilibrium in polar co-ordinates are 
l 


9A o, Te, 9. : | = 0 
(A+2u) sind — Zu E “ə sin 6) 


| | 
Dad 1:95 94 3 ta- ==] = 0 





sın Ó Oe Qr 
and | 
I Ë 9-99 | = o 14 
(N+ 2u)r sın n. 2u T Les, sin 6) T (1 4) 
where | 
A= — [2 (r*u sin 6) E " sin += (n |, 
1 [2 AES | 
ə. et 
2; "urs lae (rw sin , Bg (rv) | 
1 au o | 
Zug = —— 6 1:5 
Ge r sin 0 Pa ar d ne) —- 
and | 


SE "lə: )7 P 


u,v, being the components of displacement in the directions m which r, 6 and 9 
increase respectively. 


As particular solubions of the equations (1.4) we take 


A =[A,rPi(cos 6) + ^: Pi (cos 6)] cos 9, 


ges [5.5 | Pl(cos 6) sin 9, 
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ous; LA, An Ps (cos 6) +( Be - 5 ye (cos 4 ə 


2 art sin 6 2 Br“ dé 
and | 
"T 1-4: dPi (cos 6) , A, dPş (cos 6) 
I 2 do ort do 
By € ) PE 008 e .. 
MUSİ İla 9 un cc owe n 1.6 
( 2 Ar? sin 6 SSC x 
where a='A+2u/u) 


Assuming the inclusion to be very small compared with the dimensions of the beam, 
our problem is to find the components of displacement %,, v,, w, that will satisfy the 
equations of equilibrium given by (1:4) end will yield stresses given by (1,8) 


when r>o., (1.7) 
Also for the inclusion the components of displacement us, v,, w, must be such that 
AL = ty, (1.8) 

0, =, 

and w, = WwW, 


when r =a, a bemg the radius of the spherical inclusion. 


Further when r = a. we must have 
[= [n]; 
[6], = [76], (1.9) 


[re]; = [79]; 
where suffixes 1 and 2 denote stresses calculated from u,, v, w, and tiz) v,, W, respectively. 


2. Solution of the Problem. 


The components of displacement satisfying (1.0) can be taken as when r œa 


| A,(8— 2E, 
u = [18-9 pa zs “51 [Pi (cos 6) cose 


+|% 3(1+ 2a) EC 
51? 


T 0D,r?-— 





| P (eos 6) cos 9 


v, [4799 2 Bari, Gy Lier? Lus, 


VE d6 cos +| 6 sın 6 


10 r 
E? 24,(8-—a) 14Pİ(cos6 
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c 1 
~ P SÉ Pi (cos) „, 


` In ọ — 
r. sin 6 


Bir, 0, 122 (cos 6) 
[ e ch — Q sın 9 (2.1) 





and when r < a 1.6. within the inclusion, 


u= [4:895] Pi (cos 0) cos o + BD? P (cos 6) COS 9 
| E 
A,(8 — a) ts 6) Bar? P3 (008 6) ' a APs (cos 6) 
=| 209 “7 e) Ze £3 008 6) Dar? š 8 \cos € 
73 [ 10 E ep nn er 


; Pl Ps (cos 6) 


, 1 
w,= [- 9 Je 9) uin ç — gin 9— Dyr Bin p (2.2) 
8 sm ` 


Bar? dP loos 0),. 
10 in 6 6 dé 


where x” = (A’+2n’)/p’; A’, v being elastic constants for the inclusion. 


The stresses calculated from the displacements (2.1) are given below: 





IF =[4 Br ran, pu İF (cos 6) [ee HM. 19, Dyr+ on İF (cos 8) 
89 r 


T -[-4 (8A + 24), Ob, SE MİPifsos 0) | 
COS 9 10 r* dé | 





+ [2 SC Paloos6) , Jh pn, 10E, 10E, A, o (81-ə) PS eos 0) 
8r* sin 0 15r* dé 


EE 6p E, OE] Pi (eos 8 P (cus 6) 


sin ç 10 rt J em 6 | 


_ „Bar 390, as (cos 6) «[- m 4 As(8A+ 7) ] Ps (cos 0) 
| 6 874 dd “m” ən 1575 sin 6 


| 
SEN 2A, 2 ) “HB a+ Ou Dap | 


T A=) euE, -p0, M 1 "CN P} (cos 6) — - [140,7 ZC ET 
SIE EEN I 


A,(A- 1 
-SLAA + Gu, po, V - 9 


COB 9 





CoB 9 .. 





Ey İpi s (eos [Sr , Mc D [P$ (eos 6) 


nl, +80D,)r "10, - A,(8— dn 90s [sin ging (2:8) 
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The stresses calculated from (2.2) are as follows: 


rr ABA änt pi -—" 
cosp 5 TPı (cos 6) + 12u'Dsr Ps (cos 6) 


ED: 2.2 A (BA + Qu’) dPi (cos 9) Lb Bar P3 (cos 9) d Pa (cos 9) 


A Au Dar 
COS 9 10 dë 6 sin6 ann” 
7 : , , 1 , Dn” 1 y l 
re — A (SX + 29) , Pi (cos 0) u’ Be , d Ps (cos 9) _ au’ Dir Pa (cos 8) 





sin 9 10 sın 6 6 do n 0 





6 Ez (BA’ + Qn’) 


: m Dil”: (eos 9) - 122” Ds rP} (cos 6) 


5 kk... ub Gu! Di ePi (cos 6) 


cos o 5 
Ang [Bs--80Ds]i sın Ó cos 9 
sin o 


The conditions (1.7), (1.1) and (1.9) will be satisfied if 
TE 
` BA+ 2m’ 


A 
B,=- 2, 
3 Zu 


14A (u! — u)a* 
8[n(19A - 201) +4u (6A - llg)] ' 


- —A(A + u)(u! — u)a* 
E.— — 5505505... 
` Ba[n(19x +26) + 4! (A + 11n)]' 


7A(A+'20) 


D;= , 
° 18[p(19A + 261) T 4 J4X + 11p)] 


E ce [EE UO 290 ER] 
` 20u Din — A) + a’ (BA + 2u!) 


, 5 / 
AQ 4, -—“ 99MM _ 
27771 Butt — M) + BA” + 2p) 


24. | p'(p — A)(8A" + 20") — ala — A')(8A + 95) + LOu*(2A' +80) _ 90A + 96] 
15 


Bo= 
“(da n) Sala’ — A) + "BAM + 2p’) 
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- Alp EM Hp] — pt MBA U 
l0u(4u +u) Safa” — A^) + u’(BA’ + 2p’) 
_ Ap (21A + 194) — 5u(8A + 2u)} 
| € 


a 


| (2.5) 
Hence on the surface the stresses become 


( 5 ) “ə TAKLE ŞA məalı Ap’ asin g- | 999 399 | Ana gin? 0 





cO89/r.«a L (8A+2u)Q P 
— | 
( r6 ) = İ Duidu(u'—A')—u'(8A'—-29')--2R (992 --86a) 
COS O/raa 10(4u + u')(31 +2u)Q (Dän + u/)(BA + 2a) 
— 170. 29) ez pu! Aa cos 6+ IOSA + 2m) ul 2u)u Aa cos? 6 
P P 
pm | 
( To ) = | bulu”(9A” + Qu’ CAu(u/ —7))—2hR (99A + 86a) 2A +24) ¿Aq 
sin 9^r«a L0(4p + &)(8A -- 21) Q 75(4j. +o) (BA + 2p) p ^ 
2R +20 8p) 850.2») Dë 862) iM Ka 
b(A4p. + u )(3A + 2u)Q P '75ldu + p')(BA + Bu) Q 


RI _ (p'— p)(596A + 4849) _ (n! + Op) RBC Bp) 
G08 9/ reca 5P 6Q(p! + 4u)(BA + 21) 


| ; 
„ 5p(BA + 2p) — U (21À + I9 4a dnde yı “-44u)(u =u) ` 1] Aa sin? 6 
15(u! + 4) (BA +2) | P 


(=?) el KE E ‘Ou’ + By) R + 5u (2a + Ən) 
COS Q/rea 5P 5Q (nl + 4u)(BA + 2p) 


ur 179) RİMA OR | a sin 0— pil Aa sin? 0 
i P 
| 





15(u’ + 4a) (BA + 2p) 
( Z ) -[- (u/— 00A + 44u) _ v | fo (2A/ + 8p") 
sin 9/r.a P 10(3A +2u)(u’+4u) Q 
_2(21A +194) ; 10u(8A 4-24) i o don 
- 9 ! 3 -— 
where P = 3[u(19X 201) + 4p (AA + 115] 
Q = Baler - V) eu 20) 
and R = y!'(y — A(8A* + 247) — p(w! — A/)(BX 2p) 


To obtain the stresses when there is a flaw instead af an inclusion we put x = 0. in 
the above stresses and get 


T EA HI) ging 700.--p) 
a DAUNE AEN ga NST H 
E plea 4(19A- 205). TDA + 267 





Áa gin? 6 


H 
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r 





ues et en 
69 — SEA +IOp) A sin 6 sos 9, 
(ss 0“ Ae) E ES 


(rr) sa = Gel zu GETT = 0 
The above results were obtained by Sen (1988). 


We now consider the case when the inclusion is perfectly rigid. Here we make 


H '—oo and get 


5 sə) [u E + WA 2A) Aa sin 9— 2 Ara 4 Aa sin? 6 
cos 9/..s L(8A- 2g) (4A lly) 4 (4A 1l) 


— Ate y o 2 





_ (228A +2174) ` VOD + 22) ]4 94.85 At24) Aq cos 0 
Ed EC 152 T 11) o P 4 qarına)”. 





(228A +2174) _ "(A Qu) | 7 55 (A+2u) _ 2(1142. 4-71) 4a sins 0 
1508 +2u) B(4A+ 11a) 12 (4X + (4X+11p) 75(8A + 2u) 


S: 


5(4A+11u)  15(8A-- 2p) 4(4A + 11g) 





pee 149A--191u 189A + 224 ] 4a sin 0-[ 51A 442 -1] Aa sin? 9 


109% + Un , 27A +18p ] DC m 


“İl5GA311a) 16(8A- 2g) 44A rl) ` 


t ~ 


If we take A = u or Poisson's ratio e of the material = 1, the above sresses become 


(=) = 84a sin 0 —14a sin? 0 
COS P/r=a 





ei 3 SS $884a COS 0c 14a cos? 8 





DE — M Aa — 18140 sin*ó 
os m 155042 — 18442 sin 





( 06 ) = 34a sin gedo sin? 0 


COS o / rea 


GH = 194asın dän sin’? ` (2.8) 
COB 9/ rea 


oB 


BS 69 ) = 14a sin 9 cos 6 : 
Bin O/rea d 


8—1840P —2. 
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Now putting A = — M/I where M is the bending moment about y-axis and I the moment 
of inertia of a section about an axis through its centre of gravity parallel to the y-axis 
we get the numerically greatest value of the above stresses given by 


(=) = 0'59 E? when 0 = —35'5° 
COB p/r= ef 


(- 5 ) = 0'75 Ln when 0 = 58° 
COS 9/r-a I 


ta zd 2. ende 
81D 0/rea I 





( 22 ) = 0'27 Za when 9 = — 88° 
yed 


COB 9 I 
(2) = 0'55 Mi when 6 = —51'8° (2.9) 
COB Ọ/raa I | 


(2) = 01 Ma, when 0 = — 45° 
SIN Qrad I. 


The ubove stresses when there is no inclusion are given by 





( d ) = 0:46 “İz, when 0 = —88'5° b 
COB D/r=G I 





( re ) = 0:25 Ma, when 0 = 58° 
yad 





COB 9 İ 
( 06 ) — 028 M, when 9 — — 38" (3.10) 
COS 9/ res I 


and = - - 
(r9), -a = (99 eng = (Op)roa = 0 
If the inclusion consists of a material having the same Poisson’s ratio + butan 
Young's Modulus one-quarter of that of the surrounding metal taken as a plausible value 
for slag in steel the stresses are given by 





( 7 ) => Aa sin eset. Aa sin? 0 











COB 9 16 
7 

a x = — Aa cos 0+ e Aa cos? 0 

CO8 Q/ roa 16 
A re ) = Aa + ida sin” 6 (2°11) 
BIN 9/rza 
( 2) = — 75 4a sin 9+ 22 Aa sin" 6 

COS Q rea 16 
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Lie IT gin 6 4- Aa gin? 0 
cos 0 7 ed 2 16 





zə TIME 
-——— O 0 
(= Q/rea 840 Ge 


Proceeding as in the previous case the numerically greatest values for the above 
stresses can be calculated to be 

















( ra ) — 0:24” a, when 6 = 40:89 (2.12) 
COS 9/ raa I 
( 2 ) 0:26 a, when = 180° 
CO8 o fot I 
( 77 ) — -020% a, when 6 = 90° 
B10 Q/ rsa I 
( ó ) = 18644, when 9 — --90" 
COS 9/ rea I 
( TT ) = 0:206. when g — 38° 2.12 
S 9/rca I 
(2) — 0.0784, when 9 =—45° 
81D 9/ rea I 
Also here 
(7535) = = 1,46 sin 6+ ! Aa sin? 6 (2°13) 
COS 9 84 4 
This will have the numerically greatest value when 0 = —z/2 and the value is 
given by 97 M 
75 2.14) 
687 an 


The numerically greatest value of the difference of the tensions are obtained at the 
point 
p = 0, 6=-7/2 Or Q = T, 6= 7/2 
and if fracture occurs, 1t commences at that point. 
In conclusion I thank Dr. B. Sen for his kind help in the preparation of this paper. 
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THE EFFECT OF TWO EQUAL CIRCULAR HOLES ON THE STRESS 
DISTRIBUTION IN A BEAM UNDER UNIFORM BENDING MOMENT 


By 
A. M. SENGUPTA, Sibpur. 
(Received —F'ebruargy 18, 1962) 
1. Introduction 
The problem of a circular hole with its centre at the neutral axis of a beam under 
uniform bending moment has been discussed and the result obtained have been verified 
approximately by photo-elastic measurement by Z. Tuzı (1928). The object of this paper 
is to study the influence on the stress distribution of two equal circular holes with their 
centres at the neutral axis of a deep beam of rectangular cross section under uniform 
bending moment M. For ths purpose bipolar coordinates are introduced. 


2. Bipolar coordinates and the expression for the stresses. 


Let us introduce the curvilinear coordinates defined by 


z-Fi(y +a) | (1) 


where z, y are cartesian coordinates and a is a positive length. 
Solving for z and y we have, 


asin D a sinh a 
T = „en = — O VP.s. RL. 2 
Y cosh x — cos 8 (2) 


cosh o — cos 6” 
i -(&) 4 E 
Rh? \da Qa 
we get, h = (cosh « —cosh 8)/a (8) 


The curves € = constant are a set of coaxial circles having the two poles (0, +a) for 
the limiting points, and 8 — constant another set of circles passing through the limiting 
points, (0, +a) and intersecting the first set orthogonally. The axis of z is given by «20. 
The radius of a circle z is g cosoch z, while the coordinates of its centre are (0, a coth a) 


From the relation 


We have 6 = 0 on the y axis, except on the segment joining the limiting points (0, + a) 
where 8 = +7. 

At infinity « = 0, ß = 0 and at the limiting points z = + oo. 

The stresses in terms of the stress function x ara given by (Jeffery 1921, p. 269) 


— ı 
aza = {(cosh «— cos Bas —sinh véi Sin e5 + cosh a) (hx), 


aßß = {(cosh a — cos £) ə — sinh en —sgin Bx + cos Bl(hx), 


axB = — (cosh « — cos B) 5 (hy). (4) 
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So that, we have 
qəs — Be) E pio - S +14. | (5) 
38? Ja? 


The necessary and sufficient conditions for a boundary a = constant, to be free from 
stress are (Jeffery 1921, p. 280) that on the boundary 


S) = constant = o (say) (6) 
and hx = o tanh «a +o (cosh a cos B —1) +7 sin B, (7). 


where o, or, 7 are Michell's three constants of the boundary. 


3. The effect of two equal circular holes on the stress distribution in a beam 
under uniform bending moment M. 


Let the holes be defined by a = c (750) and x =-c. The origin of coordinates 
and the centres of the holes being at the neutral axis. 


For a beam ot rectangular cross-section of thickness 2b and depth 2d subjected to 
a constant bending moment M, we have, the stress function 


. ud | 

X,——z 8) 
E 

giving te = xy = 0, yy = Az, where 


AB : 


a? sın? 8. (9) 


So that, 7 A 
6 (cosha-cos'ß)? 


and this has different expansions on different sides of (helinea=0. Thus 
for o > 0 


| hX, = $A a? 677“ gin B --$A a? > (2 cosh« — n sinha) e "e gin n f, (10) 
nc ' ` 
and for o < 0 I 
hX, = $A a*e?- sin B+4A a? ` (2 cosha+nsinha) es gin nf. (11) 
nm? 


For the complete stress function we have to add to X, a stress function x, which 
gives no stress at infinity and is such that the complete stress function (x, +X,) will give 
no stress over a = +c. 

We assume 


hx, = C, cosh 22 sin B+ > [Cn cosh (n + 1)« + D, cosh 1n —1)«] ein nf. (12) 
n=2 


so that the stresses due to X, at infinity vanish. 
The complete stress function X = x, +x, given by (10) and (12) satisfies the boundary 


conditions (6) and (7) for the boundary o = c. So that if o, o, 7 be the Michell’s three 
constants of the boundary o = c, we find without difficulty that o = o = 0, 
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C, cosh 2c + 44 oeh = r, 
“2C, sinh 2c — A a? 67% = 0. (13) 
For » > 2 . l 
On cosh (n +1) c + D, cosh (n — 1) c = —#Á a?(2 cosh c —n sinh c) 67%, (14) 
(n+1) C, sinh (n+1)c+(n—1)D, sinh (n — 1)e = —4A a? | (n? + 2) sinh c —8n cosh c]e-**. 
For «— —o, the stress function calculated from (11) and (12) satisfies the boundary condi- 


tions (6) and (7), so that we get the same system of equations (18) and (14) where p, c, 
7 are now replaced by three other constants plo. 


Solving equations (18) and (14), we get 
: C; = $A a? e—* cosech 2c (15) 
Forn>2 _ 

C, [sinh 2nc +n sinh 20] = 14 a®[n{n + 2) 6-26 — (n — 1)(n —2) +(n+2)e-?re], (16) 
and — D,[sinh 2nc +n sinh 2c] = 44 a? [n(n — 9) er — (n * 1)(n 4-2) — (n - 2)e7?]. 


On the hole « — c, an = O so that BB can be very easily calculated from (b). 
We have, D 
BB, = 24 a(cosh c — cos f) [cosech 2c sin 8 —3 2, B, sinnß], (17) 
where 
İRA +n(n—2)e-@-N0) (sinh Inc +n sinh 2c) 
—n[n(n —2)677 — (n + L)(n +2) — (n — 2)e “79 İcosh (n — 1)c 
The series in (17) is slowly convergent for small values ofc. To make it rapidly 
convergent, we put 
fy, = 2(n? sinh c —8n? cosh c +n sinh o)e-" + M, (18) 
where, 
M, [sinh 2nc +n sinh 2c] 
= ne "ler + a 737^ —9n sinh 2c] [(n? +2) sinh c — 8n cosh c] — (n — Din —2)e^ sinh 2c], 
Remembering that 


4( cosh c — cos ay 2n(n* sinh c —8n cosh c +2 sinh c)e “10 sin ng 


ne] 


SÉ sin D 


so co [cosh 2c cos*8 —2 cosh c cos 8 +1— sinh*c |, 


we have, 


25— sin 8 (cosh c — cos 8) 
Rose Saa? esinhQc 
sin 8 


PRU: [sinh?c +2 cosh o cos 8 — cosh Oo cos? 8-11 
cosh c — cos 


— $ (cosh c — cos 8) > M, sin np (19) 
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The table below gives the values of the stress BB on the boundary of the hole a = € 
for positive values of 8 between O and s when c = H From the same table the values 


of Bß for negative values of B are obtained by changing the sign of BB. 
Table. 


g = E 


—0' 1636 


—0'0886 


—0'0298 





For the purpose of comparison with the corresponding case of a single hole we 
introduce the angle 6 between the radii of the points (c, 0) and (o, 8,) on the boundary 
of the hole « = o. So that 
cosh c cos 8 —1 


d es i 
H cosh c—cos 8 


It is found after calculation that the maximum stress + 1'901 RA (R being the 
radius of the hole) takes place at 8 = +42° approximately i.e., at 0 = 90°36" approx. 
as compared with +2 RA at 7/2 in the case with a single hole "There are six zero values 
of this stress at the points 8 = 0, m, +18°10’ and +86°28° approx. ie, at 0 = 0, 
m, +45988' and +186° approx. respectively, instead of at 0 = 0, m, +4r and Län in 
the case with a single hole. 


In conclusion I thank Dr. B. Sen, for his kind help 1n the preparation of this paper 
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SOME INVERSION FORMULAE FOR THE GENERALISED 
HANKEL TRANSFORM 


By 
R. P. Acanwar, Lucknow, India. 


(Communicated by the Seoretary— Received February 24, 1952) 


1. Introduction: In a recent paper I gave a generalisation of the wellknown 
Hankel transform by solving the integral equation 


ite) = à» Í (ey sif" atıldı 0) 
where à 
H - (—a)* 
(x) SS > riL A+ pr)’ H > 0. 


I showed that the reciprocal kernel was expressible in terms of a function of the 
type J#(z). In fact ib was proved that 


g(y) = a 9220-3143) fa f (ay BAA TR [al (Qty?) tr] fd. 
0 


In particular, when x = 1, the above integral equation reduces to the Hankel 
integral equation. We shall call f(x) the generalised Hankel transform of g(y) in (1). 

The object of the present note is to find two more inversion formulae for (1) in 
which the reciprocal kernels are not functions of the type Ji(z) but are different. 
In proving these inversion formulae we shall use some known results which, for the sake 
of elegance, I cite in $2. 

2. Gupta (1948) has proved that 
f 2 IP etes (—)'a 104189) fada 

= wc ba uw) (2.1) 
valid Tor Ausl, Rea>-1 and in the case ps 1 an additional condition 
Re(4k+4l-2A+4m+1)>0. Also 


rw) 


J ten) Se CO ane -8-r-e-1dg = QwieK (2008) (2.2) 
Te 
provided 0 < u < 1, ReA>-—1and Re(4p+2A +8) > 0 in the case when u = 1. 
Meijer (1933) has proved the following inversion formula for his generalisation of the 
Laplace transform : 


If 
ug = [etsi əllə 
0 


4—1840P—2 
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then T(1-k+m) 1 


Bi 
g(s ) = lim ve T(1-2m) As Irid et My Ls (t) (st) -*- f (tdt 
A 


| 


provided tbat — 


(i) J e-&* | g(a) | ds is convergent for a > 0, 8,> a, 
0 | | 
Gi) g(s) is defined for positive 8 and is integrable in the Riemann sense in every 


interval 0 < T, < t < T,, | 

(iii) g(s) is a continuous function of 8, | 

(ivy kama-k, B> a. 

The last result that we shall need is ene to Boas (1942), who has proved the 
following result : 

Define the operators W and Qk as 


W[g(s)] = ei" [21-7 "p (ƏV 


and 1 
Qıl/(u)1 = (2k) | (2k | 9*3 Wb [2"-4f(«)] t=?kju 


where thë dashes denote differentiation with respect io g. 


If 
fis) = aries f HE, ay E 


is convergent and if ¢(t) is integrabie in every finite interval of t, then for almost all t 
= Well 


when —4 € Rey < 4. 
8. The first inversion formula: | 
Let us consider the integral equation (1). Multiplying both sides of this by 
iz, i= 9X1 / part) [ n-X—8 [218g [4+4pA—A—2 X | 
| 
Ta-k+m- än ah +247) —k— m= $n 4p 2*7) yep tote] le, 
TrID(A—2k—3u — uA 2 cr A x 7)/ id 


and integrating from 0 to oo with respect to x, we have ` 


f Hehite, t)de = (1/2) f (a, t)da il (eyy ie y*)g(y)dy. (8.1) 
0 0 0 


If the change in the order of integration is permissible! we get 


x(t) = f fite, Ge = (4) fy Mtg (y)dy f HRG) Ye, Ode, 
0 


Integrating the inner - integral- with help of (a. 1) we get, after a slight change 
in the variable, 


x(t) = Í (6077 5e 9 Lua mot )g lon) ds (3.2) 
0 
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valid for +Sus1, ReA>-1 and in the case a=1 with the additional condition 
Re(da + 2uX — 2A + 4k +4m — B) < 0. 
Now applying the result (2.8) to (8.2) we have got the required inversion formula 


sos 55” Xə 
Lay ha əə ce 48t(9f)-k-4M,_4 n (st)ytü)dt 
g(sta) 752 Got Gen? a ett(st) 74m (8) x(E)¢ 


where x(t) is defined above. This gives 


ig 2 Tü -k+m) f" yoqa ts?  ul(ta? fa)i-ky()d 
— —— MRCUUEM. 3 tn s 4 i 3.8 
24 M M 
zə 

provided that g(s*#) satisfies the conditions of Meijer’s inversion formula mentioned in 
$2 and x(t) exista. 

It now only remains to justify the change in the order of integration ın (3.1). To do 
so we see that (Wright, 1985) 

die, t) ~ gi (65-2484 - 283) [1g Bet nh ah Al exp ( — karte) 
for k > O and |z | large, and 
di, t) = Ole / 28) QA +42 — Ba) {4 (Bat 2nd — 42-8) 

for | a | small. 

The change in the order of integration is now valid provided 

3 g(y) = O(y*), Ref <-1 for | y | large 
an 
I = O(y*), Rea, > —1 for | y | small, 

with3<us1l, Re(z +À+8)>2> 0 and Bei > —1. 

Thus our inversion formula (8.8) will be valid under the conditions 

tü) Ra>-1; k —<u = 1 and in the case p =1 with an additional condition 
Re(8p +2 — 2A +4k 44m — 5) <0; 

(ii) gly) = O(y£), Reg <—1 for İyl — co 


= O(y*), Rea >—1 for İyl—0 
where Beie + À + 8) > 0, 
Gi) g(yt*) is a continuous function of yiz and is defined for y > 0 and is mtegrable 
in the Riemann sense in every-intervul 0 < T, < t < fa 
(v) x(t) exists. 
4. The second inversion formula: 


As in §8, let us next multiply both sides of (1) by [Gupta (1948) 1. 


p^ 1p? 79779 CEET NEE 


eo EY, 
+ where Gic m rectc, a € 1 is the generalised parabolic cylinder function given by 


Gupta In particular 
oe 97 T (0X)D Als) /9)ef* . 
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Integrating both sides with respect to z from zero to infinity we get 


p [ftat aba lp? at) 
I | 


= ($) pmp?” J Olli 79 -1da J yt*ia*yt)gw)dy — (4.1) 
0 | 0 
Changing the order of integration, if permissible, and integrating the z-ıntegral on 
the right with the help of (2.2) we get for O «Zu < 1, Bei >-—1 and Re(4v+2A+8) > 0 
that ! 


9X vk 


| 
p pri J eecht a ip" (x?) da Si y il (yp)tK,(yp).yAtrg(y)dy. (4.2) 
0 0 





Now, using the inversion formula (2.4) we get, for almost all y 
nell = lim Q,[y(p)] (4.8) 


where 
Qr Y 


(p) = 
and Qk is the operator defined in (2.4). The above result holds provided that 
(i) the integral on the right of (4.2) is absolute] y convergent, 
(ii) dun) exists and y"*^g(y) is integrable in any finite interval of t, 
Gü) —4<Bev<t. | 


In order to justify the change in the order of: integration we see on studying the 





ıl +a Tee 
: | 


behaviour of Giel for large and small values of [z |, that for 0 «uz 1 
Gf(z) = O(z-^I») for |z]— co it A/p sb an integer 
= O(1) for |x| 0. . 
Hence, the change in the order of integration wiil be permissible by absolute 
convergence if | 


(i) ReA2-1, 
with | 
(11) g(x) 


O(z8), for | 2 — co 
O(z*), for |z |-»0 


ll 


where Re(2v+ z+ À +8) > 0. | 
Thus our second inversion formula (4.8) is valid under the conditions 
0) —4<Rev <š, ReA>-1, 0 < < 1 and Re (4v+2A +8) > 0, 
(i) y**"g(y) is integrable in any finite internal of t and 


g(y) = O(y*), for İy | — co 
= O(y*), for İy l-50 


1 
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where Re(Zv+a+A+ 8) > 0, 
(ui) din) exists. | 
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ANHARMONIC OSCILLATIONS OF A PARTICULAR MODEL 


By 
R. 8. Kusuwawa, Allahabad, India 


(Communicated by Prof. A. C. Banerji— Received January 5, 1952) 
Summary 


The anharmonio radial oscillations of a stellar model consisting of a point mass, equal to one third of the 
total mass of the star, at the centre and the rest homogen-ously distributed throughout the star, have been 
considered. The skewness and the period increase is calculated. Also the effects of the variation of the 
serni-amplitude and the variation of the effective ratio of specific heats on the skewness and the period 


are studied. d 


1. The theory of anharmonic oscillations was pub forward by Rosseland (1943) 
and was well discussed by Bhatnagar and Kothari (1944). These oscillations for different 
models were studied only during the last three years. Schwarzschild (1949), C. Prasad 
(1949) and L. D. Chatterji (1951) investigated most of the models considering the effects 
of the higher modes also. The author (1951) studied the Roche-model (having an 
extreme concentration of mass) and also considered (1951) the effect of Increasing 
semi-amplitude on the period and skewness for that model. In the present nole the 
author has investigated a model for which one third of the whole mass of the star is 
concentrated in a point at the centre and the remaining two thirds is distributed 
throughout the star homogeneously. This gives a lesser central concentration of mass 
than the Roche-model. 


9. The pulsation equation for the small adiabatic radial oscillations with finite 
amplitude is 


Poron = — (1 +7)’ ə İP,(1 5)?" (1-9 7,9)" ] ^ goes (1 tal: (1) 
0 


where P,, p, and r, are the equilibrium values of pressure, density and the distance 
from the centre, y the effective ratio of specific heats, the relative displacement àr/7,, 
and dashes denote differentiation with respect to r, and dots with respect to time. 

For infinitesimally small oscillations the above equation can be solved (1948) 
giving a set of harmonic solutions in terms of a finite series, with periods 2m/oc& and 
amplitudes nx'r,). We express the present displacement 7 in terms of an infinite series of 
these characteristic functions and unknown time coefficients. We shall take these time 
factors as known in terms of a Fourier series and the Fourier coefficients will ultimately 
be determined. 

Thus we have 


1 = F mlrdga). (2) 


Normalısing n,’s to unily at the surface of the star we get the relative displacement al the 
surface as 
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"n = > q(t) 


(8) 


We substitute this value of q from (2) in (1) and expand the right side in powers 


of are neglecting those greater than two. Due to the orthogonal property of 7, 8 iE we 


multiply the resulting equation by T,'n, and integrate for the whole star we get 





d? 
m + Bide = —— | Du? 22 Di ZAI 


where 
T = gil, 
Br = OO, 
R 
he f Poto Ny. dro, 
0 
and 


1 4 Ro, R z 
Dat = — 5(8- =) (By +1) f Porno mdr, * 3(8y -1) f Pyrolnen ym 
A ü 


R 
Tone naogg)dro + aly +1) f Prog ni dr, 
0 


8. Putting z = r,/E we have for this model 


Mr) = 2. p, (1-227), 


2; 1 
Ja = £ Gp R(À +22), 


P, = x Gp dE el 


Normalised o, and n, are 


nod | 
kä 


-. 9 (q+ 8)(q + 5)—« 
where 
= /(l+e)-1 and «-8- f. 
Y 


Substituting these values of P, and g, trom (10) and (11) in (4) we get 


1 9 
FL Iz vadi +22 Dinar, |; 





where 


1 
Iu | zömda, 
0 


(4) 


(7) 


(8) 


(12) 


(18) 


(14) 
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1 1 . 
Dig = 5 (3- *) @y+n / (1 + Zeletdaide + Y(8y—1) f (1— z* z^ (gat 
Y 5 0 
£ , dé , 1 d ÁF D 
+ nk + yda + Ay +1) J (L—2x*je*nnyede, (16) 
0 
and 
F, = q(q + 5) 424. (16) 
Hence vve get 
Bel l 
a, = (q*89)g*8) 2» f” SEN 
š q(q +5) +2e 
Be... (18) 
` 2q+5 
I = (ran | 1 _ 6 _ 86(q + 4)° } (19) 
` 9 2445 (q+8(q+5)—% §(q+8)(q+ 5) — «P (2g +11) 


A For x 06 we evaluate the coefficients. Consequently equations (18) 
reduce Lo: 








Ga, 2:8727q,"— 0'8127949,4, 
T 

2 (20) 

EE + 10:8028q, = — 0 326024q,? + 48°1602q,q,. 

Q 

Here we subatitute 
q, = > a, COS (! nr) 
2 | 
| (21) 


q, sss 2 b, cos {r.nr) 


where 2z/n is the period of oscillation, 


Equating the coefficients of cos (knr) on both sides of (20) we get infinite number of 
infinite equations. We solve them approximately as explained in one of the author's 
previous papers (1951) choosing one arbitrary coefficient a, = 0'05 which gives the 
relative mean surface amplitude of oscillation nearly 1/20 of the radius of the star. 


The solutions of (20) are as follows 


q, = 0'003648 + 0°05 cos nr — 0001218 cos 2nr + 07000022 cos dnr... | 
(22) 
q, = — 0" 000042 — 0:000099 cos nr — 0'000069 cos 2nr — 0:000084 cos Inr... ( 
where 
n? = 0°98252. (28) 


Thus we see that by the effect of Lbe first mode the amplitude of the fundamental 


16 decreased at the surface. ‘lhe total surface displacement for the star is 
5 —1840P—2 


78 R. 8. KUSHWAHA 


ge = 0°008C06 + 0°049971 cos nr — 0°001287 cos 2nr — 070000118 cos Bnr... (24) 


The skewness for the time-velocity curve is found to be 1.133 for this model and 
the period increase 18 0°80%. The table 1 shows the values for Roche-model and the 
present one for comparision, where z, and z, denote the external and internal ampiitudes 
rerpeotively and semi-amplitude is the numerical mean of x, and aş. 


TABLE 1 












Bemi- | Period 
©; Zs amplitude | Skewness | increase 


— A | YV a TAAA ee À 


0°0523 00476 004990 1'133 0 89 


NEE 


0:0444 0°0507 1'086 









From the above table it appears that the concentration of mass at the centre of the 
star has little effect on the skewness and the period increase. ‘he effect on the two 
quantities is in the opposite directions, i.e. if one increases the other decreases. The 
semi-amplitude in both the cases remains almost the same, although the external and 
internal amplitudes differ. 


8. Now we calculate the coefficients in (18) for x = 0'4 and x = 0'2 and compare 
the results to study the effect of the change of adiabatic exponent cn the skewness 
and the period. 


(i) For a = 0'4 equations (18) become 


d'q 
— LO, = 2/716018,* —0:406854q,q, 


dr? 
: (25) 
SE +15'846g, = —0°480844q,? + 60°507q,q, 
+ 
Solutions of this are 
q, = 0°008428 + 0'05 cos nr — 00011495 ccs 2nr + O'CCCOLES cos Bnr... 
(26) 
q, = — 0000058 — 0°000022 cos nr -- 00000574 cos 2nr — O'CCCCÜEB cos Bnr... 
where 
n? = (0798448 (27) 


The displacement at the surface of the star is 


1p = 0°00889 +0°049978 cos nr —0°001201 cos 2nr + 0°0000116 cos 8nr.... (28) 
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(u) Fora = 0:2 the equations (18) take the form 


P: 
al +q, = 2 59187q,’—0 535991q,q, | 
: 
š po zu" (29) 
2 
d qa 4.28 723q, = — 0:577768q,“ + 99 7566q:q, |. | 
Solving these equations we get | 
q, = 0°008269 + 0°05 aos nr — 0001095 cos 2n7 + 0'000018 cos 8n7... I | 
(30) 
q, = — 0000027 — 0000018 cos n7 — 0000051 cos 2n« — 0:0000028 cos Bnr... | 
v here ; 
n? = 0 98587. (91) 
The displacement for the stellar surface is given by 
ge = 0 008242 + 0°0499866 cos nr — 0:001126 cos 2n7 + 0000018 cos Bnr. (32) 


Ihe values for skewness and period for these two cases are also obtained from the 
respective time-velocity curves. The results are collected in Tabie 2 for comparative 
study. From there it seems that by ıncreasing the value of y the period and the skewness 
both increase although the increments are hardly appreciable. Hence obviously by 
varying the value of y between the physically possible limits the skewness and the period 
cannot be brought in agreement with the observed values. 


TABLE 2 


bemi- Period 
m a mphtudel skewness | ıncrease 
= (z +z.) 


0:05228 0:04764 0°01996 


0:05218 004780 0:04999 


0:05212 0' 04789 0'05000 





6. Now we shall vary the semi-amplitude keepmg y = $ and study the effect 
of this variation on the skewness and period. The range of variation js chosen trom 
ink to eh, R being the radius of the star. For the higher values our approximations 
become very crude. The results are collected in Table 3, in which the corresponding 
values for the Roche-model have also been given from the previous paper of the 
uthor (1951). 


~ 
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TABLE 3 


Semi-amplitudes 


Period-ine1ea8e 


2 | 





Present Roche- | Present 
model model model 





001000 | 0:010000 1028 
- 
0:019997 | 002005 | 1:05 








0029995 0 0802 





0089984 00404 
0049967 00507 


0:0618 


0070866 


— — €— 


m 0°08981 | 0:0948 
| 
| 


009974 0°106 





In the end the author records his most respectiul thanks to Prof. A. C. Banerji 
for his very kind guidance in the preparation of this paper. 
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STRESSES IN AN INFINITE STRIP DUE TO A NUCLEUS OF 
THERMO-ELASTIC STRAIN INSIDE IT 


By 


GuNADHAR Paria, Midnapore, | 
(Received — AÁppril, 4, 1952) 


Introduotion— The stresses in an infinite strip of a thin plate of isotropic elastic 
material due to various distribution of loads have been found by Filon (1908), Melan 
(1926), Howland (1929) and others. The direct method of slress determination given by 
Ben (1938) has been utilized by the same author (1942) to find the stresses when an 
isolated couple acts at a point midway between the edges of the infimte strip. In the 
present paper, a nucleus of Thermo-elastic strain, instead of an isolated couple, is 
assumed to exist ; the nucleus ıs such as would be produced by heating to a certain 
temperature an elemeut of area of the plate while the remaining part of the plate is kept 
at zero temparature. 

Solution of the Problem.—Let the axes ol z and y be in a plane parallel to the faces 
of the strip of which the edges: are given by z +a. If an element of area dQ surrounding 
the origin be heated to the temperature T in the mannor as described above, the average 
stressesin a thin plate evtending in all directions are given by, (Goodier, 1937, Ben, 1950). 


— HaTdÜ wii 
qu = 





Të (ey! 
HaTdO , x3—u? 
ür” S 2 (1) 


2: (2° 3)? 


~ _ EeTdQ zy. 


MET, (ey) 


where E is the Young's modulus and z, the coefficient of linear expansion. . 





Putting 
B2TdQ 
q =— (2) 
7T 
we have, on the edge z = a "T 2. qa 
50 
(a? + y*)* 
(3) 
(ry ena — Mr 
and on the edge z — —a - “sn. 
(pq q — qr la zi 
(aè? y*y | 
(4) 


ay | 


(zU le a = Almen, 
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To nullify the stresses (8) and (4) on tbe edges, we are to superpose a siress system 


—N  — 


LL», TA xy, such that on the edge z = a 


(ez Ju ə o. = ny) | 


(a? ra 
t (5) 
.. _ ay _ 
(rä lee = 20 (a? +)" fal y) | 
and on the edge x = —a 
Lee a = fh 
E ; (6) 
Lite, az —faly) 


The superposed stresses must satisfy the equations of equilibrium and the consis- 
tency relations, These are given by, (Sen, 1938) 








" 
quu. Ae E +9,(8,%) 
Jj. Oat 22 2 ele H 


— 


o 
TH, = SC + 9, (z, y) ) 


where (9,522, + ZU, is a plane harmonic function and 9,(z, y) and e,(z,y) are also plane 
harmonic functions satisfying the relatinns 


öz or OY 








əy oy öz 


— -a m 
— 
go 
— 


We assume that 


qız, y) = J Am co h mz cos my dm | 


° | (9) 
(T, y) = f Bm einh mz sin my dm 
0 


where Am and B, are function of m. 


Substituting the values of 9, and e, in the relations (8) 


we get ag a 
i = m(A,,+ Bm) sinh mz cos my dm | 
0 


\ (10) 
2 — il m(A,,+ Bm) cosh mz sin my dm | 
0 


[^a] 
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irom which we have 
®©, = 2 f (Am + Bm) cosh mz cos my dm (11) 
0 
Substituting the values of ©,, ?, and e, in we get 


zz, = f (Aal cosh me- mge sinh mz)—B,, mz sinh mz]cos my dm] | 
(12) 
zy, = f [Am mz cosh mz + B4(mz cosh mz) +sınh mz)|sin my dm 1 
0 


From the boundary conditions (5) and (6), we have, by tbe help of relations (12). 


Js İ Fen) eos my dm = AN | 
t (18) 
ala İ Tum F ,(m) sin my dm = f (y) | 


- 


where (2/z)* Ë (m) = Am (cosh ma — ma sinh ma) — Bu ma sinh ma | 
(14) 
(2/z)* F,(n) = A4, ma cosh ma + B,,(ma cosh ma + sinh ma) 
By the ‘transform theorems’, we have from (13), 
I' (m) = m" fi(y) cos my dy = — Q(2/a)* dE [T ; Cos my dy | 
L (15) 


= (2/m)! zm sin my ay = 200 (cedes ce SCHO | 


Evaluating infinite integrals ın (15) and substituting the values of F,(m) and Fm) 
from (14), we get 


A m (cosh ma— ma sinh ma) — Bm ma sinh ma = Gate mg | 
(16) 
A, ma cosh ma + B,,(ma cosh ma + sinh ma) = Qine ma 
These results are satisfied if 
A = 2Qm(ma+e”"@ sinh ma) | 
m 2ma +sınh 2ma | 
(17) 


2Qm(e ma cosh ma — ma) | 


2ma + sinh 2ma ] 


Thus, from (11) and (12) 


£ 2: m cosh mz cos my dm 
” 2am + sinh 2ma - 


er, = 29 f ` mooemy _ [(ma + em sinh ma) cosh ma — mz sinh mz dm 
2am + sinh Zma 
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m sin mı : m f | 
20 = 90 Gar nah ma sin TJ mz cosh mz + leng cosh ma — ma)sinh mzjdm 
2ma + sinh 2ma 


— — m COS m ' 
yy=©,-12,=20 f LESEN, cosh mz + mx sinh mz—cosh mz 
Oma + sinh 2ma | 


i (ma+e-"4 sinh ma) ]dm ) 


The resultant stresses can now be written as | 


mee [MA 


zr = TT, + zz, 


yy=YyıtyYs E 
: gy = TY, iy. 


m cosh mz cs m cosh mz c^s my 
= T2 + 
2 = yy = = f  9ma-&inh2ma m | 


where the two component stress systems are given by (1) and (19) — 


Since zz = 0 onthe edge 2 = a, we have 


eo | 
(YWaua = (O)asa = 4 f M CORE NAGO LEM NE 
Í Dam -+ sınh,2ma 


Putting ma = u and y = an, we get, on substiluting the values of Q from 12). 


2BaT40 


(YYy)asa = nit) 


where | 
_ (^ ucosh u eds un 1 
.: / Zu + einhiZu ` 
When n 1s small, ! 


l 
i 


i “J mes wcoshu gy m u° cosh u ER 
Qu + sinb 2u > 2u + sinh 2u 


Tbe values of the infinite integrals, as given " Filon (1903), are 


f u cosl u Ju = 0018. f wu? cosh u du = 5'750 


: 2u + erh Qu A Du -+ sinh Qu 


Hence, when n 18 small 
I(g) ~ 0'918 — 2°875n* 


= l 
similarly, {he value yy on the other edge x = —a can be calculated. 


(19) 


(20) 


(22) 


(28) 


(24) 


(25) 


In conclusion, I cffer my grateful thanks to Dr. B. Sen for his .. help in the 


preparation of this paper. | 


MIDNAPORE COLLLEGE, 
MIDNAPORE, WEST BENGAL. 
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ON THE EFFECT OF A SMALL SPHERICAL CAVITY IN A SEMI- 
INFINITE ELASTIC SOLID UNDER STRESSES PRODUCED 
BY A COUPLE ON THE PLANE BOUNDARY 


Br 


SISIR CHANDRA Das, Chandernagore 
(Recewed —April, 21, 1952) 


Introduction: Dipolar co-ordinates have recently been used by Sternberg and 
Sadowsky (1952) in solving some three dimensional problems of elasticity. In the 
problems discussed by the authors the symmetrical displacement in which the azimuthal 
component was absent was only considered. In the present paper the solution of a three 
dimensional problem involving the azimuthal component only has been obtained. 


Dipolar co-ordinates have been used here to find the concentration of stresses due 
to a spherical cavity situated near the plane face of a semi-infinfte solid when the stresses 


have been produced by a couple acting at the plane surface with its axis perpendicular 
to it. 


Nomenclature : 

x, y, £ = Cartesian co-ordinates. 

0, 9, z = Cylindrical co-ordinates. 

£, n, ? = Dipolar co-ordinates. 
a =A positive length. 

h,, h,, h, = Stretch ratios. 

r = Radius of the sphere. 
d = Distance of the centre of the sphere from the plane surface. 

u. v, w = Components of displacements in dipolar co-ordinates.- 
A = Dilatation. 

2a¢, Ze, Ze, = Components of rotation in dipolar co-ordinates. 


Tt, Ca, To = Normal stress-components in dipolar co-ordinates. 


g? 
Ce. Gə, Ges = Shearing stress-components in dipolar co-ordinates. 


M = Moment of the couple. 


Ay, Bn, Am Bn = Constants, 
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1. The Coordinates. 
We take the curvilmear co-ordinates defined by 


ZEN - 0 + ilz "E a) | : 1 1 
£ in — log od ig — a) ” 4 ( = ) 
where g = (a? + yt, x, y, £ the cartesian co-ordinates and a, a positive length. 
If o, 9, z be the cylindrical co-ordinates, we get 
_ 4 Sin n cos 9 
cosh £ — cos ?) 
— asin yn sin 0 
2000000050 (1.2) 


cosh £ — cos y 


and | 








The relations 





č 
H- uc - 7 (=) 1.8 
dy (a tn +a 57 
and i. _ ()+(2)- (2) 
h,? (Nöş Oo Öp 
yield 
1 
hehe x (co:h £ — cos n) | 
and (1.4) 
1 1 
h. = ` = 
e bü” (cosh £ — cos n) 


The surfaces é = constants give a serios of coaxial spheres having a common radical 
plane £ = 0 (or g = 0). We take the plane surface to be given by C = 0 and £ is taken 
positive towards the interior of the body. Now u sphere of radius r having its centre at a 
distance d from the plane along the s-axis can be obtained by properly adjusting a and a 


by the relations 


T = acosech a I 


and (1.5) 


d=acotha | 
£ = =, giving the spherical surface. 
3. The equations of equilibrium and the required solutions. 


Let u, v, w represent the displacements ue, u,, u,. Here we take u and v to be 
zero and w independant of ¢. 
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Then 5” D ane SUS 9d See ee 


(2.1) 
10 
“um ee): | | 
and | 
40 mU vus xax deg 
Two of the equations of equlibrium (Love, 1944) are ıdentically satisfied and the third 
takes the form 
OE Uo ag Ou te ön 
The solution of this equation can be written as 


= se x (As cosh (n+4)¢+ B, sinh (n+ AEN P (cos q) — Pax (cos n)] (2.3) 
ji wa) 


when Pa is a Legendre's function of the n-th degree, and An, B, are constants. 


Since n(n+1) 


(u.27—D)P a(n.) = x 





(Pry ii) — Pha) (2.4) 


we can take as a solution suitable for our problem 


(cosh £ — cos n)* ena 
a 


Y {An cosh (n + 3)£ + By sinh (n+ WEP (2.5) 


ne2 
PL denoting differentiation with respect to cos y. The stresses that do not vanish are 
given by 


məni SEN 
(reg) EE ő 


= — [S(cosh £ — cos 1)* sinh £ E {An cosh (n +4)£ + Bn sinh (n+ HEP» 


10, = 


+ (cosh ¿ — cos 1)? /? > (n+4){A,, sinh (n + 3)£ + Bn cosh (n + HE, 
and SS (2.6) 
= = 9 fo; 
(74), = H EMN Ən 2a) 


il 


-2 


7 7 Tö(ecsh £ — cos n)! í? > {An cosh (n +4)€ + Ba sinh (n + HELP» 
a 3 


— (cosh £ — cos n)5 12 s cosh (n + 4)£+ B, sinh (n+ š) CU, 
p, being the mgidity of the solid and P, is the second derivative of P, with respect 
to cos 7. 
If a couple of moment M with its axis normal to the plane surface act at the origin 
on this plane, we have the displacemer t due to it in the form (Love, loc. cil.) 
9—1840P—9 
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u = =0 e 
and M 


(Q1 — "7 


0 
iu rn (2.7) 


Alternatively, we can put 


— M sin (cosh £— cos d 2/2 > (= 1)?e - HOC P., (2.8) 
Arpa" = 


From this result we get 


My sin a oo aie J 
bah = an L(eosh € —eos q) P 3 (UH (n Be t BEP, 





+8 (cosh £ — cos ail sinh £ $ (— 1)*e- ** CP; ] 
nel 
and (2.9) 


an ” 
(54), es Misin"? I [ — (cosh £ — cos m5 /* 7 ( - 1)^e- Ə PEP,, 
“HW Onpa® pa n=1 


+ $(cosh É — cos mi IC 1)^e- (-DEP/] 


Let {he complete displacement be given by 
w= w, + 10, (2.10) 
Now w, is to be so chosen as to make 
(rg); = 0, on £ = 0 (2.11) 
(Teo) m İrəli = 0, on £ = & 
where suffixes 1 and 2 denote stresses calculated from w, and wy. 


On using the recurrence formulae 
! 


| PL) — u Palu) = m +1)P,(u, 
and . “Üz ) K Hi ( ) H ) l (3.12) 


| Pre) Philu) = Pals) N 
we find that the above conditions are satisfied if we take 
e B, = 0 
_ Mv 26-52 


| S ua sinh ğa 
i 4 sinh Aa A, = [3 sinh a cosh ğz + 5 cosh e sinh Ja] A, 


-VAM e- [2446-7 [2) (2.13) 
mpd 
and when n > 2 


[(n + 2) sinh (n+ 8«]4,,, = [B smb o cosh (n + 4) | 
+ (2n + l) cosh 2 sgmh(n+4)a]A, +(n—1)A, ., sinh (n — })« | 
—2M,[(n —1)e “710 + (n + 2)e- (t3 1996) 
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93 
where 
(—1)"M 
SS . 2. 
"VQ. apa (2.14) 
The complete stress at € = « is obtained from (2.5), (2.9) and (2.14) as follows: 
(r4 tə ES bir. ER 
14542 o 
= — (cosh a — cos m” 2, LA, cosh (n+ Ain — Mpe- tHe] 
x (8P,—(cosh a-cosn)P»}. (2.15) 
where A, = 0. 


The series evidently converges. 
When « = 2, we obtain on retaining first few terms 
ma’ 


M 73 Tw sini (277622 — cos n)* x [0°0B784 + 0°0083(7'5 cos a — 11'2866) 
+ 0°0366(26°25 cos??? — 56:488 ccs q —2:25)1 (2.16) 
approximately. 
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ON CERTAIN FORMULAE IN OPERATIONAL CALCULUS 


By 
N. N. Boss, Lucknow 


(Communicated by Dr. S.C Mitra—Received May 6, 1952) 


The object of the present note is to derive certain formulae in operational calculus 
and to obtain certain contour integrals for Bessel Functions and Weber’s Parabolic 
Cylinder Functions. 


1. We have the formulae: 


If e(p) = f(t); Yp) = g(t), 
then , 


J LU (9 a, = J fn © dt (Goldstein, 1982). (1) 
Let š : 
y(p) = Er then g(t) = exp (— at). 


The relation (1) becomes 


f(t) HE d am 9(z) 
dt = j exp ( —az) P) da, (2) 


o t t 
Let us write 1? for t, a? for a in (2) and multiply both sides by J,(ab)a'"*!, Integrating 
with respect to a between the limits zero and infinity, we have i 


qm m SS 
2 | "yd | Zelt ` dz J s) de J exp (—a®x)J,,(ab)a™*da- 
0 0 


where we hara — that change in the order of integration is permissible. 


Putting b = p and z = 1/95, we have on integration 


İ ims = s İ v?" *1 oxp (— 1p?v?)g(1/v?)dv 
0 


0 





= dus f exp (—ip*v)s(1/v)dv. (2) 
0 
If we put m = — $, we fall back upon the result (Humbert, 1984, p. 12). 


f p (—pt)f(t?)dt = 1f exp ( — 2p?v*)9(1/ v?)dv. 
0 


0 
Again writing 4p for p?, we find that 


zie. J emt piat + rend 
0 
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Bo that if ç(p) + f(t), then 


ap | fX, (pb)4t «0/9. 4) 
0 


2. Let us now write 2p for p and n—4 for m in (8). We obtain 


à [ eti” KİL (opa)da = pi | v*-1 exp (— p'v)e(L/v)dv 
0 9 
Let us multiply both sides by p""texp (—3p")D,,(p), and integrate between the limits 
zero and infinity. We have on changing the order of integration (assuming that is 
permissible) 


4 J a" t if(e*)da J p^- t exp (— 1p?) K, -,(2pg) D, (p)dp 
0 0 


=f vr-to(l/v)dv fə exp [ -2*(v +) |Dan(p)dp 
0 0 
On using a result due to Mitra (1046, p. 171) the left side becomes 


elt — LHT (9n) il gün exp (s3)/(s?)D....(2x)da 
0 
and the right hand side is 


Tan) f vk (m+, n; i; —2v)o(1/v)dv. 
0 
Hence we have 


| am exp (a°)f(2*)D A(2a)da 
o 
Le | ont F (nk, n; 4; —2v)9(1/v)dv (5) 


0 
Again let us put m = n+4 in (8). Multiplying both sides by p"+r3/2? exp (— 1p?) D4.,.(p) 
and integrate with respect to p between the limits zero and infinity, we get on using a 
result due to Mitra (1946, p. 172) and simplifying 


— 
— 








Hə exp (8*)f(2?) DL, ,(29)ds = ur f ro, n+2; 8; —2v)o(1/v)dv. (5) 
0 ü 


8.1. We shall now give some simple applications of the above formulae. 


In (8) let us put m = 0, p = 2c cosht and integrate with respect to t between the 
limils zero and infinity. We get 


fJ Eear = à [exp Cei ete ée (6) 
0 0 


Let 
(p) = p" exp (3c?/p); f(t) = (2t/c?)#*I,(c  (2t)) 
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we get 





f mi (V 26s)E, (call = 27 (Dn DF, Bln) > —1. (7) 


: C 
8.3. We have the formula (Humbert 1984, p. 12) 


v —” a e(p - 1/p) 
Í J. [27 fa(t—at|f(xjda + “0 
Appiying this to (4) vve get 
x i 
af (GK [2vilp+1ip)lat = [ n visti iTe( add. 
0 


0 


f(t) = dell (at), (p) = exp (-a/p). 


Let 
We get 
ob i 
2 J J (27 (at))E (2 {t(p +1/p)hdt = d exp (—az)J,[2 {z(t — z)]]dz. 
0 0 


The left side is p/(p?+ap+1). 


We therefore have 
1 
il tJ 120 y 1261 s)11 exp(—atz)dz = p/(p? 4 ap+1). 
0 


Let a e2. We know that p/(p 1)” + texp(-t); whence by  Lerch's theorem 
(Lerch, 1908) 


1 
J exp (-2tz)J [2tyjell-z)]de = exp (-t) (8) 
0 


f(t) = t-t na 2, (2t); o(p) = (zp)! exp (~2/p). 


Let 
We obtain 


2 fe cos 2y (20) K [2 v ft(p + 1/p)l]dt = rt f ate?) [2 {z(t — z)]]da. 
0 0 


The left hand side ıs 
$ 
ozü = 4i(D* liy (Qt) — DA.(- i (2t)) (Mitra, 1988). 
Hence again by Lerch’s Theorem, we have 


0 
More generally, when 


f(t) = tJ,(Q(20), olp) = 2p" exp (—2/p), 


1 
J 2-t exp (—2tz)J [21 let) -2,}]de = dat [D (iv 2t))—D2(-iv (2t)] (8°) 


we have 





1 
J z"J [2t V (z(1—2)]] exp (— 2ta)dz = 
0 


15777 n+2; —t) (8°) 
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8.3. It is known that 


avtt oxp (42?)D_2,-n(x) and z"-+ exp (1:”)D al are R,. 


J unn exp (1z*)J4(zy) D -s4-o(2)d2- = y^ exp (1y*)D - 24 -a() 


s | 
Now | 


I 1 
—f} -4n = 
Writing $2 for 2 and putting | 


f(z) = y *s-vJ,(2y 2) and 9(p) = 


That is 


Fi(1;n+1; —y?/p). 


Fin DS rp: F1; nl; —y* p) 


we get 


f sra ntl; —y%w),F n+, n+2; 8; —2v)dv 
f = (—1)*trtI\(n +1) exp (ty?) D_an—sly). (9) 


In a similar manner, with the same values of f(s) and 9(p), we get 
found: n+1; NALE n; k; —2v)dv 
i = (—1)PrtT (n 4 Dy7 exp Gy1D. ag). (10) 
X. Certain Contour Integrals. | 


If 
- In in 
= ii 


with £ real and positive, it follows from an application of the Bore! sum formula that ' 
| 


g(u) = uf exp (—ut) f(t) dt. (11) 
0 
If the integral (11) is absolutely convergent when u is real and equal to x, then p(u)/u 


is an analytic function provided R(A) > « and we have 


since Ä 

IO ne | p NEN IUE 
ë | 2” : exp ( —az) zo de, 
. we get by an applieation of this theorem (Mellın, 1902) 


9(a) 1 
ü Ont Eé 


exp (aA) əf faal, (c > 0) (12) 


&.1. Let f(t) = t^* cos bit. Then p(p) = (rp exp ( — 25*/p). 
Therefore e(a) 63300 


eti 
1 


=o exp (aA —bAY) ACTIA 


ço € — t @ 
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Substituting for p(a), we get 
Ë 1 é + 100 
( = "exp (- göra) = 2 | exp (ad—DAHA-HAA 
e — (OD 
Let a = 3. Therefore 


(2:)F exp (— 45?) = ij exp (4A — bA3)A FAX 


o- 27) 


Differentiating n times, we obtaın 


ctro 
Det asa pa f exp (4A — bak)Jarn-Dax. (18) 
2 (2) ue 
The path of integration is curved (if necessary) with cuts so that the integrand is a 


single-valued analytic function. 


4.9. Let us multiply both sides of (18) by exp (—#b*)b™ and integrate with respect 
to b between the |: mits zero and infinity. 


Inverting the order of integration,* we get 
[ex (—#b?)b™D,,‘b)db = raai J exp aan f exp (— 357 — bAt)b™db 
0 24 (2n): e—10 0 


Lm) id Ata- And 
= a] exp (3A) - m-1(A)dA. 


But the left band side is equal to 


miğin-m-iT(m--1) 


“ann (b d perdo gm nn 


Hence 
1 0 + í 
5x. exp (3A) 1(A*) 
Odn-m—1 
= sch of us +3; —3n--l, $). 14 
4.3. Again 
1 -. > Ad(mtn-0gm 
D (b D, | A—bAt eqjm--- E 
N (-1)” ` m 15 
- 5 D,(bj)Dman(b)c™, (15) 


on term by term integration which can be easily justified. 


* Inversion can be easily justified. 
8—1840P—3J3 
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By applying (1) we can easily prove that if 
e(p) = f(t), then p/(p+z)t = exp (—12)/ y (xt) and so 


"f(z)de 1 


120770 = — —1t)0-/5f3) 
; a m 77 to)v e(v)dv. 


oe, 


Hence if we write 





then 


o+ia 
tot) = .. E. ` e mal Wd (16) 


f(z) = 2-tK,(Qaxt), 9(t) = at: exp (3a? /[t) K, (a? 29). 
Then (16) reduces to 


1 eko E (az) 
ges | trA)d J 7 
ül a... 


Let 


m /2 064-100 i 
E | (J (aM) +y (aX)) exp (tA)dA. 
Hence we have J 
a? erim I 
t^ exp (4a7/t)K (5 ;) = = rd I (Ja (GA) + y,  (aA)) exp (tA)AA. (17) 


We may remark here that the condition for the. validity of the changed order of 
integrations is the same as that for integration under the integral sign when one of the 
limits is infinite, since the processes are similar. When both the limits are infinite, the 
conditions are the same as those given in Hobson’s Theory of Functions of a Real 
variable and the Theory of Fourier series, (Hobson, 1907). 


My best thanks are due to Dr. 8. C. Mitra for his help and guidance in the 
preparation of this paper. 
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GENERALIZATIONS OF STIELTJES TRANSFORM 


BY 
K. M. SAKSENA, Sagar, M. P. 


(Communicated by the Secretary — Received May 22, 1952) 


1. If we iterate the Laplace transform 1.0., if we take 


a» 


fe) = f eolia 


where : = 
o£) = J eoa 
0 
then 
— İ YÜ) 
f(z) = ef (1.1) 


and the last equation is referred to as Stieltjes transform. 
If we replace Y(tjdt by da(t) we get the more general case of the equation (1.1) in 
the form 


_ f“ dalt) 
NT 


Dr. R. S. Varma (1952) has given a generalization of Laplace integral in the form 


ao 


f(s) = | (sme i YY) (atda) (1.9) 
0 
where Ve, (s) denotes Whittaker funclion, and in the theory cf this generalized Laplace 
integral he gives two forms for the an alogue of Stieltjes transform. In either of these he 
takes a(t) to be absolutely continuous so that the generalised Laplace integral (1.2) can be 
written as 


f(s) = ik (at) tem, (st)o(tjde (1.8) 
0 
If f(8) be the generalized Laplace transform of 9(u) in the sense of (1.8) and o(u) be 
the ordinary Laplace transform of y(t) then 
I'(2m-r1) T İs: 1 
alm — k + 8) R S CIm-kae8^ 


where R(2m+1)>0. This is the first analogue of Stieltjes transform. 


t 
f(s) = | ynat (1.4) 


If f(e) be the ordinary Laplace transform of ef) and 9(u) be the generalized Laplace 
transform of W(t), then 


_ Tem+1) [^1 2m cl, 1. a 
(8) meer | mung mı Pos o 


where R(2m + 1) > 0, 
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This is the second analogue. 


In this paper I have given a few more generalizations of the Stieltjes transform by 
extending Dr.. Varma’s results. | di 


2. The first generalization is obtained by iterating the generalized Laplace trans- 
form (1.8) and the result can be stated 1n the form ol tlıe followıng theorem. 


Theorem 4. I | 
f(a) =| (sur tet], „(su)plu)du (2.1) 
and ` j a i 
sa) = f (ut7o HH, ,tul)y(i)di (2.2) 
0 


"ww Läit + Dm Han D) f pma 


L'(à—k-—m)I'(2m-n-—l4&) 
x F [ie mer. “üsu 2 ; ad 
1(2m)T(2n +1) $ = 5 1, 2n «1 Ed 

2. 0000900 0. t)dt (2.8 

I4—k+m)I\(n—1+ 3) 3 311—2m,n I+ : y) ‘ 
provided that Hiel > 0 R(p+ntn+1)> S R(m+*m +n +x3A3: 1) > 0 and 2m Æ T 
and W(t) = Olte) for small t, 

= O(e t”) for large t, R(v) > 0. 
Proof ı Substituting the value of ş(u) in (2.1) we have 


Hs) il (sum "fer: Are Wir, q (gu) f (ut^ ke AW, (ut) (tdt du 
0 0 | 


= ami í in—tu(t)dt) d umh-1g- bu, (su) Ws (ut)du (2.4) 
0 0 | 
on changing the order of integration. 
Substituting the following expression, in terms of Kummer’s function, for Wx, (su), 


Wr, m (su) ee ‚F[ld+m+k;2m+1; -su] 


I(—2m) a Sa | 
m0 (su) ++ F[4-m+k;1-2n; —8u] 


i 
4 


+ 


2m not being an integer, we have 


= amt L. 7 1) Zen, (jo mk), amin 


= A -m+r+4 
x (ut)2m+n+r—-+ e iut W, (ut) + . I(2m) QG-mak), sten 
ant) I{#-k+m) (1—2m|, intr 


x (ut) t7 eit W, n (ut) ] t du 
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Integrating term by term, which is valid since 


(i) the ,F, series is unitormly convergent in the arbitrary interval (0, A) of u, 

(Gi) (ui)?**»-tg-ivtTY, (ut) and (ut)*7*e-1* W; (ut) are continuous functions of u 
for u > 0 provided that E(m) > 0, R(n) > 0, 

and (ui) the series obtained upon integration is uniformly convergent in (0, R) 
with the help of Goldstein's integral (Goldstein, 1932). 


İ səmi, ea D FP Ch, F,[ mig, m3; l-k1; a 
i A. Dim + 1) > 0, Blo? 4-1) > 0, 
we get 


: I(-2m)  (@tm+h), emt 
= gü-—tk n— + Mc cur Mem 000 o i 22 = 
.. | .. 2c rl Ir (&-k—m) (2m+1),. nn. 


X L(2m +r+DI@m+2n+r+1 T(2m) (A-m+k), amr 


8 —— — ———— ç -.— nn 


I(m+n-Il+r+3) T(y—k+m) (1-2m), intrt+t 


y T (r* DIQn- r1) | 
f I(n-!I+r+3) i 
provided that R(9m-t 2n +1) > 0, R(2m +1) > 0, R(2n +1) > 0. 
In order to justify the change in the order of integration in (2.4) let 
X(t) = yaer- f umth-1g-du6 HOW, (au) Wi, (ut)du 
0 
and 


m | 
Blu) = unto, (su) | tnckecistet m, utat 
0 


where A is small. 
Now x(t) is uniformly convergent in t > 0, if R(e+n+n) > 0 and 
R(mimtin+n4+i1)>0. 


Similarly B(u) also converges uniformly in u > 0 provided hat 2555 € >0 
and R(n +n +o +1) > 0. 
Next let us consider the integral 


r= Í duni, sn) du f eoim mi tut | dt, 
T T’ 
where T and T’ are large mə a 

It is easy to see that I does not exceed a constant multiple of 


il |atultmtntk-1g- u | du f: ]t5-ig-nt-:i | di 
T T C 


which tends to zero provided that R(,) > 0, R(v) > 0. 

Hence the change in the order of integration in justified when R(w)>0, R(s)>0, 
Rin) > 0, R(m) > 0, R(p) > 0. 

By analytic continuation, these conditions can be relaxed to those mentioned in 
the theorem. 
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If, however, we substitute the value of VVLAfut) in terms of Kummer’s function in 
(2.4) and proceed on similar lines, we shull get 


f(s) I (7 2n)I'(2n + 1)T(2m +2n + 1)s-29-1 


ao 


< f Ps [2n 2m +1, nb i mein-ke 2; - = | yat 


n 


I(2n)I(9m--iject (^ , [2m-c1,1,l—n41 t] f 
Breed Pas —9n+1' SIE 
where R(v) > 0, R(s) > 0, R(g--n +n +1) > 0, R(mtm+n+£n+1)>0. 
Corollary 1: If we put k = 4—m in (2.8), we get 
f(8) Tap) ake E -:]wo dt, 
R(8n 4-1) > 0 


This is the same as the result (1.5) given by Dr. Varma. 
— Colloraıy 2: Tf we put l = $ -n in (2.5), we get 


— T(2m-+1) m ave puna. sh 
/(s) sl Ja 1. tim k+3; € |va)as, 


R(2m +1) > 0. ` 
which is the result (1.4). 

3. ‘The second generalization of Stieltjes transform can be obtained by taking fis) 
to be the generalized Laplace transform of e(u) as given by (1.2) and 9(u) to be the 
generalization of Laplace transform of ¥(t)in another form given by Dr. Varma (1947). 
The result can be stated in the following form. 

Theorem 2: If 


fs) = | ("ie (oujetu)du 
| 0 | 
and e ; 
eil = | Qtuy-im,.ai)y (tdt 
0 


then I(-2m) (f+m-k), amir 


BS J 0902 2 nq Tm) (m 7 rl imn 


, Um * nr Tm n1) p hə 2m—n+r+h i-z] 
I'(Zm-l-r- 1) Tus  2m-I+r+4 2t 


. rm) ` 4-m-k), srtH Trt nt Ql (r—n + 5) 
° + 2T bem) (1-2m),rl (209164 Dien" 


xa F rens r-n+3;r-I+]7; EA) (8.1) 
where Reexn+3)>0, Rím +m +n + $) > 0, Ris) > 0, RU) > 0, 
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and the behaviour of y (t) is given by 
y(t) = OliP) for small t 
O(e-!") for large t, R(v) > 0. 


Proof: Substituting the value of ọ(u) in f(s), we have 


f(8) = il (su)™-ke- tT, „(su) J / (Qtu)-# 07; ët V(t)dt du 
0 0 


= gn-i f GIN di d vn Ze Ben, (att) Wis (Qut) du (8.3) 
0 0 


on changing the order of integration, 
4+m-k ) 
; 8u 


EE MM ii _ (2m) = 
= gm-t 2t)- *y(t at f "| mtt. L 2 um- F ( 
i / -. : ` : Dd-k-m) t NIm +L 


= I'(2m) u | 
+ g7mtt_ + —k—m: — : 
8 FGK T Pilg- k-m; —2m+1; eil W,,,(Qut)du 


on substituting the expression for Wy, m(8u) in terms of Kummer’s functions. 


Writing the equivalent series for ,F, functions and integrating term by term, which 
is easily justified, we get 


f(g) = eni f anal 


0 


a ası Saa ¿SOS t-r] 
Die —1+ r+ 1) ° L2m-l+r+3 S 
T(2m) ame. 00. Pere Sn + 


am (1-2m), nl (2) ` DTr-1-12 


re) 


$ I(—2m) (raum, amir 


L(à—-k—m) (Qm+1), ri (2tymiei 


r= f) 


x Fi ren f, r—-nti;r-lei; i =| 


In order to justify the change in the order of integration in (8.2) let 


a 


X(t) = yt) if unig-tnW, (su, Hin. „(Qut)du 
0 


and 


Glu) = ec iun- Wy, (su) | "VET, ut) g(t)dt 
0 


Now X(t) is uniformly convergent in £ > 0, if 
R(m+m+mn+33)> 0, R(s)>0 and R(e +n + 1) > 0. 


And 6(u) converges uniformly in u > 0 provided that 
Resn+3)>0, Rm +m +n +) > 0. 
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Next if we consider the integral 


I =f um-te-i W, „(8U) du Í 


T T 


t+ Walut) y(t) | dt 











where T and T’ are large, we find that the integral does not exceed a constant multiple of 
| «f 
T T”, 
which tends to sero provided that R(a) > 0, R(v) > 0. 














Hence the change in the order of integration is justified under the conditions 
R(s) > 0, RW) > 0, Ro+n+h>0, R(m +m +n + 1) > 0, 


which can be relaxed to those stated in the thoorem by the principle of analytic 


continuation. Ki 


Oorollary: If we put w = 4-min (8.1) we gel 
“1 +m Dn zl I; Él hd 
= +n, ġ-n; f- — 
fis) | Gr DO 7) 4 -ng i Y. 
which gives another analogue of Stieltjes transform: 
A The third generalization is obtained by taking f(s} to be the generalized Laplace 


transform of ş(t) in the form given by Meijer and Greenwood (1940) and 9(u) to be the 
generalized Laplace transform of Y(t). The result is as follows. 


Theorem 3: If I 
jlo) = j f uii uou) du 
T 0 i 
and | 
. plu) = İ (utn-iecimuzstutyi (dt 
0 | 
then 
I(—2m) | (4t+m—k), I(2mctr-v43)I(2m-rtvt 3), 
ee il $n -m-k) (2m + 1), rte | gam r4- 8/2 
, POR id, Tüm). (mk), .L(-vegI(rrve$, Prip (t/a) ] 
omm iti T(axm-k) (—2m +1), rb grt 2/2 (£? 187 — — {jit 
| x y(t)dt (4.1) 
provided that — E(g-- m +m + 1) > 0, Rim + +v+]) > 0, R(B) > 0, R(s) > 0, 
and š Y(t) = O(t^) for small t, 


= O(e-"*) for large t, R() > 0. 
Proof: The proof is similar to that of theorem II, 
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If we put k = &—m in (4.1) we have 
= a 2 z ` 
f(s) = İTg-əT6-ə) J (5 -1) m utpat 
0 


which reduces to the ordinary Stieltjes transform if we further put v = A 


Corollary : 


I am grateful to Dr. R. S. Varma for his kind guidance and help in the preparation 
of this paper. 
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TWO THEOREMS ON SELF-RECIPROCAL FUNCTIONS AND 
A NEW TRANSFORM 


By 
K. P. BHATNAGAR, Lucknow 


(Communicated by the Secrelary— Received July 7, 1952) 


4. Gupta (1945) has given two theorems which enable us to derive by processes 
of differentiation and integration new self-reciprocal functions from the known ones, 
which are self-reciprocal ın the Hankel Transform. 

The object of this paper is to generalise these two theorems and to obtain by the 
help of these, new self-recıprocal functions from the known ones, which are self- 
reciprocal in the new transforms. In recent papers (Bhatnagar) we have proved in 


detail that the kernel 
d 
m, (zs) = (sa)i f J, (01, (25) 4€ 
0 


plays the role of a transform. We shall call a function Ra» if it is self-reciprocal in the 
new transform for which the kernel is m, (v4). 
It can be easily shown that 
n, Q... ne (za) 


defined as 


“Tao m... w. (62) = (za) ZEN f Tallid alba... da, atn- x da, E 
Ú 0 t, ..... LE" tit. . 
can also be taken as a kernel in the most general transform. In proving our theorems, 


we shall use this kernel and obtain the generalisations of Gupta’s theorems. 

2. Hardy and Titehmarsh (1980) have defined a class of functions A(ə, a), where 
"0 < o =< r, G < t, as fonctions which (i) are analytic functions of z = reif, regular in the 
angle A defined by Tr 550, |6| € o and (ii) are OT |a|-2-8] for small x and O[ |z ]e-1*?] 
for large x, for every positive 6 and uniformly in any angle (gl woneem, Concerning 
this class of functions, they have proved a theorem. In the paper referred to above, we 
have given a similar theorem. Our theorem is 








ur 


A necessary and sufficient condition that a function f(x) of A(w, a) should be BA, 
ig that it should be of the form. 


f(z) = zi f 2T(1- An rhel Dt iv + 48)z 5 2(8)dəs (1.1) 


where (8) is regular in 
gege Lü: (1.2) 


satisfies the identical equation (s) = V(1—5) and is 
O eir əə) tel] (1.8) 
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for every positive n and uniformly in any strip interior to (1.2) and c is any value of o 
in (1.2). 

This theorem can be extended when f(x) is HR, a" the proof being on the same 
lines. 

Theorem I. Let E and E, stand for the operators 


0 

respectively, R(A) lying between a and 1—a. Then E"E,--ET^E," is an operator 
tiunsforming an Ra, a, an, into the Rain... s, fransform*, where E'"E,^ stands for 
the operation on E performed m times and the operation on E, performed n times in all, 
in any order. 


x < 
a1 | z-*dx and a | nide, 
0 


Proof: Let f(z) belong to the class A(w, a) and be Has... an Then by the 
theorem stated in Art. 2. 


o+ 100 
f(x) = J P(8)z *da 
em 160 
where 


P(s) = ATE An, Hd) + an, tds)... E+ dna + do C). 


and Y(8) satisfies the condition 


del = V(1— 8) 
and 
dia? = O [eG - ets) Iti]. 
Now 
Ef(z) = a | səda | P(8)a7tds 
0 e-12 
c€ + LO x ctio . 
— 2a f P(ds f gözdə == f ee (3.1) 
(b 0 6-00 


The change in the order in (8.1) is valid for R(A+8) < 1 by de la Valleo Poussin”a test, 
on account of the absolute convergence of both the æ- and s-integrals; the latter by virtue 
of the asymptotic behaviour 
a +iß+iyt) = O(e” be ivita-Y) ag | t 1-ə co, 
| P(8)e7* | = O(e(- t? It ) for all q > 0. 
Repeating the process it i8 easy to see that in general 
e + to E 
seo (7 peni, 


om im 


Similarly 16 inay be shown that provided Ris — A) € 0, 


E,"f(z) = il ən 


omi 


* That 13 into another Ha, s,... s. transform. 
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16 follows therefore that == 


¿+ v - P(a)z "° - > 
E"E," = 00 —4 
provided R(s—A)«0 and R(e+A) <1, ie. a< R(A) <1-a. the order of operations 
being indifferent Hence 
(E mn + E E) (a) 
otio 
... opaet A dn, + he). G+ Ana +48) x s (8) "de (8.2) 


Zei 
e- ı m 


where 
y (8) x ((1—a—2A) "At + (178 A) "(8 - A) "h(s) 


since V,(8) satisfies all the conditions as ¥(s) in Art. 1, the theorem follows on comparing 
(8.2) with (1.1). 
Cor: When m =n, we have the theorem: 


If f(x) is Raum... n, then B™B,™f(x) is also Ram, na provided a < RO) < L-a. 
Theorem II. Let A= X +a and A, = A-i-2 t; then A™A,"+A"A,™ 13 an 
z z 


operator transforming an a... n, into another Ran, ..n, transform. 


Ihe proof of the theorem follows the same lines as those of the previous theorem. 
It may be noted that 


d d d 
AA, = | =| | — (eae =] — = -( d 
7 ME A—1 er = AA 1) Ta 


since the sum of two or more Fa, n ,,,n, transforms 18 itself a Rn. n, transform, the 
operator AA, —A(A —1) and consequently any polynomial in (d/da) (z*(d /dx)) will transform 
an Rainy. oem, transform into another Es, a... n, transform. 


A. Example: 
We know that z"tu++K,(z) is Ry; v = ut 2n, n > 0. 
We can easily prove that 


d 


a r 
ər ("tata (2) | = ater AE, (a) 


+ (2n +2u+3)arters I3K (2) + (n i-r) n n F8) et E). 
Hence ; , 
[artats (2K, (2) + (2n + än + 8)omt 2, (z)] is Run. 
This can also be verified directly. 
My best thanks are due to Dr. B. C. Mitra for his help and guidance in the 
preparation of this paper. 
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FINITE DEFORMATION IN THE INTERIOR OF THE EARTH 


By 
P. P. CHATTARJI, Kharagpur 


(Communicated by Prof. B. R Seth— Received July, 10, 1952) 


Abstract 


Application of the finite strain theory has been made to find out the deformation 
n the interior of the earth considered as a self-gravitating, isotropic, heterogeneous 
sphere, density within which varies according to Roche’s law. Small strain theory was 
applied by Love to the same problem, density in which was, however, assumed to be 
uniform. Results of both the theories Lave been compared. 


Introduction 


The theory of finite strain has recently received the attention of many 
workers. It has been applied to various problems which cannot be treated by the 
classical theory of small deformation, (Seth, 1935, 1950; Murnaghan, 1937, 1951; Trues- 
dell, 1052). It gives good agreement with Bridgman’s (1988, 1985) experiments on the 
compressibility of media under hydrostatic pressure upto 10° atmospheres. Its success 
in solving various problems rejating to finite deformation has naturally given rise to the 
possibility of applying the theory to the discussion of problems relating to the earth. 


Applications io seismological problems have been made by Birch (1939) and by 
Biot (1940). Assuming that the initial stress is hydrostatic and pressure gradient uniform 
Biot has shown that the influence of the second order effects on seismological waves is 
negligible, though they may be significant in the case of the earth’s tidal movements. 
Jeffreys (1942) pointed out in discussing elastic instability that the relative importance 
of rotation and strain need main consideration; departures from lınear-strain theory 
ean be significant only when the form of the body is such that there exists possible 

displacements in which the rotations are large compared with the strains; within the 
earth such displacements are not possible at places of very great stress. 


- However. the application of the classical theory to the calculation of strains in the 
interior of the earth is, as pointed out by Chree (1891), beset with certain difficulties, 
even if we neglect the effects of rotation and other disturbing forces, as the internal 
stress is loo great to permit the direct application of the superposable small strain 
theory. One way to avoid this difficulty is to treat the material of the earth as 
homogeneous and incompressible which, however, does not correspond to physical facts. 

The problem of the sphere was developed by Lord Kelvin (1863) who sought to 
utilise it for the purpose of investigating the rigidity of the earth (Kelvin, 187@) and also 
by Darwin (1879, 1882) in connection with other problems of cosmic physics. Love 
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(1944) solved the problem of a homogeneous isotropic sphere held strained by the mutual 
gravitation of its own parts and showed that the radial strain is contraction within the 
surface r = r, {(8—n)/(8+8n)}*, but it is extension outside this surface, where 7, is the 
radius of the sphere and n Poisson's ratio. That this result cannot be applied to the case 
of the earth was pointed out by him, the strain inside the earth being too great to admit 
the application of the olassical theory. 


In the present paper it is proposed to apply the finite strain theory to find the 
deformation in the interior of the earth, assuming it to be a delf-gravitating isotropic 
sphere and taking into consideration its heterogeneity of structure, the density varying 
according to Roche’s law (Jeffreys, 1924). As done by Seth (1986) the stress-strain 
tensor relation is taken as linear. 


The body stress equations of eqilibrium 
i, jt Fi = 0 
refer to the condition in the stramed state. In any actual case the boundary conditions 


are also given in the strained state. Hence, as stressed by Seth (1985) the Eulerian 
components of strain should be used. 


As experimental evidence is lacking, results cannot be verified, but as we can ` 
anticipate, we get the dividing zone nearer the centre than what is given by Lhe classical 
theory, the difference being of the order of 10% and 13% respectively according as we 
= a particular or a series solution of fhe differential equation involving ß, zn 

= f(r) a function to be defined later on. 


Self-gravitating, Heterogeneous, Isotropic sphere. 


Assuming the earth to be spherical and at rest only under the mutual gravitation 
of its own parts, r = a being free from traction, where a is the radius in the strained 
state, the symmetrical radial displacement can be written as 


u, = r(1—$) (1) 
where - 
r= git +a? and 8 = f(r) only. 


The finite components of strain are given by 


“ 


e = 417 go) -r° (sa 288) e 
Cop = Gen = $/1 - 87) (2.2) 
8,9 = €g, = Der = 0 (2.3) 


where BI = dß/dr. 


The eumponenís of stress are given by 


en RSS (pr -. . [a-e m (^ EA (3.1) 


66 = pp = À [3-89 - in Kal .— ge 


it 
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rÜ = 69 = gr = 0. 
where ğ” = dß/dr, A, u being elastic constants. 


The stresses have to satisfy the body-stress equations. The equation : 


Ser), 101150. d 000). 100. ee UR 
- 36 Leer 21577 gp +76 cot 6) of, 








where p = o(1—kr?) [Roche's law, o = 10:10, ka? = 0753], F, = -9- 


gives 


S Ton + Gap? Q2) (9° +) rt ay / arar] = — 29 e Ae (4) 


The other two body-stress equations are identically satisfied. 

Henee 8 which is a function cf r alone can be determined from the differentia! 
equation — " i 
(BA + Bag? + (A + 94) 7? De + - i5 J pdr = 29T (5 ZC k, 

r a 
where k, 18 a constant of Integration. 


If we put 





the above equation becomes 
, 26 š Z go 2 krt 
88? ( 2 sən 2c f 27 -- 8 _( 1. +k 
P + BT — r” 2e J g7rdr “əv a 
where k, is also a constant. 
It may be remarked that since $ > n > 0, we have 1 > c > 0. 
Let us again put 


r = ta, p = t 
sə À + 2u 





and we get ka? 
28? + (219" -- B)? + 4c J ip"di = —P(t——— t*) +k, (5) 
where 8' = dg/dt. 

Determination of 8. (a) Particular solution. 

The differential equation involving 8 is a non-hnear one, It has not been possible 
to find an exact solution of the equation. Before proceeding to find a solution in the 
form of an infinite series, we note that we can get a particular solution in the form 

Box A,+A,t 
Substituting this value in (5) we get 


z 
9(A, + Art)? (A, + BA, T)? + 4c [tiat = pie £2) +I, 
Comparing the coeffiecients we get 84. = k; l (6.1) 
104,4, = —P (6.2) 
: i ku? 


5—184)P—3 
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Let us now turn to the boundary condition viz. rr = 0 ab r = 0, where a is the 
radius of the earth in the strained state. Since all the unknown constants in the solution 
are definitely determined, the satisfaction of the boudary condition will give us a relation 
between A and u. By (8 1) 


Gunu... 











or A,?(11—2c) + A.4,(10 — 46) + (8 —26)(4,-1) = 0 at T= a 
or goa 11—2c ka” geg 10- 4c + (8 —2c) | 2ge(11 +2c) -1} —0. (7) 
A+2ull+2c 2 A+2p 10 ka”A + 2u)10° 


Taking 7=0.5(u=0), c=0, g=VS1.262 cm./sec?. (value at Potsdam) a= 6871km., 
c = 10.10, ka? = 0.753 (mean density of the earth = 5.526), we get 


A = 0.53 x 1017 dynes/em.” (8.1) 
Taking y = 0.33, c = 0.507, we get 
A+ 2p = 0.49 x 1015 dynes/cm.” (8.2) 


According to seismological observations within the earth the value of A varies from 
0.74 x 101? dynes /em.? to 12.6 x 1012 dynes/cm.? from the surface upto a depth of 2898km., 
below which is assumed the existence of fluid layer, mainly from the fact that S-waves 
are not propagated through it; p varies from 0.09 x 1017 dynes/em.* to 8.08 x 101” 
dynes/em.? and n from 0.269 to 0.5. (Bullen, 1947). 


At the surface the thenretical value obtained here is nearly one-fourth of the value 
obtained by seismological observations. 


Case I. n = 05(u-0)c-0 


By (6) | 
e = 1'867— 0:698 7, (9.1) 


By (2.1) the radial strain is given by 
e, = 311— (6 7877) 


2 \2 
er I1-(4.+34,5) \ 
2 a 


2 
-ih- (Da: 1:914 7 E 
2 a 


radial strain is contraction within the surface r 0'678a and extension outside this 
surface. 
Case II. , = 0'88; c = 0'507 


ge . 
B = 2029 —0'628 —," (9:2) 


The dividing zone is given by r = 0'785a, 
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In Case I the classical theory gives the dividing zone as 7 = 0"74$a and in Case II 
as + = 0'818a. 


The particular solution of the problem by the application of the finite strain theory, 
therefore, gives the dividing zone nearer the centre, the order of difference being 10%. 


(b) Solution in Series, 
Assume 8 = A, “ Alt - AQU-E AD... tH Ant” + 
Substituting this value in 


282 + (218! +8)? + 4c J tps dt = — P(- t) th, (5) 


wheie 8” = d8/dt, we get 
2(ÀA, tA tA HACHA ltt ...)? 
+(4,+34,t+54,0°+ 74, 19A UM .. .)2 


+40 [14, 424, 9 84, +440 + ... dt = — P((— ən k, 


Comparing the coeffecients we get 


34. * = ky 
104,4, = —P 
Pka? 
2 _ 
14A,4,+(11+2c)A,’? = — (10) 


I 
184,4, (84 +165 A.A, = | 
0 


(4c +27)A,?+22A,A,+ (46 + 6c)4,A, = 


and so on. 
We express all the coefiecients in terme of A,. By (10) 














A, P 1 | 

A, na 

"r 7 qa (+20) xs ga) | 

4 = BIH (24+ 16 aN = Geh —. m às] (11) 
T Z x İ(46-- 6o) — (666165 Är is — (11 +20) ast 








1 E a 1 papa (11526) 1 , (11-7 20)" a 


+ (4c +27) — 14: 1 7, — 107 ^A, Toe A, 
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All the constants except one are then determined. Ihe unknown constant 18 Lo be 
obtained from the boundary condition. For that we have to take some tentative values 
of the elastic cunstants. The variation of these constants within the earth have already 
been mentioned. However, we make an agsumptien about the mean value of 
A +2u(= 9 x 1012 dynes/cm.?) from the known values of A and x obtained by seismological 
observations. TL 


Substituting then numerieal values in (11) we get (corresponding to Case I, 
q = D'Biu = 0), c = 0) retaining figures correct to four decimal places, as the mean 
density of the earth is known correctly upto three significant figures, 


A, 1 
— = — 0702 —. 
A, un 
| 
A, 1 mi 
—3 — “0189 — —" E 
A, ER k | 
(lla) 
A, xL öl ` 1 
—3 = “0025 —— — “0005. — 
A, ui d 
A 1 “0006 “0001 
“4 = — 0004 — + — - — 
0 A A A, J- 5 
and so on. 
Therefore 





1 1: | 0001 
A pes 0510001855 ee, 
d CR ARS Aş 


when t= 1. The numericul value ın the last term is doubtful as other terms containing 
1/A,° have not been taken into account. 
By (8:1) the boundary conditian can be written in the form 
8—2c = (2tB' + B + 2(1— c)8? l:-i. 
Substituting numerical values and retaining terms upto 1/49“ 
ARCU 4 A 
l 


By the method of iteration, taking positive value of square root 


A, = 1'070, A, =—0'066, A, = 0'015, A, = 0'002 








1 _ 01406 00008, 


3 4 
8 = (1:070 —0:006 + 00157 + — (12.1) 
neglecting higher order terms. 


The dividing zone is given by r = 0°648a. The negative value of square root leads 
to the same result. | 
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Proceeding in the same way we get corresponding to Case II 
A, 1:080, A, = —0:066, A, = 0'014, A, = 0'002 


2 4 
Bz (17080 — 07066 = 0014 + (12.2) 


The dividing zone ıs given by r = 0'707a 


As io «1 and the coeffisients Ao, .A,, . . . rapidly diminish and also from phy-ıcal 
considerations we may expecl the 8 series to be convergent. 


The series solution gives the dividing zone 4% nearer the centre compared to the 
particular solution and 18% nearer the centre compared to the solution obtained from the 
classical theory. 


The depth from the surface of the earth below which the fluid layer is assumed 
to be situated is 2898km. and the radius of the earth consdired as a sphere is 6371km. 
Therefore, the fluid layer extends from fhe centre upto a distance r = 0'545a. The 
dividing zone in any case lies beyond this distance from the centre. 


In conclusion, I wish to thank Pıofessor B. R. Seth for his helpful guidance during 
the course of the work. 
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KHARAGPUR 
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ON AN EXTENSION OF BLASIUS THEOREM 


By 
Kamint Kumar De, Calcutta 


(Communicated by Dr. S, Ghosh—Reoeived July 22, 1952) 


It is known that in the case of two-dimensional steady irrotational motion of an 
infinite liquid round a solid cylinder, the resultant thrust of the liquid on the cylinder 
can be calculated easily with the help of Blasius Theorem (Milne-Thomson, p. 168). In 
this note it is proposed to show that in the case of rotational motion also, it is possible to 
calculate the resultant thrust by an extension of Blasius Theorem. 


-Let X, Y be the components of the resultant thrust along the axes and M the 
moment about the origin. Ifl, m be the direction cosines of the outward drawn normal 
to the boundary of the cylinder and p be the pressure then 


X = - [plds — — (pay, Y = - İpmds = İpda (1) 


X-iY = -i [pda (2) 


and 


M = f p(ma—ly)ds = f p(zdz + ydy) (8) 


the integrals being all taken round the contour of the cylinder. 
In steady rotational motion, the pressuse at any point is given by 


p = c—1oq* 


where c is a constant along a stream line. Since the boundary of the cross section is a 
stream-line, c is constant on the boundary. Substituting in (2) and (8), we get easily 


Zn f q”ds (4) 


M=-% f q*(zde + ydy) = —4eR / q*adg (b) 
C " e 


where R denotes the real part. 
If y be the stream function of the liquid motion, 


u=- pa 
Ou Or 


Therefore 


(u+iv)ds = ES de: SL dy e dy 
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Hideo di 
(u —1v)da = 2777 idy 


On the boundary C of the cylinder dy = 0, so that on C 


(u+iv)ds = (u—1v)de 
and therefore 8 2 
g*ds = (ut+iv)(u—ivjdz = (u—iv)*da. 
Hence 


X —(Y = tig | (u —1v)! de (6) 
| 
x —toR | (u — iv)!zds. (7) 
| 


The formulae (() and (7) are extensions to the case of rotational motion of the correspond- 
ing formulae given by Blasius for ırrotational motion ‘These formulae have not the same 
advantage as those of Blasius, for the integrands in (6) and (7) are not analylıc functions 
of z. But this difficulty can be overcome by replacing the integrands by analytic 
functions which reduce (o the integiands on the boundary C. 
Let the boundary of the cylinder be given by the curve £ = £, in the transformation 
ec, C= £t in 


When the liquid motion outside the cylinder is known, u —iv is known as a function 
of z and z, and therefore as a function of C and Ç Now, on the boundary of the cylinder 


ç = 2£6,—t, so that on the boundary u —iv can be replaced by an analytic function of ¢ 
and Lherefore of s. 

As a simple application let us consider the motion of a liquid of uniform vorticity 
round an elliptic cylinder. in this case, we use the transformation 

2 = ceosh¢ 
and let £ = É, be the boundary of the ellipse 11 the velocity at infinity i8 due to shear 
motion parallel to z-axis, given by 
u = — U —2wy, v20 

and if we have in addition, a circulation « round the cylinder, then 


1 


u—iv = — U +iog —iwe cosh &— |Ae-t- 2iBe- + x | 7 
c sinh ¢ 


2n 
where b 
A = —Uce sinh £, = — U z (a+b) 
B = —iec'e*6(1— cosh 2£,) 
a = c cosh és, b = camh ğ. 


This complex velocity satisfies the equation of continuity, gives 2w as the vorticity 
at any point cf the liquid, satisfies the condition at infinity and makes the boundary of the 
ellipse a stream line. 
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On the boundary £ = és, Š = 26. —Ç, so that on the boundary 


u -iv =—U + twa - təc ecsh (26, — 9 - |4e- (— S Boa- 2 + x] 1. 
i Cot he d 
Now, 1 1 
ez = cosh C—sinh & = — (2 - (27 — c?)}), e72 = "a IO 
C 


Henee, on the boundary 


u—iv = -U-ioz — iw cosh 2£,.2 + io sinh 2¢,(2? —c2)t 


5o 71) “x sem ə tl t- ei ir zs |. 


Substituting in (6) and (7), we observe that the integrals can be taken over a large 
circle 2, with centre as origin, as there is no sıngularitiy of the integrands between 
C and 3. We have 


5 Br Duo" á 1 
tat, Weimar 
Expandıng in powers of 1/z, 


u—iv = — U 102 —iws cosh 2€, +to sinh 26, (2 = 


LI | — 
al 





1 s i K Ac 
= (ol - cosh 2£, + sinh 368 U, (we sinh 2£, + Jj gor a di 
5-75 U- (u2ab + * = ae (8) 
Gr 22 mu: 09 
(u-v ? = ... —QAciw 5: +12 (aab + el x real terms. 
CR 2 T z 
Now Ie di 
X -1Y = hip x Qni[iU (2 wab + - an m 2Ab(a-b)] 
= —trp[U (ch + E ) +2000] = —irp E + U2ub(a+ dd 
"ra TT 
Hence 


X 20, Y = pUk+rpoU2b'a +b). 


We now observe that citculation immediately round the cylinder is given by 


Jude vay) JI Se de + St dz) = In ée 


From (8) it is seen that the integrand has a pole at the origin with residue 
8—1840P- 3 
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- iens 2) 


circulation immediately round the cylinder 


— ani} - 5 (2unb + &) = 2radib +x = x’ (say) 


T 


Then 
Y = oUk!+2r0wUD* 


and lastly 
M= -HR [tdi = — ko Rani x roal terms = 0. 
0 


Thus we see that the force exerted is a uft force and the m ment about the centre 


of the eylınder ıs zero. 
In conclusion, I wish to express my thanks to Dr. 8. Ghosh D Se, for his 


continuous interest and kind advice regarding the presenta.ıon of this communication. 
Reference 


Mılne-1hemson, 1940) Theoretical Hydrodynamics, 9nd ed. 
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ON THE GENERALISED HYPERGEOMETRIC FUNCTION 
OF APPELL’S TYPE 


By 
SHANTI DARAN, Lucknow, India 


(Communicated by Dr. S. C Mttra—Recetved May 6, 1952) 


Tbe object of the present note is to obtain certain generalisations of Hypergeometric 
Functions of Appells type (Appell ot Kampé de Fériet, 1926) by increasing the number 
of parameters (Burchnall and Chaundy, 1940; Bailey, 1929, 1934), the number of 
variables being three. We have obtained relations betwoen them and Appell’s and 
ordinary Hypergeometric Functions and have also obtained partial differential equations 


satisfied by them. 


1. Let us write 
Pi (zı, Zar 93; Bi; Bs, Bas y» Ya: Ya» T, Y, g) 


SE Sr BB ary ti) 
en (yrly)ayderisiti 


t=0 s=0 r=0 
P(t, Ae, Œs; B; Ba Ds; Yi» T, H y, g) 
N >>>! gll 0. AR (a..) Pus (Ba) (Bake xTU t (2^) 
DET T &ltl 


län, 94, Fy; Bus 6. ; yo Ya, Z, Y, z) 
TÑ (a a), (& (a a5): (81) 182), (By): t gy sgt (87) 
(yi r(yalare Tl 81 6l 


tel s=0 rel) 
Pali 75, Bal Bi; yi; 2, Y, 2) 
co (er arytat (4^) 
(Yılrzarı 7 tel t! 





i20 s=0 r=0 


g;(,, 72) Üş) B; Yu Ya») €, V, &) 


UNUS (xı) (a). (5007 a qa i e 
= ll tE eme (57 
222 (Yıl'Yalart risit! 


Polti» Lo, Ha: Bi 5 Yr Yə Ys? T, Y» g) 


= > o cc DM dE en z^y'*s! (67) 
120 s=0 r0 (Yi), ck (ysl? lal £t 


g,(*,, Gan Da 5 Bi, Bos Yo Ya © H z) 


a 


- 222: (o), (24 ) do, (B, js (8,) rət z* sgt (7°) 


<—  (vidrlyalare rl at tl 
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BER V2» P By; Yı Ya> T, Y, 2) 


5 ə 0. İİ İİ aa | 
. (yu (ya) UL -- = 


{=0 s=( r=0 


\ AR (B, İn £4 gr Q’ 

-2226 Cice etary e) 

The series (1), (2), (8), (4), (7) and (9) are absolutely convergent if la] <1,/y|<1 and 
Jal «1l. series (5) is absolutely convergent if İzl-İyl <1 aud İsi « [y] or 


lz|+]|z| < 1and|y]< 1. The series (3) and (8) are absolutely convergent i£ İs |< 1 
and |y | | 2| < 1. | 


Polt; Bis yú X, Y, 2) 


We use the symbolic ope: ators (Burehnall and Chaundy, 1940) 


ThI + 8’ +h) Ah) = E@+h)T(8' +h) 


vü) = s PG hy JI DD +8 +h) 





We shall simply state the results without proof. ‘The results oan be verified easily. 


Pal, a, tas Bi, Ba, Bi, yi; €. y, 2) 


- Ş S Cay es Bde Bala 
r= s==Ü (Yi) ras T l 8! 


x FtUla „+8; 8,+8, B, +83 yy Tr+28; y,r)* (1) 
95(2,, 1, 4): B, Bas Bi: Yo Yağ 2, V, 2) 


-5 S (—)*(z ) (a, (8, )i (8), Bla yars 
0.) (yi) lya) z T! 8: ; 
x POla, +s; B,+8, B,--8; yi F8, ya FT+8; z, z] (2) 
Palto di, 445 B, Bi, Bi Yr Ya; T) Y, £) 
=> b (z). (za) (8.) Ba) — 1 — 8. 1) 
en (YılrlYalrra rie! 


x Flaş tr; By Hr; yu r, Ya 8r; zum, z — ez] (3) 


xz y's” 


guia, Ba, Ss; Bi) Yu © Ys s) | 
Il (a Xa), EN (81) eit ara 
ns (Yılrıas Tİ 8 | 
x FO [2,--8, 2,48; Bares, B,47+8; yy+r+2s; y, 2] (4) 
Pa tis Zi, Fe; By; yi) Ya; Z, Y, 8 a) 
-$ suz ) (2) AM iar” 
xq (Ya age 3 | 


—.. 2r; atr; yı T, ya T--8? m, 2] (5) 


+ F's are Appel”s Fun ctions. 
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=Y Y ken EE 


(yi); 7 8l 
x FO[B,- 7*8; a,41; YT, Yatr+s; y—yz, 2— ye] (6) 
p. (di, Xat 2a; A; yi? £, y, g) 


KE EEN ARCA DER T! 8! 


TÜ swa) 

x FÜ[B, 728; 2,8, 3,8, y, T r+ 28, y;t!rt2s, y, 2] (7) 
Pal, a, C35 Bu Bas; Yi) Ya: Ya» €, Y, 2) 
(Yi): (ys):(ys)+ ris 
x FP [a +r; Bir, B, r8; yit, y tr, m 2] (8) 
SS all Tr ese 
(Yılslyalelys)r rial 
XPOTB AT; aşmır, y, T, ya +T; y—ya, z—ya] (6) 
Pellis dz, €; Bes Bad yas Yü yas €, Yo 8) 

-5 > & (a, IRCH ) (a VHB; )r (Ba) CA? auser 
(yal) Yılar ris! 


1 Ü) am) 


c! ya" 


r=0 s mi 


T mü Arca 


x POIs +r, Aa +T; B +7, B.,+r+a, y +2r; a 2] (10) 
L$ $ e elus ayni 
Yılrl Yı )e (y) T! 81 
xFO[8,c2r; sr, aşır, y, t2r; y, z] (11) 


T —Ü s 0 


Pelti, Ai, As; "nu öz Yo Yi: Yas Z, Y, 2) 
- > > abl x.) (B, e (Bs ) al? 


—yı) r zu gT 
(Yı Mey, arlys)s ria! d 
xPOla +r; Bitr, Ba+r+s; y,*2r; z z] (12) 
GC Ga, 133 Bi, B25 Yu Yn Vas 2, V, g) 
= BS $c LY GPs Bs ra 
MERIT (Vids 


Tw) 3 —Ü 


T =Ü s» 


x FOB, - 2r; Za FT; yat iüuya TT; V, sl (13) 
Pelti» Ta, Ha: e 5 Yr Ya» c, y, g) 


(a, IRCH Jl, (z ) (2, )ə(z,) (8; (6,), TaT 
 (y(yar ris. rlial —. 


dé 3-0 
xF®[8 +r, Patri ayri G tT; Yet 27) y, sl (14) 
.. (I (6.25 1:6 zöu"z" 


(y1)s(Yalar [ya +r—1), rla! 
x FO[B,+2r; a+r. as +T; y, 2r, yat2r; y, 2] (15) 


Tal sm 
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geän Ti, Hai By, Bas Yo Yas Z, Y, 2) | 
(ie dB Boles 8 at y'a 
re s=0 (Yırlyalr,arl 8 
x FOfa, +r; Bitr. B,+r+8; yıt nyatrts; T, z | 
l 


Palto 955 Br, B3; ue Ya; Z, Y, 2) | 
>» E (42). 8 )4(83), 82), (Bah gyar 


nva rÍ 8! rls! 
x Fo fa, +2r; Batr, Betr; ya +2T; y, 2] 


-S3 7 (= do: Jer (By)s )a(B art ye” 
= 2 (Yalar T! 8! 


T xÜ s=0 
x pt) e, ta Bat 2r; ya 2r, ya 2v; y, al 


rel sof 


s P, in r ... 


x FO la, + Or; Batt, Batts y; 2r, ya 273 y, z] 


rTA0 s 


3. Relation between $-function & simple Hyper-geometrio function : 


94(2,, Qo, Ha: B, ) Yi» T, Y, z) 
-55 (—) ur HE (7 E E 


( risi 
1030 yı) r+28 


x F(B, Ar äs, 0,4 8; y,-Fr--28; y+e) 


(x), (a.)4 8... iy RN öz), 


x’ Vie? 
(Yılrıza rls : 


TÜ ges 
RI, kima, 0,8; y, 128, y +z- yz) 
Gala, Xe, Qa, Bi3 yu Yas 2, Yr 2) | 


-S S edet s 


en 0 (yidely hast) 


x F(5,-tr*8, x, T, vit ni g)F(B, vro, Os +T; ya ETH 


EE EN toy Za: Bi Pa Yis Yo CR T, Y, £ s) 


wm 


- N («,),(24),. "oh (8)).483) (Ba)y ot 
b THE (y), Tİ ei r! 8i S Y = 


x F(B,- r, Zune Yat"; y)F(p, +7, Gott, ya n; 2) 


KÜ s —0 


Pelti Zar 23; B, Pa; Yu Yas Ya; Z, V, s) 
= > f) (-)r (a ARC „(B)s(B;)arl 297 as), zarar 
eA (a !eYa) ar rla! 
xF(ß,+2r, a+r; ya+2r; y 2) 


T we() pret 


(18) 


(20) 


(22) 


(28) 


(24) 
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Ə, (dı, Za Ha, Bis Bz; Yi) Ya? Z, Y, 2) 

EK 

(Ya T r—-1),(YısrlYıla rl s! 
x F(B,- r, XT; Ya tar; y)F(B44 1, A, HT; Ya 2r; s) 

ƏT Hi, a> Bu Öz, Tr: Ya, T, Y, 8) 

-S S Create ula, s, 

ly alYalor 71 81 
x İB, 2r, a+r; y, 2r; y+z) 


T ss gx 


T m göəmü 


- VN fa ana (B1) B5) Cy — Bə), göyə" 
(Yılslyalar rl 8! 


x F(Bt T, a, +T; yet 2r; yz — yz) 


reali s=0 


Pali a3; Bi, Öz Yo yai 2, Y, 2) 
SS (ala zuer 
aya) 71 8! 
x Rio, kr, Bott; yat T; y4 2—y2) 


rm) s =0 


Polar; By; Yis Z, Y, 8) 


ki 


= (Bde bt, 47: B TI yt ^; y 8) 
d (yi), T! 


fo 


"AZ Een 
x Pla, +tr+s; B, r8; y T8; y+2—ya) 
3. . Partial differential Equations: 
Let A, B, C, L, M, N, P, S, R stand for 
Sw Fw Fw Fw Fw Fw Ow Ov Qu 
da!" Oy?’ Gs?’ Gedy’ Tyde Garde Əz” Gy’ Oe 
— £(1—2z)À t iy, — (xi + B, +l)a}P—2,8,w = 0 
ei yYl—y)B+tys—(ar+ Pst Dyl0 — a 8, = 0 
_ £(1—2)C + fy,— (a +B, 1)g]P—a,B4w = 0 ` 
— e(1—z)A* y L czN (y, — a, 8,“1)zi P-a, pw = 0 


T est) 5 mm0) 


şil y(l—y)B+eM+aL+fy,—(a,+8,+1)y}Q—a,8,w = 0 
| a(l—s)C+aN+yM +}y,—-(a,+8,+1)a}R—a,8,w = 0 
— e(1—-2)A * yL + fy, —(o,+8,+La}P—a,8,w =0 ` 

Pa| y(1—9)B-* 84M y, — (ss 847 1)y] Q — aw = 0 


sl —2)C +aN t fy: — (a3 84 + 1)2] E — aBa w = 0 
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(25) 


(26) 


(80) 


(81) 


(32) 


(58) 
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z(1—2)A +y -gL t {y — (a, +B, +1} P — sayi — aR + 2(l—2)N -u,ß,w = 0 | 


e| y(1—39)B +a -y)L +a- y)+ fy,- (aa + Bi + 1Y} -a ab sg = 0 | (64) 


= z(1— a)0 - &(1—2)N 4 y(1—2)M + fy, — (aa +8, + 1)a}P — a,€ P — «58 ,w = 0 
T a(1—2)A — ayL — zaN + fy, — (ai tB kel P—o,yQ au af =0 © 


el y(1—3)B-yeM —zyL + fy, — (aa + 81 1)y10 -aP — a8 R aß = 0 x (85) 
_ a(1—2)0—zzN -yeM +Íy,— (aş +B, +1)z2]R—a P —a yQ — By = O _ 


— a(l—a)A tiya +B +IP- agt =D f 

pa) y(l—y)B—yeM + 44, — (a, t B,-- Dy]Q -aR — apat = 0 (36) 

at — alt -yM+ly- (aş + By + ahh — aş) — abw = 0 _ 

z(1—z)Á + ba + By + DEP — çın = 0 

o,| 1(1- y)B x 66 —y2M + fy,- (a: + Bç : DylQ — a3 R — a8, w = 0 (87) 
— s(1—s)0 4 yB —yeM tee (a, 8+1} R- Qa, = 0 _ 

— z(l—z)4 Ey, (aa +8, + 1)zİP — ğın = 0 | 


Pa) yl—y)B+2(1—2)C—2y2M dy, laa + Ba Dy] Q — (a, +B, + D) R — opto =(0 | (88) 


— 2(1-2)O+y(1-y)B-2yeM + fy — (as B, + y)e}R — (o4 8, +1)yQ -a,B,W zü. 


| z(1-2)A +y(1—y)B 4 2(1—2)C — 2zyL —2y2M —2z2N 
| (a, B, + DíaP + yQ + Rz) - y,P— a d.t — 0 


| z(1-2)A+y(1-y)B+2(1-2)C-22yL—2ysM —2z:N 
Po | 
— (a, +8,+1)(@P + YQ + Hz) yığ — ağı = 0 


z(1—2)À -- yf1—y)B Hr 3(1— 2)C -2zyL —2y2M — 2ga N 





— (a, 8, + D) (£P * yQ + Ra) — 4, R —a,B 0 = Ü 
4. We know that 


Galti; Ha Dax 81, B5; Yı Ya; T, y, g) 


| 
' DEE) 
Dy, a 0.1) (Ya m z, )1 (Ys KR nal (aT ies (o) 


1 pl 1 

x J f J un-l(] —u)-a-lyes-l(] —v)n-a-hpges-l(1 — 10)f:—a—1 
0 0 0 l 

Xx (1—uz)-&(1— yv — sw) Adudvdw. 


U 
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Let us make the following substitutions š 
(1) u=1-t, = 1-—s, u = p 
(2) üz 1—t, 2 = s, tp = p 
(3 u= t, v=1-8s, w=1-p 
(4) u= t, v = 1-8, w = p 
(5) u= t, 0 —8, w=1-p 
so that Pelar, aş, ga" Bi, Bas Yo Yas 2, Y, 2) 
—r — 1 2 
= (1—2)-5 (1 =Y) ^9, (n — dı, Ya dy, Gaz Bi, 8:5 Yu Ta; 1-2 o y) 
ES (1-2) Bolyı-a,, Ze, Us, Bı, B4; Yu Ya? — y, z) 
— 2 
= (1— y~z) So, = Ys Ga: 8, Bs; Yu Yi» 8 mem =) 
g 
= Droe, (ay, Ya 72) Da: Bi, 85; Yu Ya, €, m +.) 
— £g 
SS (1-2) 78:9, (s. Ga Ya — 035 Bi, Bi) Yu Yas Z, per. =) (40) 


B. Relation between 6, and 3, function. 


Play, Ga 95; Bo Yu; Z, Y, z) 


I'(y "i : x 
uñ, I (1— u)r- B, 1 1-— — m) 4 l—u — 04 iS uz adu 
noms ah ( Iris ) 


— 





Put u = 1/t, we have 


mco hr V. a, tostery (4 — | v1 Bi -1t — seg (taat te 
Pen 78)J, ' n (E—1)n7&-1(E— z)7« (ty) (t— 2) "dt 


(Sen) 


and 
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substituting we get | 


GEN he; ys Big Yi» T, H, g) 


EAW z y : z z 
= üə e 95 Bi Hax X35 da Toa taz; a, Has tas, Yi; T, te 


SEN (2) 9 İs, Ga; B1; EEE Yi; jak. y, bes Y) 
z z a S 


4 


g a £ ay ! a 8 
ə E E e (ns 43, By, a, Hag aş, a, baa tas, yı: 1— sə ps y? d (41) 


In conclusion, I wish to express my indebtedness to Dr. B. C. Mitra who kindly 
suggested the problem to me and took keen interest in the preparation of this paper. 
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VERSUCH EINER LOGISCHEN QUANTENDYNAMIK DES 
ELEKTRONS. (II) 


VON 
S. C, Kan, Calcutta 


(Communicated by the Secretary—Received August 18, .1958) 


Inhaltsübersicht. Es erfolgen unten in Fortsetzung des ersten Teils unserer unter 


dem gleichen Titel erschienenen Arbeit die Anwendungen der dort vorgebrachten Theorie 
auf einige wichtige Beispiele. 


1. Worbemerkungen 


Im Anschluss am ersten Teil unserer unter dem gleichen "Titel erschienenen Arbeit 
wollen wir im Nachstehenden die Anwendungen der dort vorgelegten Theorie auf (Kar, 
1952) einige wichtige Beispiele behandeln. Als dem wichtigsten Beispiel wenden wir uns 
zunächst dem Fall des Elektrons im atomaren Felde zu. Danach wollen wir nacheinander 
noch folgende Fälle ın Betracht ziehen, nämlich, den Fall des’ Elektrons im konstanten 
elektrischen Felde, den Fall des Elektrons im konstanten magnetischen Felde und den 
Fall des feldfreien Elektrons. Jedes Beispiel wird auf zwei Weise bebandelt, nämlich, 
nach den quantendynamischen Gleichungen [B,, A,| bzw. [B,(a)j], A,(a)] des ersten 
Teils und möglichst unter den Rubriken: (1) dıe unverkürzte Quantenbahn, (5) die 


klassische Bahn u. (iti) die Quantenbahn bei nullter Näherung bzw. bis auf die erste 
N&herung. 


Die Numerierung.der im vorliegenden Teile (II) auftretenden Gleichungen bzw. 
Sätze ist als unabhängig von der im Teile (I) geübten zu verstehen. Nötigenfalls 
geschieht ein Verweis auf irgend welchen Satz des Teils (I) explizite. 

2. Das Elektron im atomaren Felde 


Aus (I) holen wir die Gleichungen (26) bzw. (26(a)) her. Mit der Setzung 


=“ y, A 0m da — 4, 
T 
































lauten sie 
y d, 
IA +6 y) t öheV)) HƏ +8, : ; 
° d, de (i 
9 + + y^ d, 
IA «ev )- ene va] Ps ; 
SS de d, 
bzw. 
id x 
erlech ev) = — 86 : 
ii "s S (2) 
== = = — £5 , 
AS - ev )- ense vo] i b 














2—1840P— 4 É 
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(a) Einige transformationstheoretische Sätse. Wir wollen nun emige bekannte 


zur Transformation auf Polarkoordinaten nützliche: Sätze vorausschicken. 


Mit den Bezeicnnungen Ä 
z-—rgncoso ] 


y = r sin 0 sin d } (3) 
g = r cos 0 | ) 
zog... 2 O N: 
Lig x (v2 e). | 
E 2), 
m (= "8s/! I M) 
. | 
o QA d 
us 9. 9) s. 
2 Toy Yar "ae 
- + i wa 
' M = L+ =<, (5) 
: 
beweist man leicht den Batz Ä 
8 ++ ə ++ 
M, | (host ev) + chen) = (+ af he V) |M, 
| (6) 


Se Je um pa me 


9 m dim o 5 ++ 
M, Ae + eV) —ch,{o vi| = Iv + ev) —ch,(o vi M. 


+ + | 
Mit der Bezeichnung: (c V)- i = K, beweist man ferner den Satz 


[Ue +eV)+ oh,teV) |= IO + ev) - el,(ev) İF, | 


| | ; (7\ 


i (red, em)- e] -[( em) ens eo]; 


Diese zunächst für die Gleichungen (1) vorgesehene Sätze gelten offenbar auch fur 


das Gleichungspaar (2), ' 


Ebenso beweist man allgemein die Sätze 


M,K=K.M, ` | (8) 
hiz Ek GG | | (9) 
GI T | (10) 
T | . 
+) +7 
Lei K+ K. rt )=0, (11) 


( 
(rV) = (22) er Mi | (19) 
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(b) Schrittweise Transformation des Gleichungspaares (1) bzw. (2) auf Polarkoor- 
dinaten, 


Auf Grund des Satzes (6) setzen wir 








véi fel 1 
(Der go) : i «| I (13) 

l d d, | 

bzw (Un 5 2 ji = o! e à j 














wo 4 eine Konstante ist Wir erwägen nun die zwei Moglichkeiten : 
(a) km-i) 
(b) p=—(m+), 


wo m eine ganze Zahl sei. 


Entsprechend der Wahl (a) findet man leicht 























d, Pi) morto, (s, pa) otn no, (14) 
enisprechend aber der Wahl (b) 
(Pa, Py) m e Orte, (uu, u^) emo, (15) 
Auf Grund des Satzes (7) setzen wir 
ul ıı” m. Yı ) 
: d a | (16) 
K [Vy 2 fred) anl” |=ik ff | 
a! del d) 














wo k eine Koustante ist Wir fuhren nun die Bezeichnungen 























by ty) = f, 10 -4)-9, | (17) 
EU, YA) = fo, MS Al = gç 
ein und erhalten 
++ h 1 ] 
[ (eL) ti] == ; | 
{la fa [ (18) 
++ 19 ü | 
[(oL) +i) =—ıh i 
0: pA ) 
Hieraus ergeben sich die Gleichungspaare 
„I 9 57: 
e [S + cot “ə m fk + (u-—1)f, | Sé 


art Il 


er E i cot İL. = {k-et} f : 


136 8. C. KAR 


Í b 
und ! 


"Ii cot ed = -k-(- pig, | - 


zu s O l. = 
—e is i cot dÉ, ik Fr (n 3) gi | 


Entsprechend der Wahl (a) obeu findet man die bə zweiter Ordnung wie 



































1 9 9 
0 k (k — 5 | _ | 
ın ae ƏM x heb 
(21) 
1 9 (m + 1)? | 
92 k(k— Z | 
[= T ET n ..- isin? 6 —. ) . 
und | 
ENT Si mö İ 8 ) 
Ls EE Deg +k (k 7 1) I xin? 0 gi = 0, | 
} (22) 
I ð d (ms) | 
pur k " " 
Ee zzgl ET +h (+1) naş |? 0, | 
entsprechend aber der Wahl (b), die Gleichungspaare 
l or, 2) (m +1)? | ) 
6 (k — 7. = 
= ¿aÍ sin 90 Tü M ı sin? 6 nr 
(23) 
1.0 >) m? ] 
m ie nen Weg k a E. = | 
l: əssin 636 FEI 1- “gin? 6 fs = 0, j 
und E 
| 
Eu 2) m +1) | "PE 
l: TE 635 + kk +1) = ig 17: 0 
| (24) 
1 393 3). 2 
Er re S ET Mk +) 


may 9-9 j 


Man hat also nun auch noch zwei Möglichkeiten in Erwägung zu ziehen, nämlıch 
entsprechend (z) k =+(l+1), (8) k= —(£+1),; wo l eine positive ganze Zahl sei. Im 
ganzen hat man demnach dıe vier Fälle zu berücksichtigen, nämlich, (a, «) 


, (b, a), (a, B), 
(b, B). | ! 
Fall: (a, a). | 
Auf Grund von (21) bzw. (22) kann man einstweilen setzen 
(fis gi) [A(r, t) Pr", Cr, t) Pam] es, | = 
(fa, gi) = [BC t) Pr, Der, t) PRED] etm, 


wo Pr", P-(m+D, usw, die bekannten Legendrefunktionen und A(r, f), B(r, t), usw. 
noch zu bestimmen sind. 
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In Ansehung aber von (19) bzw. (20) mit Rücksicht auf Rekursionsgleichungen der Art 


LG Š + (m+1) cot ot ppt). pr» (26) 
findet man leicht 
Bir, t) = Alr, t) . L— m), 
| (27) 
Dir, t) = —C(r,t) . (I+ m +2). 


Fall: (b, a). 
Setzt man auf Grund von (23) bzw. (24) zunächst 


(f, 9) (An, t) Pp l) Cr, t) pa tD] e-em tne, | 

(fa, 9a) = [B(r, t) Pr", Dis, t) Pri] eins, | (28) 
so ergibt sich nach (19) bzw. (20) mit Rücksicht auf (26) 
| Á(r, t) — — B(r, t) . (I-m) 

C(r, t) = D(r, t).(I+m+2) | SS 


Fall: (a, 8). 
Ähnliches Verfahren ergibt für diesen Fall 
(f... gx) = [Atr, t) Pul, Cr, t) Pr] em | - 
(fa, Ja) = [—A(r, t). (I+ m2) Prat, Cr, t). (Lom) Pp "+ lernen» 
und für den 
Fall: (b, 8) 
(f, 91) = [B(r, t) . (Lom 9) PREY, —Dfr, t). (15m) Py 091g» | m 
(fa, 9a) = [B(r, t) Pixi, D(r, t) Pr" ]e mo 
Mit Hilfe von (12) schreiben wir nun die Gleichungen (1), wie folgt, hin 














DE DOW 2 |- 
a DE DEDI I | 








Daraus folgen mit Hilfe von (16) 


> 
tar) (ə? Ne + 2h] 
ee h Ae We erm 
| (hom +e? one 27” 














1 ue Ob 0m 
(erem) DS SD M C D 111 
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daraufhin aber mittles (17) 
Tİ: 
( + ch H+ 
Tİ: 
79; 
(hs — eV- 2 
TGs 


Blickt man auf den Fall 


(fi, fa) = Á(r, HIP 
(935 9s) = Clr, gemeng, _ 


Le +eV— 2 








|+ chef © 





so hat man 


CC rg, 
O SE 
r 


cos 0, sin 667" 


reg 
w 














rg; sin Geng, — cos al 


mittels der Rekursionsgleichungen 


zez MI 
cos 0 PU 


p 0 Pryit (b+ m + 2) cos 8P (m+) _ —(I—mp)P/ tU. 


Fr 507 
C(r, t) 


— (Lom 4 9)P P *D gone De 


—(t+m +2) sin 6 Pr et = prt, 


AR 


igi 


"CH. 
“o 








83 
Däi i 


— — —— nn 








` (a, a) zuruck und schreibt man statt Pt 


(I+ m +2) P TP ertm+n. | 





Pr motto 


rC(r, t) (34) 





(L—m) Pp "tD emt)? | 


| 


(85) 


Danach lauten die radialen Gleichungen fur diesen Fall 


US +eV— o) (rA)— eh. ı 


dr 


ho 


Déi 


d 


+eV+ to) (rC) — hi 2.” 


Analoges Verfahren liefert für den Fall: 


——- 
0 


In gleicher Weise findet man fur die Fälle 
gen als 


( 


zı. d 


Q 
how tov + to 


9 
udu 


Ə 
soy 
lox +e +e 


.) (r0) - ch (+ 


) (rD) + TE 


: (a, 8) bzw. (b, B) die radıalen Gleichun- 


) (rA) — ai: 


(36) 
111 (A) = o 


——— E E ina, e ER enn Srt 


(b, a 


EE 


T 


)eni = Ü 
(37) 
EA 


rl \(rB) = 0. 


T 


— — À mmm -— 


ft) o) = 0, 
T 


Ig 


Ww "7 omnem 1“ 
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bzw. 


9 ) (4.- 1-1 
hy-— eV —= = = 
( Gre el (r B)+ ch, ER E1 ) (D) 0, 
(89) 


(enzel +e 2) D) + ch + tt Vor Be 


Die Transformation der Gleichungen (1) auf Polarkoordinaten ist somit beendet. 
Will man die gleiche Transformation der Gleichungen (2) ausgeführt haben und führt man 
an Stelle von (17) die Bezeichnungen: 


AL. 4X) = fü 20 eX) — di | 
IXa tX) = fa, YX, Kal = g, 


ein, so findet man die entsprechenden radınlen Gleichungen von gleicher Gestalt wie oben 
mit dem einzigen Unterschied, dass durchweg e bzw. & durch —e bzw. —s, ersetzt 


emgeet omg "O "w rn belegt D 


(40) 


werden. 
(c) Abspaltung des Zeitfaklors und Normierung der radialen Gleichungen.. 
In den Gleichungen (1) setzen wir den Zeitfaktor zu 


ap(-#) =Z, (41) 
in den Gleichungen (2) aber, zu 
exp + 2: = A (42) 


Fall: (a, x). Aus (36) hat man also nach Abspaltung von Z 


s. lu 0), | 


(43) 
(ann £ +) (10) = oh (= -11 (A), | 
r d ) 
nach Abspaltung von Z aber das Paar 
4- 
Ze? d ed ) 
Lo —. +e) (rA) = chd +; + ) TO), 3 
T dr T | (44) 
b 44 
(s - 2 =u) en) = h (2. ) (ra), | 
Mit der Setzung: 
TA = U(r), | d 
$C = Tun), 
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wird nun (43) auf die Form | 


2 ( ) 
(esə nu = ohf ptt W. : 
T dr T 


2 
(ee oje ang Ho. | 


normiet. Man bringt (44) auf die gleiche Form hin mit der Setzung 


TA = U(r), ` | 


TÜ — V(r), | 
Fall: (b, a). Aus (87) hat man nach Abspaltung von Z 


? (7** Ze -so )(rB) = SE Jen, 
m 


2 ' 
(aan ” +s, )rD) = -eh ƏLİ 


1 


— —_ exem = —Y 


von Z aber 
+ 








T ör 
(+. Ze -,)(rD) = E Je; ' 
T dr r ) 
Mit der Setzung | 
TB = U(r), | | 
rD= —YV(r), 
in (47), aber | 
rB = U(r), ' | 
rD = V(r), : 


in (48) bringt man beide Gleichungspaare auf die:Form (A) hın. 
Fall: (a, 8). Aus (88) folgt nach Abspaltung von Z 


von Z aber 
+ 


4 


(s.- Zei GE = (T LH yo) 
+ d 


l 
TT | 
LG Ze? lg = oh( H1 yea. | 


(46) 


(51) 


(52) 
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Setzt man nun 


TA = vr), 
| (53) 
rC = u(r), m 
n (61), aber 
TA = v(r), _ 
| (64) 
1C = —u(r), 
n (52), so normiert man die beiden Gleichungspaare auf die en 
H 22°) eho = a (2-1 41 | 
T d 
nü R 
Məsə de ağu = che} E- b+ | I 
r dr r 
hin. 
Pall: (b, 8). Aus (39) folgt nach Abspaltung von Z 
) 
(aan TE -s)em =- chd -tE ) eD), i 
T dr r 
| (65 
2 
(aan Z dəə) el a let, | 
T dr ^ ) 
von Z cber 
+ 
) 
(n= EE sen) = el ien, i 
T ar r 
(56) 
; f À 
(1. 22 JD) = “eh + )em. | 
r Mrt ] 
Setzb man nun 
rB = v(r), 
| (57) 
D=-ulr), ) 
in (55), aber 
rB = v(r), 
| - (58) 
rD = u(r) 


in (50), so ergibt sich für die beiden Gleichungspaare die normierte Form (B) oben. 
(d) Lösung der normierten Gleichungspaare (A) baw. (B). 
Das Gleichungspaar (A) 
Uberlegt man das Verhalten der Funktionen U bzw. V im Unendlichen, so kann 
man ansetzen 
U = ə Bis, ! 
m | (69) 
VF = = eS C. pet, | 


r=] 


8—1840P—4 
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wobei 
N (2 = gay, (60) 
= + 
Das Einsetzen dieser Ausdrücke in (A) liefert 


7 (42:6, ete O,retv-1 = ch,[ -a E Bret» 


.— 
ŞB, +v)ret=-1—(I+1) Z Bue], 
` y=1 vl 


g 8 | š 
== (€n + 60) >Ü B,retr+ Zo» B,ret»-1 = ch [ — a> C TRY 
v—l »—1 


(61) 
"m | 
5 + 2,007 y)yretv7i + (l+ U rr] 
Das Zusammenbestehen dieser beiden eege erfordert den Ansatz 
Ban = Denis zs +1) HH DE len so) 7 YZ (sa + 80], | n 
Cri = D, [yZ V (en 9) lov £1) 7 01) V (Ent 99] 
wobei e? 
TE (68) 
Daher ergibt sich aus (01) die Rekursionagleichung 
2[ —yZ.tnt+ V (e$ — 82). (u +v)]. D, = ch ol- —yZ (uv 1) HUET, Di. (64) 
Aus der Bedingung: B, = 0 = : O, = D,, folgt 
P (Dc (+1)? = 0, 
oder u= —14-Vİ(4-1) y. (65) 
Aus Bua = 0 = C1 = Dun, folgt aber 
-yZ . eat V (63-52) . (u+8) = 0 
oder yz A 
condidit (66) 


Das Gleichungspaar (B). 


In Ansehung des Verhaltens von u, v im Unendliehen macht man auch hier die 
Ansätze 


i 
vs . 
u = e719 baste, 
vl d 


t (67) 
vb 
v= e c ret, | 
=] 
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wobei, wie ım vorigen Falle, 


1 
ch, 


Nach Hinsetzen dieser Ausdrucke in (B) hat man zunächst Gleichungen vom gleichen 
Ausschen wie (61) und dann die Koeffizientengleichungen 





V (i si) 


q = 


(ew —&)b,—Choa.c, = Ze?.b,41—Chof(utv+1) F (L--1)] e, | 
(68) 
(En + £5)C; — cho fb, = Zeg o i — Cholu +v+ 1) = (I + 1)}b i. 
Dea Zusammenklingen dieser beiden erfordert den Ansatz 
br = due rv t Yt (Lt 1) len 99 —)yZ V (En~ cə)l, 
(69) 
641 = d. [yZ (ea t) — {(utvt1)—(L+)} (A — ll 


wobei, wie früher, 





əə ch, , 
ist. 

Aus (68) ergibt sich dann die Rekursionsgleichung 

21-72 . ent y (s2— 5) . (u +v)]d, = cho [-y2?+f(ut+y+1)?—(l+1)%]d,41, (70) 
Aus 5,=0=c,=d,, findet man, wie im vorigen Falle 
— y323 (kc 13 — (I+ 1)2 = 0, 

oder p= lt Vilt ++ Z$, l (71) 
und aus ` Din = O = c,+1 = dəki, 


-y2 . ent le) (uty) = 0, 
oder 272 1-1 
En = tof "wa . (72) 


(e) Dis unverkurzten Quantenbahnen. 
Fall: (a, a). Nimmt man Bezug auf (A), (45), (43), (41), (86), (27), (25), (17), so 


(det man 


— ) + | up; + VPA İen 


[va — m)Pr (m+1) _ Vil + m + 2) Pac | gi(mt ie, 
(73) 


— 
ko 
II 

^ 
is 

€” FT MOO 

| 
d 

zər” 

Alm 


: En 1 -m —m |2ım 
V; = exp 21). [upr pam lems, 


p, = oxp( 56) 1 [va-m)p; enə Ve m 2)P; een [teen | 


Danach berechnet man mit Hilfe von [(89); I] die Bahngleichungen als 


144 


app 
gt yı 


— 


gin Ó cos 9, 





2UV 
Ce U'—y? 


— 


sin 6 sin 9. 


Mit Bezugnahme auf (A), (46), (44). (42), (86), (27), (25), (40) findet man aber 


X, = exp[ + 


Sarl, 3 - 
h, T Ss 


Gëss VP 5 en I 


5 1 wel ~m iUm 
X, = oxp( +526) d [ut P, MHD Y(- m 2) p te | med, 
(75) 
X, = exp) . 1 [ver epp Jom, 
h, T 
i 1 " 
X = exp( (rt) d [ut -m)P; (m+1) — V(L+ m 2)pz teen tere 
Danach ergeben sich nach [(43), I] die Bahngleichungen als 
dt U*4 V? 
d: Ui—y’ 
ds _ 2UV 9 
od DIS ps 95 
(74a) 
de _ 2UV . 9 
cın — U: Ps sin U COS 9, 
= Ps sin 6 sin gp. 
Fall: (b, a). Bezieht man sich hier auf (A), (49), (47), (9n), (28), (28), (17) zurüok, 
so hat man 
ý, = exp( 21) i R | - UU -m)Py mt) -V(l+m+2) pi ert) | itay, 
Sn 1 2 
Ya = exp( - si). Elupi"-TPin lem", 
(76) 
EE wer E 1 m ane -í(m O) D-(m -1(m 
"ES exp( nt) al Ul- m)Pp +D + V(L+-m +2) PS 0 el, (m41)9, 
nt I UP-om-..Vp-mla-i:m 
d = exp( - 5^ ) i| pps Pax Je e, 
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und die Bahngleichungen nach [(89), I] von genau gleicher Gestalt wie (74). 


Mit Rückbezug auf (4), (50), (48), (87), (29), (28), (40) findet man aber 


TM ) 
Aq exp + ) . -[| -Utt- mp; oe» + VC m 2)p; (rr Joins, | 
E 1 
X, = exp( +5" ) L|upr=+ pr ‚le emp, 
bo 
(77) 


X = exp( T. 5 t) 1|- U(l—m)Pr + —Vil+m + DPD less) 


X. En | — 11M. m e ime 
X, = exp ux t) "lur, VP "Ë ; 
und demnächst die zugehörigen Bahngleichungen nach [(48), I] von derselben Gestalt 


wie (74a). 


Fall: (a, 8). Bezieht man sich auf (B), (53), (51), (88), (80) (17), zurück, so findet 


m en(- io): 1 [pata uni er | | 
ap tmo 2)p5O D + u(l--m) Py ma) [tnt ae, 
(78) 


— T —. =m t 
vP Yb, Je "ə, 


Va = oxp( Zi) : 
; 55 
W = e(t): 


š I|-vt-m- 2p Gn 1) — 4 ul m)p; otn jer, 


d 


İl 

© 

Gel 

ve 
” 

| 

is 

— 


Demnach berechnet man mit Hilfe von [(39) I] die Bahngleichungen als 











d 2 3 cos 6, 

| 79 
— EX "a Sin 0 cos 9 T 
cdr x : 
ay = — 2vu sin Ó sin 9, 
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Mit Bezugnahme auf (B), (54), (52), (88), (80), (40), findet man aber 


upr” leng : 


-olim +P P+D ulm p (+) Joint Dr, 


Ba 
|! 
E 
ATN 
+ 
|a 
+ 
Net 
~ | = 
rm rm 
a 
c 
ane l 
Es 
E 
Wal 


—— —-—.—- — — —  —— — ——” 





80 
X, = exp( +524) i L [oppg tup” Jems, (80) 
X; = exp( +t) a | - v(l+ m, + 25570 + ull-ın)Pr (m-F 1) Jomo, 
Demzufolge ergeben sich nach [(48), I| die Bahngleichungen als 
di  v'ru* 
dr 02—u?” 
de 2vu 
cd” = er COS 0, 
(79a) 
dz 2Uu . 
ej “xain Ó cos 9, 
dn. Se 





U . 2 
„sin sin 9, 
-u 


Fall: (b, 8). Eine Bezugnahme auf (B), (57), (55), (89), (81), (17), hefert für diesen 
Fall zunächst 


d = exp -tnt ° [dem +a)P (m+) +uli-m)P; nl... 


(81) 


y, = oxp( ai) ap: “uv” le. 


? 1 -m — im —IUmn 
" Séil dE +1) Ali «e (neo, 


lI 
E 


Ya 


İl 
© 
“3 
| 
Ei 
"at 


. 1 [Pas e pz]; a j 


und dann nach [(89), I] die Bahngleichungen von genau gleicher Gestalt wie (79). 
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Mit Rückbezug auf (B), (58), (56), (39), (81, (40), findet man aber 
Ki = oxp( ent) ` - h + m + SEH — u(l == m)Pr (mtn 1 onan, 
0 


Xa = en + 52t) =| oPn tu Pe -i(m4 1g, 

| (82) 
ef 1 
X = ex m6): 2 ptem mese stor) emen 


d En 1 — -m kel img 
X3 = oxp( +524) > - [vPas- up, |e t ; 
und die zugehórigen Bahngleichungen nach [(48), I] von derselben Gestalt, wie (79a). 


(f) Die entsprechenden klassischen Bahnen. 


Wir beziehen uns auf [(B), I], setzen dort 
Ze 


e=—=V, As z 0 = Ayk: (88) 
und so bekommen 
EHE EN 
ə C ° cot, ro S rainddg/ J’ er 


Mit her Setzung 





KE —et+w(n), (85) 
hat man daraus 


ege (e zy "Se (86) 
und nach [(A), I ] die Bahngleichungen als 


2 
pu 95 +e V =(~2+52 ) 
dr at 


(87) 
= a — >= cS’ sin 6 cos 9, 
— £j E = oS = 08" sin 0 sin ə, 


Mit gleichzeitigem Umsehlagen des Vorzeichens von e bzw. a, hat man aus (84) 


2 2 2 
(ET N eg (88) 
C coi c g O zT 
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und nach [(A), I] die Bahngleichungen als 


di 98 gen EE 
Um maa” 
da oS , 
— — 0— == 0 
orm ox cS! cos 
(8v) 
00:00 RT 
dy 98 I g: inu 
0 = — 0 , 
5 "E ^5 cS! sin 6 sin 9 
wobei 
S = et-wfr), (90) 


zu setzen ist. 
Setzt man ın (87) gemäss (85) u (86) 


O! = + +2) «a= +w’, (91) 
aber in (89) gemäss (90) u (86) 
——. : 


so ergibt sich die gleiche Gestalt für die klassischen Bahngleichungen (37) bzw. (89), 


nämlich, 





, dt, Ze 
° edr PC 
d d d | D 
A ; I : ! 
a U Lo (cos 0, sin 0 cos 9, sin 6 sin 9, | 


(oi Die Quantenbahnen bis auf die erste Náherung. 
(i) Reduktion von (A) bis auf die erste Nàherung. 


Leitet man aus (A) Gleichungen zweiter Ordnung in U bzw, V her, so lauten sie 
el $, -2 Je an 22 Jo | 
0 3 r ' 0 


dr? T 
2 
Ze LE +1 İr, 


m. — En + Ze?jr) + e]... dr m 
- / - ` 


(94) 
an| 7. - Cen 2) İr-İ-e- Zr) - 3] 
T r ° _ 
Ze? i [+ -— 
“428——7.———.65h5—b— 
r? [(— €, + Zo* [r) — s, | dk dr” T = 
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Wir schreiben sie dann. wie es folgt, um 
I 
r? r or 


ana AO | | 
r*[( — En 4- Ze* [ry * əş) LU rd 





, (95) 
ozu. 
°L\V L r? T E 
| Ze? : u zə 
_ ne Ze p= l: 
| gt est Z6?/r) — s] Gage V o T 
Dabei bedeutet der Strich eine Differentiation nach r. 
Zwecks Näherung machen wir nun in (95) die Ansätze 
Lim(h, In Uy = Bit how, Ela d uas 
h,-ko 
5 (hə In V = 8’+ lg, + h38, + ee | (96) 
` Lim [k(l 1)] = a0. | 
h.-*o J 
Man erhalt dann 
o* [I (S" + hoa, + hia, +...) + (I cha + hig 4... -01 
= de ZN "dun ah T(S’ + ha +h? 
-İ( “ür” ) E P[C-«*Ze]r-«] d — 
5 m (97) 
c^ [h,(8" + hj, + Die + ...]) (87 + hg. + həə” 0]. 
Z Ze’ \ _ ee ah Ter + hog, + h28; +0 | 
Jon =) E r°[(— En + Ze*[) ^ o] .— SR | 


Daher entnimmt man bei nullter Näherung 


eS = İ( "m Ze’ y- al, (98) 
: T 


bei erster Naherung aber 


a en 


Q İS” Ze" 
B, = [S zen 


(99) 





mn en 
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Greift man nun wieder zu (A) zurück, so erhält man im Limes (h,—0) mittels (96) 











(os SÉIER =| en] + i+] TEN] = cfs’ hag, + Mie 0), | 
À L "Sp (100) 
C En + Ze? heliz) =[en (ə H = CS’ + how, + ho, ...) — 0]. | 
r U U T J 
L L 
Bis auf die erste Näherung hat man daher 
D , 
Z U V U. V ` 7 d , I 
(m — U bə v) = 20S" + chol, ba), - 
6520022 
L 
"AR Z8" “< 
Aus (99) folgt : BL” 4 Giga (102) 
und mit Hilfe von (98) noch 
Demnächst hat man 
U y 2 ey , oh, Ze? EN | 
Le: el "ga te 7 yes dido } (104) 
U YN 2 | ‚_ch, Ze? s,(—s, au | 
(v yar... 7 çiğ" | 


(4) Reduktion von (By bis auf die erste Náherung. 


Folgert man, wie früher, die Gleichungen zweiter Ordnung aus (B), schreibt sie als 


vi) «(or tte- ET 


8 Zei el E-H | 
rT (—s_+Ze3/r)—«,| Lù r J | (105) 


ef) «(y - mentem [Cus te y ca] 
v + r* T 





RES CNN (GE + i+] | 
Pl, Zerre] “Lə ol 
und zwecks Nöherung darin setzt 


Lim(h, In u)! 2 8' + hoy, + hy’ 
h,>0 


S Se? 


(106) 


Lim(h, In v)’=9’ + hyd} + A38; + .... 
R.—0 

Lim[h{l+1)]= a0, | 
h, +0 J 
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so erhält man 


c* İR (ST + hoyi + hey, +...) + (8/4 hoy, + böy + ...)?] 


0 








xı. Ze? )- ] ET Ze? : 
I ent S e Ile Ze] xs] c*h,|S "ep? Moya us] 


" . ; : (107) 
c? [h(S +8, + hiği +...) + (S + hod, + h28, + el 
a | _ Ze? ) j]- Ze? ; D 
(rm ze ar, 


I 
Daher findet man bei nullter Naherung dieselbe Bedingung (98) oben, bei erster 
Näherung aber : 


dE EEN Zei Ce 
| = ə ls 


; S" Ze? 
SÉ ran] = | 
= Us m{(—e,+Zetjr) rell 4 


(108) 


———— —— — 


Aus (B) erhält man im Limes (h,—0) mittels (106) 
2 
(7s + - )-«(2) = (S° + hoy; k RI, +...) — 0], | | 
L 
} (109) 
Ze U € H 
dE = pes) = e[( + Iu, + 138! +...) + OF. | I 
T D L : 


Bis auf die erste Näherung hat man daher 


Ze? d d v d n i 
(- «+ e. E + JA i: vər z) = 2¢5 +ch,(y, + 6; ), | 





L 
(110) 
2 
25xX5508 
T u v d 0 1 1 
L UL 
Nach (108), (102), (103) gilt aber 
13; * 9i =(, 
_y—_ Zei s (111) 
Ty 1 7? 0392 
Demnach ermittelt man 
van. 2 f/_, Zeche h Ze: oe, ] 
zal + Je ml 


(112) 
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(iii) Die Bahngleichungen bis auf die erste Nüherung. 
Die Fälle: (a, «) bzw. (b, a). 


Ist der Zeitfaktor oxp( —i8t) =Z, so sind die unverkürzten Batngleichungen die 
e 


durch (74) gegebenen. Bis auf die erste Näherung haben sie also nach (104) die Gestalt 





dt : Eus 
dr (y V e , h, Ze* &( — En + Ze*/r) 
EB 485. 
y Ü : 9 p? ew 
dz da dy )= en 

Lë, leng egen on 

f , : ci? 

(cos 0, sin 6 cos 9, sin 6 sin 9) - ‚_h, Z6? e{—8,+ Ze'[r), 


ctw’? 


indem gemäss (91) S’= +w’ zu setzen ist, 


— 
” 


İst aber der Zeitfaktor oxp( +/2t)=2, so sind die unverkürzten Bahngleichungen 
0 + 


die durch (74a) gegebenen. Bis auf die erste Näherung haben sie nach (104) also die 
Gestalt | 


? x 
(> 2 ES SE — u^ 
dt 4 








H U/, 293 hy 
de U zl °: E h, Ze? "İsə + Zen) 
So = + (cos 0, sin 6 cos 9, sin Ó sın 9). — r (114) 
Chr car T 65” 
Le U /L 
. ; cw 3 
=—(cos 8, sin # cos 9, sin 0 sin 9) 7 000 
WU Kir Ter 


weil gemäss (92) S’= — w gesetzt wird. 

In nullter Näherung stimmen also (118) und (114) miteinander und mit (93) überein. 
In erster N&herung weichen sie aber voneinander ab. 

Die Fälle: (a, 8) bzw. (b, 8). 


Im Falle des Zeitfaktors : oxp( - et) =Z sind die unverkürzten Bahngleichungen 
0 
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die durch (79) gegebenen. Bis auf die erate Naherung haben sie nach (112) also die 
Gestalt 


; (29 5) Lan 2 
de x0 en a i 
T (ss) ande t ote tein 
u v ° 2. f rw" 
L | 
TES )- | 7 2 \ (115) 
Ka ee | 
d V/L 
. m M gu 
= — (cos 6, sin 6 cos 9, sin 6 sın al, ho Ze? &(—&,+ Ze”İr) 
| ular 
) 


wobei gemäss (01) = +w gesetzt worden ist. 


Im Falle aber des Zeitfaktors : oxp( 720) —Z. sind die unverkürzten Bahngleichun- 
0 + . ` I 


gen die dureh (79a) gegebenen. Bis auf die erst Naberung haben sie nach (112) die Gestalt 


2 3 2 
(242) (za = u əəə 
T 2 P gur” 
30 1-0 MN 
dr (2 Si ‚ho Ze? e(—:,+Ze?/r) 
DE E ee 
9 T? c^w'? 


( 2, SSE )= + os, sin # cos 9, sin 6 sin 9) Š 5 
Ë 1 v/L 


= — (eos 0, sin 6 cos 9, sin 0 sin 9) 


wobei nach (92) S'— — w' gesetzt wird. 
In nullter Näherung stımmen also (115) und (116) miteinander und mit (93) überein. 
In erster Näherung weichen sie aber voneinander ab. 


fak?) 


` Das Elektron im konstanten elektrischen Felde nach der Fassung (i) 1B,, A,l bzw. 
(ii) [B,(a), A,(a)] des Teils I. | | 
“ () Fassung [B;, 4,]. 


Mit dem Ansatz: ] 
p = p- Es = V, A, D = À, = Ae (117) 
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hat.man' nach [(85), I] 














y ijv 
IG 2 27 EM Cep I = <2 i , 
cot C C C Ze 
A _ n Və ; M (118) 
+ ANE A 
[eg Tch] A | | 
SS d Ya S Və 
Mit dea Setzung 
, , TES — et 
(dy, Va; Wo V) = (f, g; f^ g^) exp —— m9, (119) 


ergeben sich daraus 


d” (eot eB LJ? )( SAMIR 
pa ül cah, ch, ch, T op İ om 


\ (120) 
d? [SE eE } eE y s” i i 
Z 345 TOA 4 65 ch Ba = 0, 
İz PE DNE ch (6 (z+ eH /17 
ef eE on, di | |- | 
al oA (a4 SE JEt heit. ip)g | f, 
(191) 
ci eE e^ |, dg ; |- , ' 
el s («+ 2E y hog, t (P1 +ipa)f d, 
wobei zur Abkürzung gesetzt werden 
6 — 69,5 = e^, 
A | (122) 
Pit p = RT, 
Mit der Setzung 
GE e,7+¢7K? PS 
24 128 
(2n) ch, eh ” (1281 
sieht man an (120) 
f= dën ($), gv bes (€), ] 
} (124) 





f eB. WM s” ) 
(et) j 
WO 94—;. Pn die bekannten- Hermite'schen Funktionen sind. Mit Hilfe bekannter Eigens- 


chaften von Hermite’schen Funktionen folgert man nun aus (121) leicht f 
f' = ale, g'= b'es .(£), (125) 


a! = ef- aN Eo TE 5— + b(p, SCH | "m 


b’ = =. vet IR) E ip.) 
H 
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Berechnet man nun die Bahngleichungen nach [(89). I], so findet man für sie 


di A/ (s,* + c* K?) Pn-1 Pn ) l 

—— — — | — + —- |, - 

dr 25, M Pn ?n-1 
- dz _ (e cR!) (fan — 95 ) ! 

cdr i UM Te Pn 5n -1 I I x. 
dn ep E 

cdr e, = CAT 56 : 


Die klassische Bahn. i 
Wir beziehen uns auf [(B), 1], machen dort den Ansatz (117) und to bekommen 


«e! (dS ey = ( BF 
PU TI AST s 
Mit der Setzung: 
S = —st+ep tun, t p.(9), (129) 


hat man daraus 





33 y = ə 87033 et o Kt. usus 
Een 


Nach | (4), I] lauten die klassischen Gleichungen dann 


dt ƏS x) 

020. V=- "m P)! 
cn ait? (e E) tB 
_ s dz 098 — . _ s dy _ OS _ | 

cdr ox — Pi> cà dr ay Pa (130) 
satr d )( — T-S. dp, — | 

c? dr Öz + C ES eH c^ — da | 


Reduktion von (120) bis auf die erste Näherung. 


Schreibt man (120) wıe 








IN?’ 2 2 2 € 
ORO 4 (Eye rn | 
(131) 
d L) «(£r LE. E -( E Tt < y= 0. | 
und setzt dort zwecks Nüherung i : l 
Limes (hy In fi” = S” Abee, that... , | 
kt J (102) 
Limes (h,In g) = S+ hogi 4 hə 5... , 
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so erhält man 


hi (S + hoa’ + hat +...) (S har BCEE mE 
-h, m ees Yə — ye 0, 
C C eE 
h(S! + hg; t RB +... + (S hop, + h38; +...) + 
did a d 
C € 


Daher entnimm man in nuliter Nüherung 


2 2 3 E s” 3 
c? C 5 eE 0, 


in erster Nüherung aber 
oe S" eE 1 | 
Her vB) 





(188) 


(184) 


(185) 


Reduktion von (121). Aus (121) hat man mit Rücksicht auf (124), (125), (126). 


hoe n gel ok (2+ T )= ız N (e! + c" KC) : CR 
eH C 9n-1 
h, d fine) + E 6E E (e+ 6” )- " V (83 -Fo*K*) Fna f 
"da eE C 94 


Daher hat man in Limes 


— Mie eR) (a ) = (8+ hoa; + bu +...) 
m 


L 


als 





2 2771 y 
+ M(t? +e?K | (acs) = (S° + hoB, t h38, +...) + en (zs - 
d Pn L 
und bis auf die erste Näherung f 
277% , 
Vo FOR) dE e ) = da. -ai + slk — ) 


c 95 $n-1 L 


Vize + c* KT) (e z $e) e de, ee) 


C Pn $n-1 L 


(186) 
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Die Bahngleichungen bis auf die erste Naheruny : 


Nach (127) lauten sıe wie 


SE ( s” )- 7 ] 
du pesar a wl 











a wl 
d 
da | , hə id 
— = —İlU ——. — 
I ch 
cdr! Se 2 Ww 
dz 7 dy — ə 
2... 0. 
: cdr ee cd x 


woLel 
ETS c. b en (+ ir - "Kb = dps 
Qus eges (ell z + LE — + ; 
Bei nullter Näherung stimmen sie, wie man es merkt, mit (130) überein. 
(ii) Fassung [B,(a), 4,(a)] 


Mit dem gleichen Ansatz, wie oben, nämlich, (117) hat man nach [85(a), I] 


























X X 
IO? ez) hiy" Phe oE | el |, 
cat C c? ix | 
Xa 2 ! 
X á | 
X I 
[a evn] |= Sy | 
š i x = Xs š X2 
Mit der Setzung 
et —25,— 
(Xu Xa? Xp Xa) = 595 f. g) exp Fs i 


ergeben rich daraus 


d? e, +c? K3 = deit: 
Lis d erh, uic C 


s” vi 
= 0, 
0 es (2+ eli ) | İ 
d* ec K? ob -( eE WÉI (2+ > )] =D 
| +Í ch eh, ch eE d i 


0 


[Et SYandl-o-mo]er, | 


0 











en ee” 


-4 = (2+ tL y- ho 39- (p, +ipaf]= g’. | 
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(189) 


(140) 


(141) 


(142) 


(148) 


(144) 


v, K haben dieselben Bedeutungen wie in (1929). Nach (143)w (144) findet man in 


Abhängigkeit von Hermiteschen Funktionen Prl), en-i(£) 
6—1840P—4 
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f= a9. (£). g = Be, i (E), 


£ 


j = Pnl gi Bond), | 


555 





} (145) 
ER C | (2 + c* K?) b 
a = ——} ———— a- (p —ıp))$ 
Bo 
` I : C 2+ eK 
Monte b 
wu gilt wie früher 1 2 
(2n) eÉ _ s tok 
ch, eh. 
Nun berechnet mau nach [(48), I] die Bahngleichungen wie 
di _ V CH) (a 1 + 94 ) 
dr | n-i | 
_ dz __ V (e + eK?) (===. -»_) 
m Ve ma^ | dd 
` | 
0 
cir ağ” cd s Pa 


Die Klassische Bahn, 


I 
Mit Bezug auf [/B), 1], aber nach gleichzeitigem Umschlagen des Vorzeichens von 
e und m, erhält man ınıt dem Ansatz (117) 





s! (08 eV y- ( as y 
ci (S (c | 2 öz / ` (147) 
Mit der Setzung | 
S = st— 2p, — yp, — p. (z) (148) 
hat man daraus | 
n )«( ər 1 77 = ga 
E : a+ E 2: m 8”, (149) 
Nach [(A), I] lauten die klassischen Gleichungen der Bahn 
dt 98 : e^ 
cp 0 expli bod " 
7 ar ° +oB(? + SS } 
eds. OS —— “dy oS —— 
Pde wc T. : ES 


| 
so dz EE MN eE ) e’ E od 
cd: 88 vi 7 (+ | c? I 
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Reduktion von (143) bis auf did erste Näherung. 
Verfährt man zwecks Nüherung in dem vorigen Fall analoger Weise, so findet man 
zunächst noch das Gelten von (134) bei nullier Näherung, bis auf die erste Näherung aber 
Limes (hln f)! = Limes (h, In Pa)! = 8’ + hg, ) 

h.—+0 h,+0 . f 


Limes (h, Ing)! = Limes (ho In g, = BS! + hor | (161) 
h,+0 hR.+0 


J 


wo a!,, 8”) dieselben Ausdrücke wie in (185) haben. 


Reduktion von (144) bis auf die erste Näherung. 


Aus (144) leitet man im vorliegenden Fa!le zunächst (130) noch als gültig ab un 
daher bis auf die erste Nüherung auch noch (188). r 


Die Bahngleichungen bis. auf dic erste Näherung. 


Nach (140) lauten sie wie 








| 
ds | c[, hə vil | 
cdr SE 2 wP | (152) 
_ de o dy) _ 
57707 bo. 
wobeı mit Rücksicht auf (184), (150) f 
"E Í i je 2 L a2 7 - dp, 1 
w’ = S= + = vic) a+ — (se +K?) p= tu (158) 


gesetzt worden ıst. 

Man sieht, dass die ach hieraus nach Abstreichen von h, bei nullter Näherung 
ergebenden Gleichungen mit (150) oder (180) übereinstimmen ; die Glieder erster Ordnung 
rechts in (152) sind aber gegenüberden Gliedern erster Ordnung rechts in (189) von 


entgegengesetztem VorZeichen. 


% Das Elektron in konstantem magnetischen Felde nach der Fassung (i) LB,,Aş1 
bzw. İB, (a), A,(a)] des Teils I. : 
(i Fassung [B,, Az] 























Mit dem Ansatz Q -0—A, 
(154) 
j A, =—4Hy, Au =+ MHz, : 
wo H die konstante Feldstürke in der e—Richtung bezeichnet, hat man nach [ (85), I] 
20? (, 0, eH y-( 9 eH )- 2 Q' , rel, Hİ yi o. yi 
k (beget E) n Gu tl n UR 
f . A n at V (155) 
3 (no + Ey). o hee SN °] - b” | 
h h—— dE )+oh,— = 2| |, 
| ” cat Hol GM . "Əy 2c 2055” | € Me | 














l 
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| 


Mit der Setzung See . 
(Yis 33; V, y, = (f, g, e WË P 


hat man hieraus zunächst 








og af 92 SE : oH (22-2) 
E zi 27 By ie (x? + y?) + he, = ER Me 














| H ih,eH | İ — "o 
uu. Pura 
| g 
! 
o cH İ Eo İ 
Lalin x leen AJ" ele 
dann aber in Zylınderkoordinaten 
gj a? 198.8, a CH aL eH 9 | 
N g lo? rar ro) Ach? | ch, ch ' g 
DEE | 
mit K` = — — Ian | 
"E us 1.0 ie __ | 
ere -( $ pa) +h ée Sëch Seen, 1 


Setzt man zwecks Lösung weiler | 




















f = au(r) exp| (m— = öl g = bvi(r) exp| (m+ iy] 


so ergeben sich aus (158) 


F° =| d 1d „ mn — AN a2H2 2E 
r 


717 
0 











K rd 1 d ,(maiy eH? = eH 
mrt igt ir Mero ch, 


und aus (150) mit dem Ansatz 


[i Br e 
C; i 


dr r 





| 
e 





(: +p,)+aK’] » NEN 3 


| 
| 
) 











wY.  — — 


N sten (l ae 2 


— 


(156) 


(157) 


(158) 


(159) 


(160) 


(161) 


(162) 


(168) 
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Aus (162) folgert man leicht 


d? 14 „(ms eil" 3 eH ( i hb 5 K^ u 164 
E dr? ə. r dr "n on. E 4cth,? ch Wen 
Somit setzen wir mit Rucksicht auf (161) 


KE: = K, (165) 
und weiterhin 


+? = d. (166) 
Zusammenstellend had man also endlich 


jeu exp| (m - Zei g = bü expl (m+ 5 ye], 


f = au exp|(m- Yel, g = b”ü exp| (m d y], ; (167) 
2 2 | 
| 
e = £[—a( £-p,)+0x], vw: +p) tak]; | 
66 C Eo e 
Berechnet man nun nach [(89), I] die Bahngleichungen, 
80 findet man 
250... 
dr £o I cdr Eo Ps» 
dz 
c mE A = _ 168 
B 2 ( exp (i6) E = exp ( ON | (168) 
Sab 8 
x. nat u 
7 P 1 exp (16) = exp ( — 18) | 


Die klassische Bahn. 


Nımmt man auf [(B), I] Bezug, so hat man mit dem Ansatz (154) 
TEE 
c? cot Oz 2c Qy 2c Qs 
27522 LL 
- (š 344 ay aoe No, “ər, 0: — 


-(8 Y-(Gy-((Sy EE elt 
cat O2 ‚Or 138 4c? c 9601 


nach Transformation auf Zylinderkoordinaten — Mit der Setzung 





= —et + gp; +4 0 +w(r), (170) 


a8\? _ 2 a oH oe (> a 171 
OB age qe” ` “əl an 


hat man aus (169) 
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Nach [(A), I] lauten die klassischen Bahpgleichungen wie 


"dr ^ Oi — 7 dr 9, 


x. 
La a 
o 
Co 
SCH 
CES 
s 
o 
TA 


= t” cos e-(* — oF) sin 0, (172) 
T 2c 
cdr Oy 2 — 
= w sin 69 ( £ — Sr) 056, 
T 2c 


Reduktion von (161) bis auf die erste Naherung. 
Schreibt man (161) mit Rucksicht auf (166), wie folgt, hin 











gp = ary «(any u, hm! — im —3) eH a eH h (m & a | 
u u r u p 4c? C 2 1 (178) 
7 
m= (JE EA əəə mb ar S Mm- i) | 
ü ü T ù ? 4c C 2 
und íuhrt dort zwecks Nüherung 
Limes (holau) = BI the, + h2a +... , | 
Limes (R, ln tt)’ = S' + haa, + as? . V (174) 
ht | 
Limes (hm) = a, | 
h,+0 
ein, so erhält man in nullter Näherung 
2FIa 
S^ P Lm p "uui Ed. (175) 
e? C 4c? 
in erster N&herung aber "nm les eH J| | 
utar S' SM 2% Il" 
; (176) 
I2 
Reduktion von (162). 
Aus (162) hat man mittels (174) 
"X". ah del, ge | 
(S'+hça + hiaat...) + x + -—. Wd 
(177) 


0 = 1a 
(8 + hoai + hag T 5 


; 
hk, ieH u 
0 pe J 
25^ de ` ü | 


VHRSUCH EINER LOGISCHEN QUANTENDYNAMIK DES ELEKTRONS 163 


Daher folgen im Limes (h,—0) bis auf die erste Naherung 


K( È exp (16) + € exp (ie) = de COB dit — = rom e| | 
ü L 
(a+ (eH /2c)r*) (® eH ) pe | 
Aa? - 207! cos 0 + S’ sin Gİ, 
(178) 
ü u SR a eH 
-ıK = = _ = 215 0+[— — — ) J 
GE (18) = exp ( al d sin (2 gg TI sos 
(a + (eH [2c)?) | j SE eH ) l | 
het DEAD] o [s cos 6 x 25 sin 6]. 
Bis auf die erste Näherung lauten also nach (188) die Bahngleichungen als 
75:55 
dr e, cdr sn 
= da = |. cos 0 — ( $ _ eii r) sin d 
Car t6 T 2c 
(a+ (eH /2c)r*) iie. "eH ) NE | : 
P 9 uL | a cos H tu EH ; (179) 
- 4 = “İrsin 6. ( £.— 2E +) cos] | 
cdr Eo T 2c | 
(a+ (eH [2c)r^) [wc -(2- 8 Jsi | 
2 w! cos 6 777 sin 6]. | 
Dabei versteht man mit nakahi auf (170). (171) 
2 2 2 
# em K? ze. aeH eH > as 
w’ = + ( = + I Ac? T + ör (180) 


Mit Abstieichen von A, bei nullter Naherung gehen also (179) im (172) uber. 
(ü) Fassung |B,(a), A,(a)] 


Bezieht man sich auf [95(a), I], so bat man mit dem gleichen Ansatz (154), wie oben 














iz O - (hoe - 22 y) - (2 9 , eH ə Cap ^ xin is 
ia Or 2c "By 9c Oz? o X, R . (181) 
X A. 
9 8 A y)+o e (hs. oF e) 9] ss 
E cot te bas Bo YTA agt Se J BEI | ç |x 5 
2 2 














Mit der Setzung 


dı 7 , t— 
D, Aan, Xi, Xe) = (f, g; f. g Joxp— Es (182) 
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hat man daraus zunächst 



























































B zi SN gë : H (ze Ba m eg Hİ f e? f | 
Ste Eege )— "ay "əz SCH 
C Öz? OJ € C d 
: } (183) 
, f Pr 
HLC E DL Pa ep İİ |= | 
gl" gl 
dann aber in ..—.—.—..—.—.— 
l 
E İİ rə LƏ. 9 eH? , eH ð eH I | 
h “alar ag" Ah Ich, Br ) 
g | 9 \ (184) 
mit K? = 9 oy | 
C 
56 Ce "E o t Si ie H E ) 
— — = | — — h 10 f = — — t8 
e ( nh: S t r 36/9 ə | 
\ (185) 
$y , 6 8 i 2) 16H 
— ° g” —İ — heil — + - — ar. 
SÉ (+ +2)0+ car ES 26 ` rf | 
Setzt man zwecks Lösung m 
f = AU0) oxp| - (m+ 2 )l g = BVG) exp| -(m- Zei | (186) 
so erhält man aus (184) | 
| 
Ke fd? ld (mag? eH i oH . | 
in? ‘ante + ma a | 
| \ (187) 
K° SE 1d (m—-di)? , eH*', eH/ 2 Iv [ 
hg ES EEN qu p x Aer 3 ch, 77 2 i 
und aus (185) mit dem Ansatz | 
VV E y- E en | 
°Ldr T 2c ; i 
t (188) 
pf 2 — “edil Ju. 2; rU = ET, | 
"dr T | 
s 5 f i ) 
arain) pl-i] | 
(189) 
-19-İB(i +p) xac] oxp| - (m- el | 
C C | 2 J 
Aus (188) folgert man leicht | 
j i 
d (1.4, (m+ eH? , eH (m .ı ) ED: 
İz 2 "r (e ux ch, 2, S i 10) 
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16 
wonach man mit Rücksicht auf (187) 
HIE (191) 
setzt. Die Gegenubersteliung von (188) mit (162) ergibt fernerhin 
Uu, 
(192) 
H cu | 
Die Losungsergebnisse slellt man demnach, wie folgt, zusammen: 
f ss AA exp] - (m + = bl, g=Bu expİ - (m - 5 a], | 
fi= Ai oxp| - (m+ 31! gi Bu exp| -(m - el | 
| (198) | 
Ar — “la(: -p,)+BK], i 
Š, C 
B’ = SE +p) +AK]. 
8, C 
Berechnet man nun nach 1(48), I] die Bahngleichungen, so findet man 
2 
dr s, cdr =s Pe 
dr ü ; d 2 
25” Es vi x exp( — 16) t z oP G9) ) t (194) 
2 — Hİ Ž exp (—10) — — exp 00) 
cdr 25, 


Die klassische Bahn 


Mit Bezugnahme auf [(B), I] aber nach gleichzeitigen Umschlagen des Vorzei- 
chens von 6 und m, erhält man mit dem Ansatz (154) 


v. -( 38 (8-8) 
c! \ cat da ge” Oy 2c Oe 





(Yen EEE] 
“(ər rer “əy P 


— (195) 
C Oy “öz öz 
_{ 88 y-(2y-[E 3 E + a, 0H 25] 
“(zər öz ar) "Lag Zei "e gel . 
nach Transformation suf Zylinderkoordinaten. 


Mit der Setzung S = et-zp,-ad-wifr), (196) 
hat man aus (195) 


oSv, c ,aeH eH? GE 
——0 





T C 4c? 
6—1840P—4 
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Nach [(A), I] lauten die klassischen Bahngleichungen wie 


di og s dg 098 
—E,— “niz (r olm E ces pag 
dr at c? dr g ! 
ge dz OS eH e ($ eH ja 
M uda oun PEN i= = 855 6, 198) 
cdr Oz 07 — r oo J t l 
€o dy — OS, eli, = ur sin 6-( 2-32 rose 
cdr öy 2c ? 


Reduktion von 2 bew. (188) bis auf die erste Näherung. 


Man sieht nun, de: wegen (192) alle Gleichungen von (178) bis (177) fur den vorlie- 
genden Fall verwendet werden dürfen ebenso wie für den früheren. Daher findet man nach 
- (177) im Limes (5,0) und bis auf die erste Näherung 


d = 7 got «9l 8! a _eH ) | | 
geb exp( —10) + = exp EAR de Gos di 2:0 sin 6 
a+ (eH /2eh" a 6H , 
— h, on ; EE -5S sin 6], 
| (199) 
[Üü ; u : = "wap. f & “eH 
iK( F exp( —18) — 7 exp((6) = de sin 0 ( 705 r cos d 
— h, t (eH [2e)r* l IER eH r) j | 
— PS — 2 € DE sin 6+ ST cos 0 
Setzt man nun mit Rucksicht auf (196), (197): 
H èH’ as 
wi = + Neo ee el = — P 
e 403 Or’ SH 
80 lauten nach (194) die Bahngleichungen bis auf die erste Näherung als 
(db s ds c | 
T 6 wd əə | 
= dr — Ze cos 9 - (2 _ 6H r sin d 
cd & T 2c ` 
` e a+(eH (ein? iH : | 
rA DE JE = E o r) cos Ü +W TIR d (201) 
_ dy _ ölə: gin dE eH r )eos d 
cdr & T C 


o | a cH 
+—} SE 
S GER ( 2. sin 0 — w’ cos d 
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Mit Abstreichen von ho, also bei unllter Náherung, gehen (201) in” (1£8) oder (127) 
über ; aber die Glieder erster Ordnung in (201) sind gegenüber den Gliedern erster Ordnung 
in (179) von entgegengesetztem Vorzeichen: à 

Das feldfreie Elektron nach der Fassung (i) 1B,, As] bzw. [B,(a), A,(a)] des Teils I. 
(i) Fassung [B,, 4,1. 

Mit dem Ansatz 




















— 
— 


f o = Ü = Ag == A, = As, (202) 
hat man nach [(85), I] S I 
| q e iv] " Se Vi 
0 c at? 0 d c? d L 
; (208) 
++ A g d, 
b. » + holy) | Eu. 
si le 
Mit der Setzung 
a , h — ti typ. 
(Yi, V. 5 dn, Va) = (a, b, a’, b’) exp Lp WP s (204) 
` 0 
hat man daher 
2 2 
s “pr para) =i (205) 
a! = °[-(> “pə Ja (p; — ipg |, | 
° \ (206) 
b’ = SIE + pa Jb + (Py + ipaa |. | 
56 C J 
Nach [(89). I] findet man die Bahngleichungen als 
di e ds __ C | 
dr = £5 cdr D 3? 
r (207) 
dr oe dy | 2 | 
cdr CN Pu cdr = e Hi ) 
Die klassische Bahn, = 
Mit Bezugnahme auf [(B), I] hat man b 
sl _ 98 Y- GI 208 
o JS Na) |: | ES 
- Mit der Setzung = — st +29, ka 3D,, (209) 
folgt hieraus 3 e b f 
=. — (p+ pat p3), 


c? C 
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und nach 1(4), 1) lauten die klassischen Gleichungen dann 


dt ` 5. dz 
de P957 775 gir" | 
__n de vədi | 
Pı 8 ee c cdr J 


Man merkt. dass (210) mit (207) gleich lauten. 
(i) Fassung | B,(a), A,la)]. 


Nach [35(a), I] hat man mit dem gleichen Ansatz (202), wie oben, 














A 
2 3 1 
bəz” hiv? | = ei , 
: l x, X, 
N AN: 
|^. a + holy) | sad d 
X, v X2 














Mit der Setzung 


(Xi X, ; Xy, Xs) = (a, B; ar, B") oxp Pc Pa əbə, 


hat man daher 


2 
Eo 
, 
c? 


— 


g? o° 9 2 
m —(pi- pat pa) = 


g = — Ile + pa JB = (Pi + pada], 


£o 


Nach [(48), I] ‚ergeben sich die Bahngleichungen als 


dt s dz ` 2 
1000007 Er mM NINE F 8 
dr ge cdr ... 
de _ p dy _ cy 
edr o edr 5” 


also von gleicher Gestalt wıe (207). 
Die klassische Balın. 
Mit Bezugnahme auf [(B), I], doch nach 


Umschlagen des Vorzeichens von e, hat man wieder 
$- [89-367] 
ei ` le z SCHÉI 


rt Na Sa 


Mit der Setzung 


6 


— — — pan S aana a e 


w et ne, 8 E E n aae 


(210) 


(211) 


(212) 


(213) 


(214 ' 


(215) 


` 
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folgt hieraus 
2 
e = ^ — (pj? + P3? t Pa h 


M 


und nach [(A), I], doch mit Umanderuny des Vorzeichens von &,, luuten die klassischen 
Bahngleichungen als 


dt _ & da ] 
"gg coe Pa Bede 
(216) 
"ur de p =° du. 
EK e cde ? e cdr 


Es lauten (216) gleich mit (214) wie auch mit (210) und (207). 


6. Sohlussbemerkungen. 


Am Schlusse wollen wir die Aufrnerksamkeit des Lesers auf diejenigen Ergebnisse 
der vorliegenden Abhandlung lenken, welche uns besonders beachtensweıt erscheinen. 
Für das atomare Elektron wissen wir, námlich, dessen Energie im n-ten Energiezustande 
durch die$Formel | (66), (72)] 


" m 373 -4 
"oem at sata 


auszudrücken. Stellt man diese der entsprechenden Dirac’ schen Formel, namlich, 


372 _ Y 
ve s ani 


gegenüber, so fallen uns die nuchstehenden Unterschiede frappant auf. Bei Dirac ist: 
En e s, bei uns aber: ee. Bei Dirao besteht die Möglichkeit eventuellen 
Verschwindens von /{(l+1)?—y77Z7}, falls | — 0, be: uns kann aber der entsprechende 
Ausdruck: y 1f14-1) + 429) nie verschwinden. Bei uns aind die Bahnen des atomaren 
Elektrons sämtlich radial gerichtet und von Ellıpsenbahnen ist keine Rede. Hierfür 
blicke man auf (74), (79) bzw. 74(a), 79(a) zurück. Darüber hinaus tritt uns das 
Eıgentumlichste der vorgelegten Theorie darın entgegen, dass sämtliche Quantenbahuen 
verdoppelt auftreten. Man stelle hierfür 74(a) bzw. 79(a) gegenuber (74) bzw. (79). Will 
man die Verdoppelungsverhältnisse bis auf die erste Naherung verfolgen, so hat man 
dazu (114) bzw. (116) an die Seite von (113).bzw. (115) zu stellen. Man sieht aber 
danach, dass bei nullter Naherung keine Doppelheit auftritt und die Bahn dann mit der 
klassischen zusammenfallt. 


Dieselbe‘ Verdoppelung von Quantenbahnen begegnet uns bei den Füllen (i) des 
Elektrons in konstantem elektrischen Felde und (it) des Elektrons in konstantem 
magnetischen Felde. Im ersten Falle stelle man (127), (146) einander gegenüber für die 
unverkürzten Quantenbahnen, aber (189), (152) für die Quantenbahnen bis auf die erste 
Nüherung. Im zweiten Falle vergleiche man (168), (194) miteinander für die unverkürzten 
Quantenbahnen, aber (179) mit (201) für die Quantenbahnen bis auf die erste Naherung. 
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In beiden Fällen verschwindet die Verdoppelung, falls man sich auf die nullte N&herung 
beschränkt und die Bahnen nehmen dann die klassische Gestalt an. Einzig beim 
feldfreien Elektron liefert die Theorie keine Verdoppelung von der Quantenbahn, welche 
auch mit der klassischen Bahn zusammenfällt. 


In der klassischen Hamiltcn-Jakobischen Dynamik erscheinen die Bahngleichungen 
des negativen Elektrons, wie man es merkt, durchweg mit einem negativen Vorzeichen 
der Eigenzeit dr versehen. Dieser Zug der klassischen Dynamik setzt sich durch die 
hier entworfene, eben auf der klassıschen Hamilton Jakobischen Dynamik gegründote 
Quantendynamik fort. Für die phonomenale Bahn, d.h. die von der phenomenalen 
Zeit t abhängige Bahn, um die nur der messende Physiker zu schaffen hat, hat es aber 
keinen Belang, da stets dt/dr auch negativ ist. Im Fulle des feldfreien Elektrons, zum 


l ‚Beispiel, sind die Gleichungen der phenomenalen Bahn naoh (210) oder (216) 


4 
` 


: dr dy Se 
: E adi = (p Do T). (218) 


Man überzeugt sich ohne allzu viele Mühe, dass in den entsprechenden klassischen 
Gleichungen des positiven Elektrons das Vorzeichen von dr positiv wird. Was sich am 
Tatbestand—das positive Vorzeichen von dr für das positive Elektron, das negative 


Vorzeichen hingegen für də negative Elektron—physikalisch äussern will, ist uns noch 
nicht klar. 


Nachtrag. Es leuchtet uns nachträglich ein, dass sich das Versagen der Kausalitat 
bei der oben ausgearbeiteten Quantendynamik als mit emer Unsicherheit der Anfangs- 
geschwindigkeit verbunden ansehen lässt. Man sıeht das am einfachsten an dıe oben 
behandelten Beispiele ein. Blickt man etwa beim Beispiele des konstanten elektrischen 
Feldes auf die Gleichungen (189) baw. (152), so merkt man, dass zu Anfang der 
quantendynamischen Laufbahn die  Gesehwindigkeitskomponenten seien entweder 
nach (139) 





dt - (a+! sı ərz Ro eH | 
gs Ze wd 
de _ e ch, w 
cdr Eo 2e, w^ | 
de 00.60.26 | 
cdr "ə cdr "e = ) 
d h (152 : 
oder nach (152) 21-ə 
dr Eo ch 28, W' 
dz C ? ch, w" 
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Ist die vorangehende Bahn klassisch erfolgt, so sind die Geschwindigkeitskomponenten 
an deren Einde nach (180) oder (150) 


dt _ eE ) ] 
Ze "IG. j 

de 00 " 

cdr e 

dz C d H 
ee sl 


Die Anfangsgeschwindigkeiten der quantendynamischen Bahn sind danach als unsicher 
zu meinen und die Unsicherheit betrifft nur die Komponenten dt/dr, da;cdr und im 


Masse von ə = bzw. puu, Eine Aussage gleicher Art gilt offenbar auch in dem ) 
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übrigen Fällen. 
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tx Genootschap” (The Mathematical Society of the Netherlands). It is the sincere hope of 

.” the ‘‘Wiskundig Genootschap”? that the Congress 1954, which will be open to all 
mathematicians from all parts of the world, will be a fertile international gathering. 

The Organizing Committee has invited a number of outstanding mathematicians to 
deliver one-hour addresses, hoping that in this way a survey or the recent development 
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d 3. Geometry and Topology. 4, Probability and Statistics. 
© D. Mathematical Physics and Applied Mathematics. 6. Logic and Foundations. 

n 7. Philosophy, History and Education. 


Those who wish to attend the Ccngress are requested to communicate their name 
(with degrees, qualifications etc.) and full address to the secretariat as soon as possible. 
They will receive a more detailed communication m due course. 
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